Grid-independent Descriptors (GRIND) Analysis and SAR Guided Molecular Docking Studies to Probe Selectivity Profiles of Inhibitors of Multidrug Resistance Transporters ABCB1 and ABCG2.

BACKGROUND: ATP-binding cassette (ABC) transporters, P-glycoprotein (P-gp, ABCB1) and breast cancer resistance protein (BCRP/ABCG2) are major determinants of pharmacokinetic, safety and efficacy profiles of drugs thereby effluxing a broad range of endogenous substances across the plasma membrane. Overexpression of these transporters in various tumors is also implicated in the development of multidrug resistance (MDR) and thus, hampers the success of cancer chemotherapy. Modulators of these efflux transporters in combination with chemotherapeutics could be a promising concept to increase the effective intracellular concentration of anticancer drugs. However, broad and overlapped specificity for substrates and modulators of ABCB1 and ABCG2, merely induce toxicity and unwanted drug-drug interactions and thus, lead to late-stage failure of drugs. OBJECTIVE: In present investigation, we aim to identify specific 3D structural requirements for selective inhibition of ABCB1 and ABCG2 transport function. METHOD: GRID Independent Molecular Descriptor (GRIND) models of selective inhibitors of both transporters have been developed, using their most probable binding conformations obtained from molecular docking protocol. RESULTS: Our results demonstrated a dominant role of molecular shape and different H-bonding patterns in drug-ABCB1/ABCG2 selective interactions. Moreover, distinct distances of different pharmacophoric features from steric hot spots of the molecules provided a strong basis of selectivity for both transporters. Additionally, our results suggested the presence of two H-bond donors at a distance of 8.4-8.8 Å in selective modulators of ABCG2. CONCLUSION: Our findings concluded that molecular shape along with three dimensional pattern of Hbonding in MDR modulators play a critical role in determining the selectivity between the two targets.

Structural basis of drugs that increase cardiac inward rectifier Kir2.1 currents.

AIMS: We hypothesize that some drugs, besides flecainide, increase the inward rectifier current (IK1) generated by Kir2.1 homotetramers (IKir2.1) and thus, exhibit pro- and/or antiarrhythmic effects particularly at the ventricular level. To test this hypothesis, we analysed the effects of propafenone, atenolol, dronedarone, and timolol on Kir2.x channels. METHODS AND RESULTS: Currents were recorded with the patch-clamp technique using whole-cell, inside-out, and cell-attached configurations. Propafenone (0.1 nM-1 µM) did not modify either IK1 recorded in human right atrial myocytes or the current generated by homo- or heterotetramers of Kir2.2 and 2.3 channels recorded in transiently transfected Chinese hamster ovary cells. On the other hand, propafenone increased IKir2.1 (EC50 = 12.0 ± 3.0 nM) as a consequence of its interaction with Cys311, an effect which decreased inward rectification of the current. Propafenone significantly increased mean open time and opening frequency at all the voltages tested, resulting in a significant increase of the mean open probability of the channel. Timolol, which interacted with Cys311, was also able to increase IKir2.1. On the contrary, neither atenolol nor dronedarone modified IKir2.1. Molecular modelling of the Kir2.1-drugs interaction allowed identification of the pharmacophore of drugs that increase IKir2.1. CONCLUSIONS: Kir2.1 channels exhibit a binding site determined by Cys311 that is responsible for drug-induced IKir2.1 increase. Drug binding decreases channel affinity for polyamines and current rectification, and can be a mechanism of drug-induced pro- and antiarrhythmic effects not considered until now.

Structure-activity relationships and in silico models of P-glycoprotein (ABCB1) inhibitors.

