Sabera Bijani1,2, Faraz Shaikh1,3, Sheefa Mirza4,5, Shirley Weng In Siu3, Nayan Jain6, Rakesh Rawal4,6, Nigel G J Richards7, Anamik Shah1,8, Ashish Radadiya1,7. 1. Center of Excellence, National Facility for Drug Discovery Complex, Department of Chemistry, Saurashtra University, Rajkot 360005, India. 2. Department of Chemistry, Marwadi University, Rajkot 360003, India. 3. Department of Computer and Information Science, University of Macau, Macau 999078, China. 4. The Gujarat Cancer & Research Institute, Ahmedabad 380009, India. 5. Department of Internal Medicine, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg 2193, South Africa. 6. Department of Life Sciences, School of Sciences, Gujarat University, Ahmedabad 380009, India. 7. School of Chemistry, Cardiff University, Cardiff CF10 3AT, UK. 8. Astha, Saurashtra University Karmachari Cooperative Society, B/H Forensic Lab., Street No. 2, University Road, Rajkot 360005, India.
Abstract
P-glycoprotein (Pgp), an ATP binding cassette (ABC) transporter, is an ATP-dependent efflux pump responsible for cancer multidrug resistance. As part of efforts to identify human Pgp (hPgp) inhibitors, we prepared a series of novel triazole-conjugated dihydropyrimidinones using a synthetic approach that is well suited for obtaining compound libraries. Several of these dihydropyrimidinone derivatives modulate human P-glycoprotein (hPgp) activity with low micromolar EC50 values. Molecular docking studies suggest that these compounds bind to the M-site of the transporter.
P-glycoprotein (Pgp), an ATP binding cassette (ABC) transporter, is an ATP-dependent efflux pump responsible for cancer multidrug resistance. As part of efforts to identify human Pgp (hPgp) inhibitors, we prepared a series of novel triazole-conjugated dihydropyrimidinones using a synthetic approach that is well suited for obtaining compound libraries. Several of these dihydropyrimidinone derivatives modulate human P-glycoprotein (hPgp) activity with low micromolar EC50 values. Molecular docking studies suggest that these compounds bind to the M-site of the transporter.
Multidrug
resistance (MDR) remains a deadly obstacle in anticancer
treatments despite advances in understanding the molecular basis of
this phenomenon. A key contributor to MDR is P-glycoprotein (Pgp),
also known as multidrug resistance 1 protein (MDR1 gene), which is
a membrane protein that couples ATP hydrolysis to the expulsion of
small molecules, such as anticancer drugs, from the cell.[1] The resulting decrease in intracellular drug
concentration leads to reduced chemosensitivity and, ultimately, drug
resistance. The ability of human P-glycoprotein (hPgp) to adopt multiple
conformations and its substrate promiscuity, however, complicate efforts
to design potent inhibitors using rational, structure-based strategies.
This problem is reinforced by the fact that several binding sites
are located within the TM helices, which are believed to form a pathway
by which small molecules cross the membrane. To the best of our knowledge,
no hPgp inhibitors have been approved for clinical use in the last
four decades, including third generation candidates such as zosuquidar,
elacridar, and tariquidar.[2] New synthetic
strategies are therefore needed to explore chemical space so as to
obtain small molecules that can inhibit or at least modulate hPgp
activity to overcome MDR.In previous work in our laboratory,
we found that functionalized
1,4-dihydropyridines had activity as hPgp inhibitors,[3] which led to the discovery of compounds that combined anticancer
activity with the ability to inhibit hPgp.[4] Encouraged by this success, we are exploring whether alternate,
synthetically available, and heterocyclic scaffolds represent novel
leads for accessing hPgp inhibitors with improved potency. We now
report that the dihydropyrimidinone (DHPM) scaffold can also be functionalized
to give hPgp inhibitors by means of chemical transformations that
can be easily adapted to prepare diverse DHPM libraries.Our
work on the DHPM scaffold is motivated by literature showing
that DHPM derivatives exhibit activity against a variety of anticancer
targets. For example, DHPM derivative I (Figure ) inhibits the sodium iodide
symporter, which is a transmembrane glycoprotein associated with many
types of cancer.[5] Oxo-monastrol analogue II and its derivatives are cytotoxic,[6] and N-phenyl DHPM derivatives of III inhibit Hsp90 at nanomolar concentrations.[7] Lipophilic DHPM variants, such as IV, exhibit antiproliferative
properties against two glioma cell lines,[8] and derivatives of DHPM V are active against lung cancer
cell lines at nanomolar concentrations.[9] Finally, aryl-α-haloacrylamide-functionalized DHPMs, such
as VI, show anticancer activity in the micromolar range
by inhibiting tubulin polymerization.[10]
Figure 1
Structures
of known DHPMs [I–VI] with anticancer
activity (see text for details).
Structures
of known DHPMs [I–VI] with anticancer
activity (see text for details).
Results and Discussion
Using a convergent synthetic
sequence that combines the Biginelli
reaction with CuAAC chemistry, we prepared a library of novel triazole-conjugated
DHPM derivatives (TRZ-DHPMs) (Scheme ). 1,2,3-Triazoles are easily prepared by copper-catalyzed
alkyne–azide cycloaddition (CuAAC), thereby allowing the rapid
construction of compounds capable of exhibiting a range of biological
activities.[11,12] As well as being straightforward
to prepare, triazole substituents are present in clinically approved
drugs[13] and fungicides[14] as well as in ligands that bind to a variety of cancer
targets.[15,16] Substituted hydroxybenzaldehydes 1a–1d were propargylated in the presence of cesium carbonate to give benzaldehydes 2a–2d,[17] which
were then reacted with acetyl acetone and urea following the Biginelli
protocol[18] to yield DHPM derivatives 3a–3d. These alkynes were then connected
to the substituted phenylazides 5a′–5h′ (prepared from the corresponding anilines 4a′–4h′[19]) in the presence of CuSO4 and sodium ascorbate
to give the target 4-aryl-5-acetyl-6-methyl-3,4-dihydropyrimidin-2(1H)-ones 6 as racemic mixtures (Figure and Table S1, Supporting Information). This coupling reaction gave the
best yields when a 2:1:2 t-butanol:DMF:water mixture
was used as solvent (Table S2, Supporting
Information).
Scheme 1
Synthetic Route for Preparing the Library of Novel DPHM Derivatives
(a) Propargyl bromide and
DMF; (b) acetylacetone, urea, and PEG-400; (c) HCl, NaNO2, and NaN3; (d) 5a′–5h′, sodium ascorbate, CuSO4.5H2O, and (2:1:2) t-butanol:DMF:H2O.
Figure 2
Structures of selected TRZ-DPHMs for which IC50 <
1.5 μM against the Caco-2 cell line. Structures and biological
data for the complete set of TRZ-DPHMs are provided in the Supporting Information.
Structures of selected TRZ-DPHMs for which IC50 <
1.5 μM against the Caco-2 cell line. Structures and biological
data for the complete set of TRZ-DPHMs are provided in the Supporting Information.
Synthetic Route for Preparing the Library of Novel DPHM Derivatives
(a) Propargyl bromide and
DMF; (b) acetylacetone, urea, and PEG-400; (c) HCl, NaNO2, and NaN3; (d) 5a′–5h′, sodium ascorbate, CuSO4.5H2O, and (2:1:2) t-butanol:DMF:H2O.With
these compounds in hand, we characterized their cytotoxicity
against a colorectal adenocarcinoma (Caco-2) cell line using an MTT-based
assay. Half maximal inhibition concentration (IC50) values
in the range of 0.6–3.7 μM were obtained for the library
of functionalized TRZ-DHPMs (Table S3,
Supporting Information), which compares favorably to those observed
for known drugs such as carboplatin (7 μM), gemcitabine (5 μM)
and daunorubicin (14 μM) under the same assay conditions (Table ). The relatively
small range of IC50 values observed for this series of
functionalized TRZ-DHPMs complicates efforts to develop any SAR models
but the presence of a halogen in the N-phenyl substitute
of the triazole ring appears to correlate with sub-micromolar activity
(e.g., 6ad′, 6bd′, and 6be′). An ortho- rather than a para-relationship between the triazole and dihydropyrimidinone
rings also appears important for increased cytotoxicity (Table S3, Supporting Information). The molecular
basis for the cytotoxic activity of these TRZ-DHPMs, however, remains
to be elucidated.
Table 1
Activities of Selected TRZ-DPHMs and
Known Drugs
compound
cytotoxicity,IC50 (μM)
inhibition of calceinextrusion, EC50 (μM)
1
6aa′
1.1 ± 0.2
12.56 ± 0.08
2
6ab′
1.3 ± 0.8
15.42
± 0.06
3
6ad′
0.57 ± 0.03
3.14 ± 0.03
5
6bd′
0.76 ± 0.06
4.8 ± 0.2
6
6be′
0.8 ± 0.1
9.25 ± 0.02
7
6bh′
0.6
± 0.2
4.2 ± 0.2
8
6ca′
0.9 ± 0.1
4.0 ± 0.2
9
6cb′
1.22 ± 0.08
3.4
± 0.3
10
6cc′
1.25 ± 0.08
8.9 ± 0.1
11
carboplatin
7 ±
2
N.A.
12
gemcitabine
5 ± 1
N.A.
13
daunorubicin
14 ± 2
N.A.
14
verapamil
N.A.
9.52 ± 0.01
15
cyclosporin A
N.A.
5.36 ± 0.02
Despite a lack of information
on the molecular mechanism of action,
these nine cytotoxic TRZ-DHPMs were also assayed for their ability
to reverse the hPgp-mediated efflux of calcein-AM, a lipophilic, non-fluorescent
dye and a known Pgp substrate, from a Caco-2-VB cell line that expresses
hPgp in much higher amounts and is resistant to vinblastine (Figure ). In this assay,
hydrolysis of intracellular calcein-AM to calcein by endogenous cellular
esterases causes the cell to become highly fluorescent (see Figure S1 in Supporting Information). As a result,
compounds that block the function of the hPgp efflux pump, which otherwise
actively exports the dye from the cell, cause Caco-2-VB cells to accumulate
calcein-AM, thereby increasing their fluorescence. The nine TRZ-DHPMs
exhibiting the highest cytotoxicity were assayed for their ability
to block calcein-AM release from Caco-2-VB cells (Table ). The results of these experiments
show that the TRZ-DHPMs 6ad′, 6bd′, 6be′, 6bh′, 6ca′, 6cb′, and 6cc′ exhibit
similar or better efficacy than known Pgp substrates, cyclosporin
A and verapamil (Table and Figure ). The
most potent TRZ-DHPM, 6ad′, has a chlorophenyl
substituent attached to the triazole although replacing chlorine by
an electron-withdrawing nitro group gives a 5-fold decrease in potency
(6ab′). A similar decrease in potency is also
observed when chlorine is replaced by hydrogen, raising the possibility
that the ligand forms a σ-hole interaction with its target site
in the ATP-dependent efflux pump. Remarkably, the presence of bromine
substitution in the phenyl substituent of the dihydropyrimidi-none
(6cb′) overcomes the effects of the nitro group
and restores hPgp inhibition. Fluorescence microscopy measurements
confirmed intracellular calcein retention when Caco-2-VB cells were
treated with 1 μM calcein-AM in the presence of TRZ-DPHMs 6ad′, 6ca′, 6cb′, cyclosporin A, and verapamil (Figure ). Statistical analysis also showed that
that 6ad′ (p < 0.0001) is
a better modulator of hPgp function than either cyclosporin A or verapamil.
Thus, these compounds have dual activity although their cytotoxicity
is unlikely to be associated with their ability to modulate hPgp efflux
activity.
Figure 3
Top: Dose-dependent dye retention in Caco-2 VB cells treated with
selected TRZ-DHPMs in the calcein-AM uptake assay. Bottom: (A) Fluorescence
microscopy imaging of intracellular calcein in (a) untreated Caco-2
VB cells and Caco-2 VB cells 30 min after being treated with (b) TRZ-DHPM 6ad′, (c) verapamil, (d) cyclosporin A, (e) TRZ-DHPM 6cb′, and (f) TRZ-DHPM 6ca′. TRZ-DHPM
concentrations correspond to the relevant EC50 value (Table ). (B) Graphical representation
of calcein retention after treatment with TRZ-DHPMs 6ad′, 6cb′, and 6ca′. Data are
represented as mean ± SEM with p value (*p < 0.001) as compared to verapamil and cyclosporin A.
Calcein retention in Caco-2 cells was used as a negative control.
Full details are provided in the Supporting Information.
Top: Dose-dependent dye retention in Caco-2 VB cells treated with
selected TRZ-DHPMs in the calcein-AM uptake assay. Bottom: (A) Fluorescence
microscopy imaging of intracellular calcein in (a) untreated Caco-2
VB cells and Caco-2 VB cells 30 min after being treated with (b) TRZ-DHPM 6ad′, (c) verapamil, (d) cyclosporin A, (e) TRZ-DHPM 6cb′, and (f) TRZ-DHPM 6ca′. TRZ-DHPM
concentrations correspond to the relevant EC50 value (Table ). (B) Graphical representation
of calcein retention after treatment with TRZ-DHPMs 6ad′, 6cb′, and 6ca′. Data are
represented as mean ± SEM with p value (*p < 0.001) as compared to verapamil and cyclosporin A.
Calcein retention in Caco-2 cells was used as a negative control.
Full details are provided in the Supporting Information.As in a recent study of small-molecule
hPgp inhibitors,[20] molecular modeling was
used to obtain insights
into how these DHPM derivatives might inhibit this transporter. Given
the unavailability of high-resolution experimental coordinates for
the “inward”-facing conformation of hPgp, we built a
homology model of hPgp based on the experimental structure of Mus musculus Pgp, which has an 87% sequence match
to hPgp,[21] using well-precedented methods
(see the Supporting Information).[22] Energy minimization and equilibration using
molecular dynamics (MD) simulations (100 ns) gave an inward-facing
structure (Figure A) that is consistent with a low-resolution (3.6 Å) cryoEM structure
of hPgp (occluded conformation).[23] Molecular
docking was used to examine how TRZ-DHPM 6ad′ might
bind to the M-, R-, H-, and ATP-binding sites[24] of the efflux pump. Thus, the (R) and (S) enantiomers of TRZ-DHPM 6ad′ were
docked into the four potential binding sites of hPgp and evaluated
using the scoring function implemented in GLIDE.[25] We observed that (S)-6ad′ exhibited little energetic preference for any of the binding sites
(Table S4, Supporting Information), suggesting
that this enantiomer might easily be transported out of the cell by
hPgp. In contrast, the (R)- 6ad′ was predicted to bind more tightly to the M-site (−8.8 kcal
mol–1) than to the H- and R-sites (−3.7 and
– 6.5 mol–1). As neither enantiomer could
be docked into the ATP-binding site, our docking calculations lead
to the hypothesis that (R)-6ad′ blocks hPgp activity, a hypothesis that will be tested in future
experimental studies. In agreement with our hypothesis concerning
σ-hole binding interactions (see above), we found that the best-scoring
ligand pose placed the chloro substituent in (R)-6ad′ within 3 Å of the side chain oxygen in Ser344
(Figure B).
Figure 4
(A) Energy-minimized,
inward-facing conformation of the hPgp homology
model in a lipid bilayer, showing the homologous domains that exhibit
pseudo-2-fold symmetry. The N-terminal transmembrane (TMD) and nucleotide
binding (NBD) domains are colored red and green, respectively, and
the C-terminal TMD and NBD are shown in blue and orange, respectively.
M-, H-, and R-binding sites are rendered as magenta, cyan, and yellow
spheres. Head and tail groups of the lipid bilayer membrane are represented
by salmon spheres and pale cyan sticks, respectively. (B) Docking
poses of the (R)- and (S)-enantiomers
of 6ad′. (C) Protein–ligand interactions
computed for (R)-6ad′ bound at the M-site
of hPgp. (D) Protein–ligand interactions computed for (S)-6ad′ at the M-site of hPgp.
(A) Energy-minimized,
inward-facing conformation of the hPgp homology
model in a lipid bilayer, showing the homologous domains that exhibit
pseudo-2-fold symmetry. The N-terminal transmembrane (TMD) and nucleotide
binding (NBD) domains are colored red and green, respectively, and
the C-terminal TMD and NBD are shown in blue and orange, respectively.
M-, H-, and R-binding sites are rendered as magenta, cyan, and yellow
spheres. Head and tail groups of the lipid bilayer membrane are represented
by salmon spheres and pale cyan sticks, respectively. (B) Docking
poses of the (R)- and (S)-enantiomers
of 6ad′. (C) Protein–ligand interactions
computed for (R)-6ad′ bound at the M-site
of hPgp. (D) Protein–ligand interactions computed for (S)-6ad′ at the M-site of hPgp.Both enantiomers of 6ad′ exhibit
π–π
interactions with Phe343 and hydrogen-bond to Tyr310 when bound within
the M-binding site of the hPgp model (Figure C,D). The origin of the stereochemical preference
for M-site binding therefore appears to arise from the ability of
(R)-6ad′ to form an additional
hydrogen bond to Gln725, a conserved residue in the access tunnel
of the transporter (Figure C).[26]
Conclusions
These experimental and computational findings provide a firm basis
for the development of functionalized DHPMs as hPgp inhibitors. Our
synthetic route also provides a rapid entry into libraries that may
yield small molecules with improved or altered selectivity for clinically
relevant efflux pumps. Given the wide range of triazole substituents
that be tolerated, our route can also introduce fluorophores for visualizing
the binding of TRZ-DHPMs to hPgp and related transporters.[27]
Experimental Section
Materials
All the starting materials,
reagents, and catalysts were purchased from Aldrich, Merck, or Lobachem
and used without further purification. Anhydrous solvents (Spectrochem)
were stored over molecular sieves. Chromatographic solvents used for
the isolation/purification of compounds were distilled prior to use.
Thin layer chromatography (TLC) was performed using 0.2 mm precoated
plates of silica gel G60 F254 (Merck), and compounds were
visualized with either with UV light (254 nm) or iodine vapor.
Apparatus
1H (400 MHz)
and 13C (100 MHz) NMR spectra were recorded on a Bruker
Avance III instrument using DMSO-d6, CDCl3, and MeOD-d4 as solvents. 1H and 13C chemical shifts are reported in ppm relative
to tetramethylsilane (0.0). The following abbreviations designate
chemical shift multiplicities: s = singlet, bs = broad singlet, d
= doublet, dd = double doublet, t = triplet, dt = doublet of triplet,
m = multiplet. All 13C NMR spectra are proton-decoupled.
Melting points were determined in open capillary tubes and are uncorrected.
Mass spectra were recorded on a Shimadzu GC–MS-QP-2010 Ultra
spectrometer by direct injection. Elemental analyses were carried
out on a Euro Vector EURO EA 3000 model.
Synthetic
Procedures
2-(Prop-2-ynyloxy)benzaldehyde (2a)[28]
2-Hydroxybenzaldehyde (10.0
g, 81.9 mmol, 1 equiv) was charged into a round-bottom flask containing
DMF (30 mL) followed by the addition of Cs2CO3 (32.0 g, 98.2 mmol, 1.2 equiv), and the reaction was stirred at
rt for 5–10 min. Propargyl bromide (8.0 mL, 90.09 mmol, 1.1
equiv) was added dropwise into the reaction mixture, and the reaction
was stirred overnight at rt. After complete consumption of the starting
material, the reaction mass was poured into crushed ice. The resulting
precipitate obtained was collected by in vacuo filtration and washed
with water. Light yellow solid; 11.5 g 88.0%; 1H NMR (DMSO-d6): δ 9.45 (s, 1H, CHO), 7.68–7.70
(m, 2H, ArH), 7.18–7.23 (m, 2H, ArH), 4.25 (s, 2H, CH2), 3.36 (s, 1H, CH).The following compounds were prepared
using the same synthetic procedure as described for 2a.
Aldehyde 2a (5.0 g, 32.0 mmol), urea (1.92 g, 32.0 mmol), acetyl acetone (3.2
mL, 32.0 mmol), and PEG (4.0 g) were mixed and heated at 100 °C
for 3–5 h. After cooling to RT, the reaction mixture was poured
into water, and the solid product was removed by filtration before
being recrystalized from EtOH. White solid; 5.4 g, 59%; 1H NMR (DMSO-d6): δ 9.19 (s, 1H,
NH), 8.96 (s, 1H, NH), 7.52 (m, 2H, ArH), 7.20 (m, 2H, ArH) 5.24 (s,
3H, 1 × CH2 and 1 × CH), 2.29 (s, 3H, COCH3), 2.09 (s, 3H, CH3). 13C NMR (DMSO-d6): δ 192.3, 155.2, 151.0, 147.3, 128.8,
127.3, 126.9, 125.6, 114.6, 108.5, 78.8, 75.8, 53.1, 51.3, 30.2, 18.8.The following compounds were prepared using the same synthetic
procedure as described for 3a.
General Procedure for the
Preparation of Phenylazides
A solution of aq. NaNO2 (5 mmol) was added dropwise
to a solution of the aniline (5 mmol) dissolved in 6 N HCl (20 mL)
so as to maintain the temperature within a range of 0–5 °C.
After stirring at 0–5 °C for 1 h, a solution of aq. NaN3 (6 mmol, 1.2 equiv) was added dropwise to the reaction mixture.
The reaction was then allowed to warm and stirred at RT for 6–8
h with progress being monitored by TLC. On completion, the reaction
mixture was poured into cold water. The aqueous solution was extracted
with EtOAc and the combined extracts were then washed with brine and
aq. NaHCO3 before being dried over Na2SO4. Concentration of the solution under vacuum gave the desired
phenylazide.
General Procedure for Synthesis
of 5-Acetyl-6-methyl-Substituted
1,2,3-Triazol-4-yl)methoxy)phenyl)-3,4-dihydropyrimidin-2(1H)-one (6aa′–6dc′)
A mixture of substituted 5-acetyl-6-methyl-4-phenyl-3,4-dihydropyrimidin-2(1H)-one 3 (1 equiv) and phenylazide 5 (1.2 equiv) was charged into a sealed tube containing a solution
of tert-butanol, DMF, and water (2:1:2, v/v, 5 mL).
To this charged assembly, copper(II) sulfate (5 mol %) and sodium
ascorbate (1 equiv) were added. The reaction was continued under stirring
for 20–24 h at rt. After completion of the reaction (confirmed
by TLC), the reaction mixture was filtered under vacuum and washed
with water.Using the above general synthetic procedure, compounds 6aa′–6dc′ were synthesized.
The colorectal adenocarcinoma
Caco-2 cell line was generously provided by the Zydus Research Centre,
Ahmedabad, Gujarat, India. These cells were maintained in DMEM (Hi-Media)
containing 25 mM HEPES buffer, 1000 mg/L glucose, l-glutamine,
sodium bicarbonate, and sodium pyruvate and supplemented with 20%
heat-inactivated FBS (cell clone), 1% (v/v) minimal essential medium
nonessential amino acids (Hi-Media), antibiotic, and antimycotic (penicillin,
streptomycin, and amphotericin B; Hi-Media) in a humidified atmosphere
of 5% CO2 at 37 °C. Cells were rinsed and split with
trypsin–EDTA at 70–80% confluency, and only those cells
exhibiting more than 90% viability were used in the biological assays.
Caco-2 VB cells, which overexpress hPGP, were generated by growing
Caco-2 cells in the presence of 1 μg/mL vinblastine. The drug
was removed from the culture medium 2–3 days prior to using
the Caco-2 VB cells in the biological assays.
MTT Assay
Caco-2 cells were cultured
in 96-well plates until they reached 70–80% confluency and
then treated with the TRZ-DHPM derivatives (0.1–100 μM).
Cytotoxicity was determined using an MTT (3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium
bromide) assay as described previously.[31] Stock solutions of the TRZ-DHPM derivatives contained DMSO, but
the final DMSO concentration in each assay did not exceed 0.1%.
hPgp-Mediated Efflux Assay
hPgp-mediated
efflux in Caco-2 VB cells was measured using the Vybrant Multidrug
Resistance Assay Kit (V-13180, Thermo Fisher Scientific). Each well
contained 100 μL of cells in culture medium, 50 μL of
a solution containing the TRZ-DPHM derivative, and 50 μL of
calcein-AM (total volume, 200 μL). The final concentration of
putative hPGP inhibitor in each assay ranged from 0.01 to 30 μg/mL.
Verapamil and cyclosporin A, both hPgp inhibitors, were used over
the same concentration range as that used for the TRZ-DPHM derivatives.
Caco-2 cells were used as a negative control (Figure S1). Cells were incubated at 37 °C for 15 min
after the addition of TRZ-DPHM derivatives before 1 μM calcein-AM
was added to each well (0.25 μM final concentration). After
an additional 15 min at 37 °C, cells were washed with cold tissue
culture medium (200 μL), and calcein retention was measured
by the increase in intracellular fluorescence. PBS was used for all
dilutions. All experiments were performed in triplicate.
Cell Viability Assay
Cell viability
was checked using a Trypan blue dye exclusion assay as described previously.[31] In this assay, the % of viable cells was determined
by Automated Cell Counter TC10 (BioRad) based on their capacity to
uptake dye that differentiates live and dead cells. Live cells exclude
dye due to their intact membrane, and dye penetrates the cell membrane
of dead cells.
Data Analysis
The effectiveness
of inhibitors on Pgp efflux activity, i.e., calcein retention was
calculated by the following equationwhere FT is the fluorescence of treated cells
and FU is the fluorescence of untreated cells.
Computational Methods
Homology Modeling
A hPGp model
was built using the multiple-template approach implemented in Modeller
9v14.[32] Experimental coordinates 3G60 and
3G5U from M. musculus and 4F4C from C. elegans having high sequence similarity with human
Pgp were used as templates for homology modeling.[33] The preparation of the hPgp model involved hydrogen addition
and the removal of water molecules, followed by minimization using
the Impact module of Schrödinger. The quality of the hPGp model
was assessed using the pdbsum server.[34] A total of 96.4% residues were found to be in the allowed region
with 86.7% in the most favored region.
Molecular
Dynamics Simulation
Molecular
dynamics (MD) simulations were performed on the final hPgp model using
the Desmond software package. Lipid bilayer position coordinates were
obtained from the Orientations of Proteins in Membranes (OPM) database.[35] The Desmond membrane insertion tab was then
used to place the transmembrane (TM) coordinates. Further solvation
of the model was done with TIP4P model in the cubic periodic boundary
box. The resulting system of 76,326 atoms was minimized and pre-equilibrated
by standard procedures implemented in Desmond. The distance between
the box wall and hPgp model was set to be greater than 10 Å,
thereby avoiding direct interactions with the periodic images. Steepest
descent energy minimization was carried out until the gradient threshold
(25 kcal/mol/Å), followed by L-BFGS (low-memory Broyden–Fletcher–Goldfarb–Shanno
quasi-Newtonian minimizer) until a convergence threshold of 1 kcal/mol/Å
was met. Further MD simulations were carried on the equilibrated systems
for a desired period at a constant temperature of 300 K and constant
pressure for the duration of 100 ns. The equations of motion were
integrated with a 2 fs time step in the NPT ensemble. The SHAKE algorithm
was applied to all hydrogen atoms, the van der Waals cutoff was 9
Å, and the temperature was maintained at 300 K, employing a Nose–Hoover
thermostat with a relaxation time of 1 ps.
Molecular
Docking
The energy-minimized
hPgp model was used in the docking study. R and S enantiomers of TRZ-DHPMs
were docked to the three drug-binding sites categorized by Ferreira
et al.[24] as well as the ATP-binding site.[36] Compound structures were prepared with the Ligprep
module of Schrödinger. All three binding site grids were calculated
using the Glide module of Schrödinger,[25] and then the Glide XP docking was performed in all four binding
sites individually.
Authors: Tamara G M Treptow; Fabrício Figueiró; Elisa H F Jandrey; Ana M O Battastini; Christianne G Salbego; Juliana B Hoppe; Priscila S Taborda; Sabrina B Rosa; Luciana A Piovesan; Caroline Da R Montes D'Oca; Dennis Russowsky; Marcelo G Montes D'Oca Journal: Eur J Med Chem Date: 2015-03-28 Impact factor: 6.514
Authors: Daniel L da Silva; Fabiano S Reis; Dandara R Muniz; Ana Lúcia T G Ruiz; João E de Carvalho; Adão A Sabino; Luzia V Modolo; Angelo de Fátima Journal: Bioorg Med Chem Date: 2012-02-22 Impact factor: 3.641
Authors: Stephen G Aller; Jodie Yu; Andrew Ward; Yue Weng; Srinivas Chittaboina; Rupeng Zhuo; Patina M Harrell; Yenphuong T Trinh; Qinghai Zhang; Ina L Urbatsch; Geoffrey Chang Journal: Science Date: 2009-03-27 Impact factor: 47.728
Authors: Mikhail A Lomize; Irina D Pogozheva; Hyeon Joo; Henry I Mosberg; Andrei L Lomize Journal: Nucleic Acids Res Date: 2011-09-02 Impact factor: 16.971