1. The efflux pump p-glycoprotein (P-gp/ABCB1) has received enormous attention in drug (xenobiotic) disposition due to its role in modulation of the drug availability and in protection of sensitive organs. 2. P-gp mediated efflux is one of main mechanisms for multidrug resistance in cancer cells. A main approach to reverse the resistance and restore the drug efficacy is to use specific inhibitors of P-gp that suppress the efflux activity. 3. This review summarizes the binding capabilities of known chemical inhibitors based on the analyses of structure-activity relationships, and computational modeling of the inhibitors as well as the binding site of P-gp protein. 4. The molecular models will facilitate the design of lead inhibitors as drug candidates. Also, it helps scientists in early drug discovery phase to synthesize chemical series with better understanding of their P-gp binding liabilities.

Phenylpropiophenone derivatives as potential anticancer agents: synthesis, biological evaluation and quantitative structure-activity relationship study.

Series of twelve chalcone and propafenone derivatives has been synthesized and evaluated for anticancer activities against HeLa, Fem-X, PC-3, MCF-7, LS174 and K562 cell lines. The 2D-QSAR and 3D-QSAR studies were performed for all compounds with cytotoxic activities against each cancer cell line. Partial least squares (PLS) regression has been applied for selection of the most relevant molecular descriptors and QSAR models building. Predictive potentials of the created 2D-QSAR and 3D-QSAR models for each cell line were compared, by use of leave-one-out cross-validation and external validation, and optimal QSAR models for each cancer cell line were selected. The QSAR studies have selected the most significant molecular descriptors and pharmacophores of the chalcone and propafenone derivatives and proposed structures of novel chalcone and propafenone derivatives with enhanced anticancer activity on the HeLa, Fem-X, PC-3, MCF-7, LS174 and K562 cells.

Structure-activity relationships, ligand efficiency, and lipophilic efficiency profiles of benzophenone-type inhibitors of the multidrug transporter P-glycoprotein.

The drug efflux pump P-glycoprotein (P-gp) has been shown to promote multidrug resistance (MDR) in tumors as well as to influence ADME properties of drug candidates. Here we synthesized and tested a series of benzophenone derivatives structurally analogous to propafenone-type inhibitors of P-gp. Some of the compounds showed ligand efficiency and lipophilic efficiency (LipE) values in the range of compounds which entered clinical trials as MDR modulators. Interestingly, although lipophilicity plays a dominant role for P-gp inhibitors, all compounds investigated showed LipE values below the threshold for promising drug candidates. Docking studies of selected analogues into a homology model of P-glycoprotein suggest that benzophenones show an interaction pattern similar to that previously identified for propafenone-type inhibitors.

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

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

Trapping and dissociation of propafenone derivatives in HERG channels.

BACKGROUND AND PURPOSE: Human ether-a-go-go related gene (HERG) channel inhibitors may be subdivided into compounds that are trapped in the closed channel conformation and others that dissociate at rest. The structural peculiarities promoting resting state dissociation from HERG channels are currently unknown. A small molecule-like propafenone is efficiently trapped in the closed HERG channel conformation. The aim of this study was to identify structural moieties that would promote dissociation of propafenone derivatives. EXPERIMENTAL APPROACH: Human ether-a-go-go related gene channels were heterologously expressed in Xenopus oocytes and potassium currents were recorded using the two-microelectrode voltage clamp technique. Recovery from block by 10 propafenone derivatives with variable side chains, but a conserved putative pharmacophore, was analysed. KEY RESULTS: We have identified structural determinants of propafenone derivatives that enable drug dissociation from the closed channel state. Propafenone and four derivatives with 'short' side chains were trapped in the closed channel. Five out of six bulky derivatives efficiently dissociated from the channel at rest. One propafenone derivative with a similar bulk but lacking an H-bond acceptor in this region was trapped. Correlations were observed between molecular weight and onset of channel block as well as between pK(a) and recovery at rest. CONCLUSION AND IMPLICATIONS: The data show that extending the size of a trapped HERG blocker-like propafenone by adding a bulky side chain may impede channel closure and thereby facilitate drug dissociation at rest. The presence of an H-bond acceptor in the bulky side chain is, however, essential.

Stereoselective glucuronidation of propafenone and its analogues by human recombinant UGT1A9.

Stereoselective glucuronidation of propafenone and its beta-blocker analogues by human recombinant UGT1A3 and UGT1A9 from the recombinant baculovirus in insect sf9 cells was studied. The glucuronides produced in incubation mixtures were assayed by HPLC equipped with UV detector, and identified by beta-glucuronidase. The stereoselective glucuronidation was measured by pre-column 2,3,4,6-tetra-O-acetyl-beta-D-glucopyranosyl isothiocynate (GITC) derivatization HPLC method for propafenone and esomolol. In all of ten beta-blocker drugs studied, six showed the glucuronidation activity with UGT1A9, while four with UGT1A3. From roughly quantitative stereoselective glucuronidation study of racemic beta-blocker analogues by UGT1A9, propranolol had a high ratio of the ratios of S- to R-isomer glucuronide (S-G/R-G), about 4.3, the ratios of terbutaline, atenolol and esomolol were 3.3, 3.1 and 2.8 respectively, sotalol and propafenone were 2.3 and 2.0. In a word, S-isomers of these drugs were glucuronidated by human UGT1A9 much faster than their antipodes.

Mechanisms of melanogenesis inhibition by propafenone.

Melanogenesis is a physiological process that results in the synthesis of melanin pigments, which play a crucial protective role against skin photocarcinogenesis. The present study was conducted to determine the inhibitory effects of propafenone on melanogenesis and to elucidate the molecular events involved in the inhibition of melanogenesis by propafenone. To accomplish this, several experiments were conducted using human epidermal melanocyte cells. The melanin content and cAMP production were evaluated, and western blots for proteins involved in melanogenesis were conducted. The melanin content was significantly inhibited by propafenone in a concentration-dependent manner. To clarify the mechanism of the depigmenting property of propafenone, we examined the involvement of propafenone in cAMP signaling. In the cAMP production assay, the intracellular cAMP level was reduced by propafenone. The level of microphthalmia-associated transcription factor (MITF) protein, the upstream transcription factor of tyrosinase, was also reduced by propafenone. In addition, propafenone inhibited the expression of tyrosinase, TRP-1, and TRP-2. Taken together, the results of our study show that propafenone inhibits melanogenesis by suppressing cAMP production, which is involved in the expression of melanogenesis-related proteins and suggests that propafenone may be an effective inhibitor of hyperpigmentation.

The hERG potassium channel and drug trapping: insight from docking studies with propafenone derivatives.

The inner cavity of the hERG potassium ion channel can accommodate large, structurally diverse compounds that can be trapped in the channel by closure of the activation gate. A small set of propafenone derivatives was synthesized, and both use-dependency and recovery from block were tested in order to gain insight into the behavior of these compounds with respect to trapping and non-trapping. Ligand-protein docking into homology models of the closed and open state of the hERG channel provides the first evidence for the molecular basis of drug trapping.

Nanostructure-initiator mass spectrometry: a protocol for preparing and applying NIMS surfaces for high-sensitivity mass analysis.

Nanostructure-initiator mass spectrometry (NIMS) is a new surface-based MS technique that uses a nanostructured surface to trap liquid ('initiator') compounds. Analyte materials adsorbed onto this 'clathrate' surface are subsequently released by laser irradiation for mass analysis. In this protocol, we describe the preparation of NIMS surfaces capable of producing low background and high-sensitivity mass spectrometric measurement using the initiator compound BisF17. Examples of analytes that adsorb to this surface are small molecules, drugs, lipids, carbohydrates and peptides. Typically, NIMS is used to analyze samples ranging from simple analytical standards and proteolytic digests to more complex samples such as tissues, cells and biofluids. Critical experimental considerations of NIMS are described. Specifically, NIMS sensitivity is examined as a function of pre-etch cleaning treatment, etching current density, etching time, initiator composition, sample concentration, sample deposition method and laser fluence. Typically, NIMS surface preparation can be completed in less than 2 h. Subsequent sample preparation requires 1-5 min, depending on sample deposition method. Mass spectrometric data acquisition typically takes 1-30 s per sample.

Interaction field based and hologram based QSAR analysis of propafenone-type modulators of multidrug resistance.

Overexpression of membrane bound, ATP-dependent transport proteins is one of the predominant mechanisms leading to multiple drug resistance in tumor therapy as well as in the treatment of bacterial and fungal infections. In tumor therapy, P-glycoprotein (P-gp, ABCB1) is responsible for transport of a wide variety of natural product toxins out of tumor cells leading to decreased accumulation of cytotoxic drugs within the cells. Inhibition of P-gp thus gives rise to a resensitization of multidrug resistant tumor cells and represents a versatile approach for modulation of multidrug resistance. Within this paper, a set of propafenone-type inhibitors of P-gp were analyzed using both interaction field based methods such as CoMFA and CoMSIA and Hologram QSAR. With both methods, highly predictive models with q2-values>0.65 were obtained. Models using logP as additional descriptor generally yielded higher predictive power. On basis of unfavorable steric and favorable electrostatic and hydrophobic interaction fields, these models were able to explain all outlayers identified in previous Hansch-analyses. For HQSAR analysis, models with q2-values up to 0.72 were obtained. Positive influences were found for electron donating groups on the aromatic systems. Highly negative influences were found for diphenylalkylamine substituents, which is a further hint for steric hindrance. The models with highest predictive power were used for screening of a small virtual library. Synthesis and pharmacological testing of a sub set of this library showed that the external predictivity of the HQSAR models generally is lower than the internal one.

P-glycoprotein substrate binding domains are located at the transmembrane domain/transmembrane domain interfaces: a combined photoaffinity labeling-protein homology modeling approach.

P-glycoprotein (P-gp) is an energy-dependent multidrug efflux pump conferring resistance to cancer chemotherapy. Characterization of the mechanism of drug transport at a molecular level represents an important prerequisite for the design of pump inhibitors, which resensitize cancer cells to standard chemotherapy. In addition, P-glycoprotein plays an important role for early absorption, distribution, metabolism, excretion, and toxicity profiling in drug development. A set of propafenonetype substrate photoaffinity ligands has been used in this study in conjunction with matrix-assisted laser desorption/ionization time-of-flight mass spectrometry to define the substrate binding domain(s) of P-gp in more detail. The highest labeling was observed in transmembrane segments 3, 5, 8, and 11. A homology model for P-gp was generated on the basis of the dimeric crystal structure of Vibrio cholerae MsbA, an essential lipid transporter. Thereafter, the labeling pattern was projected onto the 3D atomic-detail model of P-gp to allow a visualization of the binding domain(s). Labeling is predicted by the model to occur at the two transmembrane domain/transmembrane domain interfaces formed between the amino- and carboxyl-terminal half of P-gp. These interfaces are formed by transmembrane (TM) segments 3 and 11 on one hand and TM segments 5 and 8 on the other hand. Available data on LmrA and AcrB, two bacterial multidrug efflux pumps, suggest that binding at domain interfaces may be a general feature of polyspecific drug efflux pumps.

Intramolecular distribution of hydrophobicity influences pharmacological activity of propafenone-type MDR modulators.

Lipophilicity is one of the major determining physicochemical descriptors for P-glycoprotein (P-gp) inhibitory activity. Recently, Pajeva and Wiese showed that in case of P-gp interaction, lipophilicity may be regarded as space-directed property. In the present study, a series of propafenone-type P-gp inhibitors with systematically varying hydrophobicity distribution within the molecules were synthesised and pharmacologically tested. QSAR studies on the basis of multiple linear regression analysis showed that with increasing lipophilicity of the substituents on the amine moiety, the statistical significance of the indicator variables, denoting the substitution pattern on the central aromatic ring system, also increases. This indicates that the distribution of hydrophobicity within the molecules influences the mode of interaction with P-gp.

MCASE study of the multidrug resistance reversal activity of propafenone analogs.

A database containing 130 propafenone type chemicals which have been tested for their multidrug resistance (MDR) reversal activity was compiled. Using the Multiple Computer-Automated Structure Evaluation (MCASE) program to analyze this database, an underlying relationship between MDR reversal activity and octanol/water partition coefficient was found. An MDR reversal model was created based on this database by the baseline activity identification algorithm (BAIA) of the MCASE program. The main phamacophores relevant to MDR reversal activity were identified.

Three-dimensional quantitative structure-activity relationship analysis of propafenone-type multidrug resistance modulators: influence of variable selection on test set predictivity.

An extended set of multidrug-resistance modulators of the propafenone type were investigated using CoMFA and CoMSIA. A number of 3D-QSAR models were derived from steric, electrostatic, and hydrophobic fields and their combinations. The hydrophobic fields alone and in combination with the steric and both (steric and electrostatic) fields yielded the models with the highest cross-validated predictivity, in agreement with a previous analysis of a smaller data set of propafenone-type multidrug-resistance (MDR) modulators. Inclusion of lipophilicity did not lead to an improvement of the models. The results point to the importance of hydrophobicity as a space-directed molecular property for MDR-modulating activity. The influence of variable selection applying the GOLPE procedure was investigated with an external test set. Variable-selection procedure was repetitively applied, keeping at each stage variables with uncertain contribution to the models. For the CoMFA-based 3D-QSAR models, an increase in external prediction quality was found. In contrast, the CoMSIA-based 3D-QSAR models were not improved by the GOLPE variable-selection procedure.

Effects of a series of dihydroanthracene derivatives on drug efflux in multidrug resistant cancer cells.

A set of 9,10-dihydro-9,10-ethano and ethenoanthracene derivatives was tested with the aim to quantify the effect observed on drug efflux. Structure activity relationships and molecular modeling studies allowed to define topological display of pharmacophoric groups for these reversal agents.

New multidrug resistance reversal agents.

Multidrug resistance (MDR) to antitumor agents represents a significant challenge to effective chemotherapy. The use of MDR modulators is a promising approach to overcome the undesired MDR phenotype. The more effective MDR modulators are urgently needed for clinical use. This review focuses on literatures published in 1998-2001.

A combined Hansch/Free-Wilson approach as predictive tool in QSAR studies on propafenone-type modulators of multidrug resistance.

A series of 48 propafenone-type modulators of multidrug resistance was synthesized and their P-glycoprotein inhibitory activity was measured using the daunomycin efflux assay. Both a Free-Wilson and a combined Hansch/Free-Wilson analysis were performed using log P, partial log P and molar refraction values as Hansch descriptors. The results of the Free-Wilson analysis show that modifications on the central aromatic ring generally influence pharmacological activity, whereby in almost all cases a decrease in MDR-modulating potency is observed (Q2cv = 0.66). The combined approach results in equations with remarkably higher predictive power (Q2cv = 0.83), specifically molar refractivity shows high significance in all equations derived. This indicates that polar interactions also contribute to protein binding.

Synthesis and pharmacological activity of the stereoisomers of GP-88, a propafenone-type modulator of multidrug resistance.

All four stereoisomers of the propafenone-type MDR-modulator GP-88 (1) were synthesized using a combined approach with chiral pool building blocks and an acetalic protective group, which allows not only diastereoseparation but also assignment of absolute configuration via NMR spectroscopy. Those isomers with different configuration on the center of chirality in the propanolamine side chain showed statistically different PGP-inhibitory activity. Generally, the (R)-configured isomers were by a factor of nearby two higher active than the (S)-isomers. No differences in activity were observed for isomers with different configuration on the benzylic center of chirality.