- Open Access
Molecular model of the outward facing state of the human P-glycoprotein (ABCB1), and comparison to a model of the human MRP5 (ABCC5)
Theoretical Biology and Medical Modellingvolume 4, Article number: 33 (2007)
Multidrug resistance is a particular limitation to cancer chemotherapy, antibiotic treatment and HIV medication. The ABC (ATP binding cassette) transporters human P-glycoprotein (ABCB1) and the human MRP5 (ABCC5) are involved in multidrug resistance.
In order to elucidate structural and molecular concepts of multidrug resistance, we have constructed a molecular model of the ATP-bound outward facing conformation of the human multidrug resistance protein ABCB1 using the Sav1866 crystal structure as a template, and compared the ABCB1 model with a previous ABCC5 model. The electrostatic potential surface (EPS) of the ABCB1 substrate translocation chamber, which transports cationic amphiphilic and lipophilic substrates, was neutral with negative and weakly positive areas. In contrast, EPS of the ABCC5 substrate translocation chamber, which transports organic anions, was generally positive. Positive-negative ratios of amino acids in the TMDs of ABCB1 and ABCC5 were also analyzed, and the positive-negative ratio of charged amino acids was higher in the ABCC5 TMDs than in the ABCB1 TMDs. In the ABCB1 model residues Leu65 (transmembrane helix 1 (TMH1)), Ile306 (TMH5), Ile340 (TMH6) and Phe343 (TMH6) may form a binding site, and this is in accordance with previous site directed mutagenesis studies.
The Sav1866 X-ray structure may serve as a suitable template for the ABCB1 model, as it did with ABCC5. The EPS in the substrate translocation chambers and the positive-negative ratio of charged amino acids were in accordance with the transport of cationic amphiphilic and lipophilic substrates by ABCB1, and the transport of organic anions by ABCC5.
The transport of small organic molecules and ions across cell membranes generally requires a transporter protein, and these transporter proteins have recognition sites that make them specific for particular substrates. Drugs can interact with these recognition sites and inhibit the transporter, or be substrates themselves. There is an increasing focus on transporters as drug targets, and the information on transporter structure and function is rapidly increasing. The number of drugs interacting with transporters will probably increase in the future.
According to the transporter classification approved by the transporter nomenclature panel of the International Union of Biochemistry and Molecular Biology [1, 2], transporters are divided into classes based on both function and phylogeny. These classes are: 1. Channels and pores, 2. Electrochemical potential-driven transporters (secondary transporters), 3. Primary active transporters, 4. Group translocators, 5. Transport electron carriers, 8. Accessory factors involved in transport, 9. Incompletely characterized transport systems.
ABC (ATP binding cassette) transporters belong to class 3 (primary active transporters), subclass A (diphosphate bond hydrolysis-driven transporters) and family 1 (ABC superfamily) . Primary active transporters use a primary source of energy to drive active transport of particles from regions of low concentration to regions of high concentration. The ABC superfamily transporters are structurally related membrane proteins sharing a common intracellular motif that exhibits ATPase activity that cleaves ATP's terminal phosphate, using the free energy from ATP (adenosine triphosphate) stored in the high-energy phosphate bond as the energy source for activating the transporter [1–4].
The human genome encodes more than 40 ABC transporters divided into five different subfamilies: ABCA, ABCB, ABCC, ABCD and ABCG, based on phylogenetic analysis (Additional file 1). According to the TCDB , these subfamilies belong to subclasses 3.A.1.201–212, ABC-type efflux permeases (mostly eukaryotic) . The ABC genes are highly conserved between species, indicating that most of these genes have been present since the beginning of eukaryotic evolution . These transporters feature both transmembrane domains (TMD) and nucleotide binding domains (NBD). In general, the domain arrangement of these transporters is TMD-NBD-TMD-NBD, but TMD0-TMD-NBD-TMD-NBD, NBD-TMD-NBD-TMD, TMD-NBD and NBD-TMD also exist [5, 6]. TMD0 is a 5 TMH amino-terminal domain present in ABCC1, ABCC2 and ABCC3. The NBD contains the Walker A and B motifs  and a signature C motif. Two further subfamilies, ABCE and ABCF, are related to ABC transporters, but they lack transmembrane domains and thus are not membrane transporters [4, 5]. The substrate specificity is provided by the TMDs, which contain 6–11 transmembrane helices (TMHs) .
Cells exposed to toxic compounds can develop resistance by a number of mechanisms, including increased excretion. The result is multidrug resistance (MDR), which is a particular limitation to cancer chemotherapy, antibiotic treatment and HIV medication. Transporters in subfamilies ABCA, ABCB, ABCC and ABCG are involved in multidrug resistance [8–11]. Development of inhibitors of drug efflux transporters has been sought for use as supplement to therapy to overcome multidrug resistance . In order to elucidate structural and molecular concepts of multidrug resistance, we have focused on the TMDs of ABCB1 and ABCC5 using molecular modeling techniques. ABCB1 and ABCC5 both have a TMD-NBD-TMD-NBD arrangement, with TMDs consisting of 6 TMHs. ABCB1 transports cationic amphiphilic and lipophilic substrates [13–16], while ABCC5 transports organic anions [17, 18]. Information about the molecular aspects of ligand interactions with these transporters can be used to design therapeutic agents that may aid to overcome multidrug resistance.
Several electron density maps of ABCB1 have been published [19–22], giving insight into ABCB1 architecture. The latest electron density map had a resolution limit of ~8 Å, and although this structure reveals the TMH packing of ABCB1, it is not possible, at this resolution, to predict the TMH numbering . In lack of an X-ray crystal structure, molecular modeling by homology may be an alternative for gaining structural insight into protein drug targets. The bacterial ABC transporter Sav1866 from Staphylococcus aureus has been crystallized in an outward-facing ATP-bound state . Sav1866 is a bacterial homologue to ABCB1 , indicating that the Sav1866 crystal structure could be used as a template for the present model building by homology. The 12 TMH arrangement of the Sav1866 crystal structure is consistent with the electron density maps of ABCB1 . The NBDs of both Sav1866 and ABCB1 are functionally equivalent; both NBDs are responsible for ATP binding and hydrolysis [23, 24].
In this study we have constructed an ABCB1 model based on the Sav1866 crystal structure  using molecular modeling techniques. Among the transporters in the ABCC subfamily (multidrug resistance proteins, MRPs), ABCC5 has a "P-gp-like" ("ABCB1-like") core domain organization (TMD1-NBD1-TMD2-NDB2) . We have previously constructed an ABCC5 model in a cGMP docking study (submitted) using the Sav1866 crystal structure  as a template, and in the present study we have performed a comparative analysis of the ABCB1 and ABCC5 models in order to understand the molecular concepts of the substrate difference between ABCB1 and ABCC5. The comparative analysis included the electrostatic potential surfaces (EPS) of the substrate translocation chambers, and the positive-negative ratios of charged amino acids of the TMDs of both models. A phylogenetic analysis of human ABC transporters has been performed in order to understand the phylogenetic relationship between ABCB1 and ABCC5. The ABCB1 model has been compared with cross-linking and site directed mutagenesis data published on ABCB1 [26–32].
Evolutionary tree of the human ABC transporters
The evolutionary tree of the human ABC transporters, together with Sav1866, is shown in Figure 1. ABCB1 and Sav1866 are localized on the same branch of the evolutionary tree (the "ABCB-branch"), while ABCC5 is localized on a different branch (the "ABCC-branch").
Amino acid sequence identities of TMDs and positive-negative ratios of charged amino acids
Table 1 shows the amino acid sequence identities between the Sav1866-TMD, ABCB1-TMD1, ABCB1-TMD2, ABCC5-TMD1, and ABCC5-TMD2. The TMD with the highest sequence identity with the Sav1866 TMD is the ABCC5-TMD1 (21 %), while the TMD with the lowest sequence identity with the Sav1866 TMD is the ABCC5-TMD2 (16 %). Both ABCB1 TMDs share a 17 % sequence identity with Sav1866. The percentages (%) of the charged amino acids aspartate (D), glutamate (E), histidine (H), lysine (K), and arginine (R), and positive-negative ratios of amino acids in the ABCB1 and ABCC5 TMDs are shown in Table 2. While the positive-negative ratio of amino acids is 1.1 (1.4 when histidine is included) in the ABCB1 TMDs, the corresponding ratio in the ABCC5 TMDs is 1.5 (1.8 when histidine is included). Thus the positive-negative ratio of amino acids is higher in the ABCC5 TMDs than in the ABCB1 TMDs. The charged amino acids were mainly localized in the substrate translocation chamber.
The refined ABCB1 and ABCC5 (submitted) models are shown in Figure 2, panels A and B. The loop connecting NBD1 and TMD2 of ABCB1 was mainly α-helical from residues 623–703, except from a parallel β-sheet formed between residues 614–618 and residues 646–650, and an extended stretch from residues 651–657. The first part of this loop was folded and covering NBD1 of ABCB1 towards the cytoplasm. A central cavity perpendicular to the cell membrane was formed by TMD1 and TMD2, and TMHs 1, 2, 3, 5, 6, 7, 8, 9, 11 and 12 contributed to the cavity lining. TMH5 and TMH2 of TMD1 were packed against TMH8 and TMH11 of TMD2, respectively, with mainly hydrophobic interactions. The substrate translocation chamber was closed towards the intracellular side, and the TMDs were twisted relative to the NBDs. The TMHs diverged into two symmetrical parts towards the extracellular side, one part consisting of TMHs 1 and 2 of TMD1 and TMHs 9–12 of TMD2, and one part consisting of TMHs 7 and 8 of TMD2 and TMHs 3–6 of TMD1 (Figure 2). Interactions between the NBDs were relatively hydrophilic, and the secondary structure of the areas of each NBD forming the contact area between the two NBDs was generally in extended conformation. The NBDs, having the same fold as the NBDs of the Sav1866 crystal structure, were tightly packed at the intracellular side of the membrane, containing the nucleotide binding sites formed by the motifs Walker A, Walker B, Q-loop and switch regions.
EPS of the substrate translocation chamber
Figure 3 shows the EPS of the substrate translocation chambers of ABCB1 (Panel A) and ABCC5 (Panel B). While the EPS of the substrate translocation chamber of ABCB1 was neutral with negative and weakly positive areas, the EPS of the ABCC5 substrate translocation chamber was generally positive.
The overall quality factor of ABCB1, as shown by the Errat option of the Savs Metaserver, was 96.2, and a value above 90 indicates a good model. According to the Ramachandran plot provided by the Procheck option, 87.1 % of the ABCB1 residues were in the most favored regions, 11.7 % were in additional allowed regions, 0.8 % were in generously allowed regions, and 0.4 % were in disallowed regions. The summary of the Whatcheck option reported that the model was satisfactory.
Several ABCB1 models have previously been published [33–36] based on MsbA X-ray crystal structures that later were retracted . The 12 TMHs of the present ABCB1 model are arranged as the TMHs of the Sav1866 crystal structure , and both are consistent with the electron density maps of ABCB1 .
As shown in Additional file 1, the human ABC efflux transporters comprise a large group of transporters featuring a wide range of functions and selectivities. ABC efflux transporters play important roles in physiological processes by transporting ligands such as bile salts/acids, conjugated steroids, cyclic nucleotides, ions, heme, lipids, antigens, retinoids, peptides, leukotrienes, organic anions, cations, and cholesterol, and many of them are involved in drug efflux.
Since ABCB1 transports cationic amphiphilic and lipophilic substrates [13–16] and ABCC5 transports organic anions [17, 18], the substrate translocation chamber localized in the TMDs of these transporters were of particular interest from a pharmacological point of view. The EPS of the substrate translocation chamber of ABCB1 was neutral with negative and weakly positive spots (Figure 3A). In contrast, the substrate translocation chamber of ABCC5 was generally positive (Figure 3B). An amino acid charge difference could also be seen between the TMDs of the two transporters (Table 2), with a lower positive-negative amino acid ratio of ABCB1 than of ABCC5. Thus, ABCB1, which transports cationic amphiphilic and lipophilic substrates, has a more neutral substrate translocation chamber than ABCC5, which has a positive chamber that transports organic anions. Substrates for these transporters bind to a binding site accessible to the intracellular side of the transporters. During the translocation process the binding site changes conformation, and the substrates are released to the extracellular side. The Sav1866 structure is captured in an outward facing conformation with the pore representing an extrusion pocket, rather than a binding pocket, and the modeled ABCB1 pore also represents an extrusion pocket. Even though the conformation changes, from a high affinity binding site (substrate recognition) to a low affinity binding site (substrate extrusion), the amino acids in the translocation area will be expected to contribute to similar ESP in both conformations.
Several cross-linking and site directed mutagenesis data have been published on ABCB1 [26–32]. These studies have indicated that TMH6 and TMH12 may take part in ligand binding [26, 27, 30, 31]. Cross-linking has also shown that TMH5 and TMH8 are near each other , and that TMH2 and TMH11 are near each other . As shown in Figure 3A, the present ABCB1 model is consistent with these experimental data; TMH6/TMH12, TMH5/TMH8 and TMH2/TMH11 are indeed adjacent. Comparing the reported residues from the experimental studies with the orientations of these residues in the present ABCB1 model verifies that the pore-lining residues of the TMHs are correctly localized, confirming that the alignment used for the ICM modeling procedure is realistic. Cross-linking studies have shown that residue pairs Asn266-Gly774, Ile299-Phe770, Ile299-Gly774, and Gly300-Phe770 (TMH5 and TMH8, respectively), are adjacent . In the present ABCB1 model, these residue pairs are in direct contact with each other. According to cross-linking studies, Val133 and Cys137 (TMH2) are close to Ala935 and Gly939 (TMH11) , and this is also in accordance with the ABCB1 model. Furthermore, experimental studies have suggested that Leu65 (TMH1) , Ile306 (TMH5) , Ile340 (TMH6) [26, 31], Phe343 (TMH6) , Phe728 (TMH7) , and Val982 (TMH12)  may participate in ligand binding. All these residues line the aqueous pore of the ABCB1 model and may indeed have ligand contact.
Site directed mutagenesis studies on ABCB1 have proposed a verapamil binding site including residues Leu65 (TMH1) , Ile306 (TMH5) , Ile340 (TMH6) [26, 31] and Phe343 (TMH6) . In the ABCB1 model these residues may form a binding site (Figure 4A). Ligand interactions between the TMH6 residues Ile340 and Phe343 and rhodamine have also been proposed in an ABCB1 modeling and docking study . The corresponding residues in ABCC5 are Gln190 (TMH1), Val410 (TMH5), Asn441 (TMH6) and Thr444 (TMH6), respectively (Figure 4B). Gln190 (TMH1), Asn441 (TMH6) and Thr444 (TMH6) of ABCC5 have previously been proposed to take part in ligand binding in a previous MRP5 modeling and cGMP docking study (submitted). Interestingly, the above mentioned ABCB1 residues are more lipophilic than the corresponding ABCC5 residues. This is in accordance with the lipophilic efflux featured by ABCB1, and with the more neutral EPS of the ABCB1 substrate translocation chambers.
Even though the EPS differences between the substrate translocation chambers of ABCB1 and ABCC5 are in accordance with their substrate specificity differences, one can not be certain that the Sav1866 crystal structure is a suitable template. According to the evolutionary tree (Figure 1), there are five main clusters of ABC efflux transporters: ABCA, ABCB, ABCC, ABCD and ABCG. Two main branches are seen, with ABCB, ABCC and ABCD in one branch, and ABCA and ABCG in the other branch. ABCB and ABCC subfamilies are closer related to each other than to the ABCD subfamily. Sav1866 is situated on the "ABCB-branch". The evolutionary tree thus indicates that an Sav1866 crystal structure may be a suitable template for at least ABCB1. The phylogeny of ABC transporters is based on homology of their NBDs , which is why the number of TMHs may differ within one subfamily, such as in subfamily ABCC5, where ABCC1 has 17 TMHs and ABCC5 has 12 TMHs. Even though ABCB1 and ABCC5 are localized on different branches in the evolutionary tree, they both have a common core domain organisation (TMD1-NBD1-TMD2-NDB2) . The identity between Sav1866 and ABCB1 is 31%. Accurate predictions can be made with an amino acid sequence similarity greater than 50 % between the target and the template protein, but even with very low homologies there may be considerable structural similarities, such as for the G-protein coupled receptors and bacteriorhodopsin, where the sequence similarities within the transmembrane regions are 6–11% . The conservation of the secondary structure elements is also relevant, since active sites and functional domains can have very similar geometries, even for distantly related proteins. The sequence identity between Sav1866 and ABCC5 is 23%, and phylogenetic analyses of ABC transporters have indicated that eukaryotic ABCB transporters (including ABCB1), ABCC transporters (including ABCC5), and bacterial ABC transporters have a common ancestor, and that they have similar domain organizations . Among the ABCC transporters, ABCC5 is most similar to ABCB1 , indicating that the Sav1866 X-ray crystal structure could also be used as a template for constructing an ABCC5 model by homology. The identity between the Sav1866-TMD and the ABCC5-TMD1 is actually higher (21%) than the identity between the Sav1866-TMD and the ABCB1-TMD1 (17%). In comparison, the sequence identity between the human serotonin transporter (SERT) and the crystal structure of the bacterial homologue from Aquifex aeolicus (LeuTAa) is ~20%, and several SERT models have been made been made using the LeuTAa as a template [41, 42]. The TMD sequence identities between the Sav1866-TMD and the ABCB1- and ABCC5-TMDs thus indicate that they have an overall similar organization and that the Sav1866-TMD may have been a suitable template for modeling the TMD segments of ABCB1 and ABCC5.
Membrane proteins may be highly flexible, metastable molecules, making them generally difficult to crystallize, and it has been suggested for the major facilitator transporter Escherichia coli lactose permease symporter (Lac Permease) that substrate binding in transporters may result in widespread conformational changes, and scissors like movements and sliding or tilting motions may occur during turnover . The crystal structure of Sav1866 indicates that domain swapping and subunit twisting takes place in the transport cycle . Thus, the substrate may be "pumped" from the inside of the membrane, binding with high affinity to the binding site, to the outside of the membrane, binding with low affinity, and thus being expelled to the extracellular space . It is therefore possible that the Sav1866 crystal structure represents a substrate expelling state where the binding site has changed drastically into a low affinity conformation through twisting and squeezing movements.
The calculations did not include water molecules or membrane phospholipids, and this omission may have influenced the model structure. The N- and C-terminals and two loops of ABCB1, the loop connecting TMH1 and TMH2, and the loop connecting NBD1 and TMD2, are relatively long and are not accounted for in the Sav1866 crystal structure. These segments are outside the limits for reliable loop generation via PDB searches and could not be predicted or modeled with confidence. Thus, The N- and C-terminals were not included in the model, but the two loops were included in order to get a more correct distribution of masses and electrostatics in the calculations than in a model with gaps were these loops are. Anyhow, it should be kept in mind that loops of such lengths modeled with computational techniques for loop modeling are relatively inaccurate, and, consequently, these were the most uncertain parts of the model. The fragment-based ab initio ROSETTA approach to the prediction of protein structure  may have been used, but the conformations of the modeled loops would still be too uncertain because of their lengths. Thus, the most certain regions of the ABCB1 model are the NBDs, because of their high level of sequence identity to the NBDs of Sav1866, and the TMD parts, which are in accordance with cross-linking and site directed mutagenesis data published on ABCB1 [26–32], confirming that porelining residues of the TMHs are correctly localized. The most uncertain parts are the loop connecting TMH1 and TMH2, and the loop connecting NBD1 and TMD2, which implies that these regions should only be considered as relatively crude approximations. Since the loop connecting NBD1 and TMD2 started 17 amino acids further towards the N-terminal, the NBD1 region had amino acids in its C-terminal end that was modeled as a loop instead of with homology to the NDB1 of Sav1866. Thus, the conformation of this 17 amino acid segment is uncertain, but this short segment does not include the Walker A and B motifs and is not a major part of NBD1. The loops are probably highly flexible, so any conformation generated by molecular modeling will only be a model of a temporary loop conformation. Anyhow, since the substrate binding area is of particular interest from a pharmacological point of view, focus was kept on the TMH area, and not the loops, in this molecular modeling study.
Making crystals of membrane proteins is in general technically difficult, and when no X-ray crystal structure is available, molecular modeling is a step forward towards structural knowledge of drug targets such as ABCB1 and ABCC5. In this study, the molecular concepts of the substrate specificity differences between ABCB1 and ABCC5 have been visualized using molecular modeling techniques. Even though there are uncertainties concerning the overall models, it seems that both site directed mutagenesis data [26, 27, 31] and the EPS in the substrate translocation chambers are in accordance with the transport of cationic amphiphilic and lipophilic substrates by ABCB1 [13–16], and the transport of organic anions by ABCC5 [17, 18]. This, and the consistency with the latest electron density map of ABCB1 , indicates that the Sav1866 X-ray structure  may serve as a suitable template for the ABCB1 and ABCC5 models. The ABCB1 model presented here is considered as a working tool to aid experimental studies. Eventually, membrane transport modulating agents may be developed, which may be used in the search for overcoming multidrug resistance.
Co-ordinates of the ABCB1 and ABCC5 models are available from the authors upon request.
Phylogenetic analysis of human ABC transporters
The Swiss-Prot Protein knowledgebase  and the TCDB  were used to retrieve fasta files of human ABC transporters, together with their Swiss-Prot accession codes, their synonyms and TCDB classification numbers. The ICM software version 3.4–4  was used to create a multiple sequence alignment and an evolutionary tree of the human ABC transporters, together with Sav1866. The ICM software creates evolutionary trees by the neighbor-joining method .
The amino acid sequence identities between the TMDs of Sav1866, ABCB1 and ABCC5 were retrieved using the ICM software. The start and endpoints of the TMDs were 14–298 (Sav1866), 177–451 (ABCC5, TMD1), 852–1147 (ABCC5, TMD2), 52–350 (ABCB1, TMD1), and 709–992 (ABCB1, TMD2). Positive-negative ratios of amino acids in the TMDs of ABCB1 and ABCC5 were also analyzed using the ICM software.
Homology modeling of ABCB1
The crystal structure of Sav1866  (pdb code 2HYD), which has a 3 Å resolution, was used as template to construct a homology model of ABCB1 (Swiss-Prot accession code P08183), using the ICM software versions 3.4–9b . T-COFFEE, Version 4.71 available at the Le Centre national de la recherche scientifiquewebsite , and ICM version 3.4–4 , were used to create multiple sequence alignments of human ABCB1, human ABCC5, human ABCC11 (SWISS-PROT accession number Q9BX80), human ABCC4 (SWISS-PROT accession number O15439), Sav1866 (SWISS-PROT accession number Q99T13), Vibrio cholerae MsbA (SWISS-PROT accession number Q9KQW9) and Escherichia coli MsbA (SWISS-PROT accession number P60752). The alignments were used as a basis, and adjusted in ICM for gaps for the input alignment in the ICM homology modeling module. To strengthen the sequence alignment, secondary structure predictions were performed to define the boundaries of the TMHs using the PredictProtein server for sequence analysis and structure prediction , and SWISS-PROT . The alignment of Sav1866, ABCB1 and ABCC5 is shown in Figure 5. The ICM homology modeling module constructs the model from a few core sections defined by the average of Cα atom positions in the conserved regions. Loops are constructed by searching within thousands of high quality structures in the PDB databank  by matching them in regard to sequence similarity and sterical interactions with the surroundings of the model. The best fitting loops are selected based on their relative energies. N- and C-terminals were not included in the models. Because of the length of the loop connecting NBD1 and TMD2, the loop was particularly difficult to model. In the generated models of the loop, the residues had a tendency to overlap with surrounding amino acids (sterical clashes), and more than 20 models was constructed before a model without sterical clashes was generated. In order to accomplish this, the start of the loop was moved one amino acid further towards the N-terminal direction, or the end of the loop was moved one amino acid further towards the C-terminal direction, in the ICM input alignment per modeling round, making the input loop longer until a model with no sterical clashes was generated. The alignment shown in figure 5 is the exact input alignment used for the final model. The construction of the ABCC5 model is described in a previous cGMP docking study (submitted).
Globally optimizing of the side-chain positions and annealing of the backbones were performed with the RefineModel macro of ICM. This macro first performs a side-chain conformational sampling using "Montecarlo fast" , a program module that samples conformational space of a molecule with the ICM global optimization procedure. Iterations of the procedure consist of a random move followed by a local energy minimization, followed by a complete energy calculation. Based on the energy and the temperature, the iteration is accepted or rejected. After the "Montecarlo fast" module, an iterative annealing of the backbone with tethers provided is performed. These tethers are harmonic restraints pulling an atom in the model to a static point in space represented by a corresponding atom in the template. Finally a second Monte Carlo side-chain sampling is performed. ECEPP3 charges  were used for the amino acids, and a surface based implicit solvation model  was included in the calculations.
The ABCB1 model was subjected to two subsequent energy minimizations by the AMBER 8.0 program package, using the leaprc.ff03 force field . The first energy minimization was performed with restrained backbone by 500 cycles of steepest descent minimization followed by 500 steps of conjugate gradient minimization, and the second energy minimization was performed with no restraints by 1000 cycles of steepest descent minimization followed by 1500 steps of conjugate gradient minimization. A 10 Å cut-off radius for nonbonded interactions and a dielectric multiplicative constant of 1.0 for the electrostatic interactions were used in these minimizations. Membrane molecules were not included in the model refinements. The electrostatic potential surface (EPS) of the ABCB1 model was calculated with the ICM program, with a potential scale from -10 to +10.
Quality validation of the ABCB1 model
The stereochemical quality of the ABCB1 model was checked using the Savs Metaserver for analyzing and validating protein structures . Programs run were Procheck , What_check , and Errat .
Saier MH: A functional-phylogenetic classification system for transmembrane solute transporters. Microbiol Mol Biol Rev. 2000, 64 (2): 354-411. 10.1128/MMBR.64.2.354-411.2000.
The transport classification database (TCDB). http://www.tcdb.org/index.php
Higgins CF, Linton KJ: Structural biology. The xyz of ABC transporters. Science. 2001, 293 (5536): 1782-1784. 10.1126/science.1065588.
Oswald C, Holland IB, L. S: The motor domains of ABC-transporters - What can structures tell us?. Naunyn-Schmiedeberg’s Arch Pharmacol. 2006, 372: 385–399-
Dean M, Rzhetsky A, Allikmets R: The human ATP-binding cassette (ABC) transporter superfamily. Genome Res. 2001, 11 (7): 1156-1166. 10.1101/gr.GR-1649R.
Dassa E, Bouige P: The ABC of ABCS: a phylogenetic and functional classification of ABC systems in living organisms. Res Microbiol. 2001, 152 (3-4): 211-229. 10.1016/S0923-2508(01)01194-9.
Walker JE, Saraste M, Runswick MJ, Gay NJ: Distantly related sequences in the alpha- and beta-subunits of ATP synthase, myosin, kinases and other ATP-requiring enzymes and a common nucleotide binding fold. Embo J. 1982, 1 (8): 945-951.
Ross DD, Doyle LA: Mining our ABCs: pharmacogenomic approach for evaluating transporter function in cancer drug resistance. Cancer Cell. 2004, 6 (2): 105-107. 10.1016/j.ccr.2004.08.003.
Leslie EM, Deeley RG, Cole SP: Multidrug resistance proteins: role of P-glycoprotein, MRP1, MRP2, and BCRP (ABCG2) in tissue defense. Toxicol Appl Pharmacol. 2005, 204 (3): 216-237. 10.1016/j.taap.2004.10.012.
Sampath J, Adachi M, Hatse S, Naesens L, Balzarini J, Flatley RM, Matherly LH, Schuetz JD: Role of MRP4 and MRP5 in biology and chemotherapy. AAPS PharmSci. 2002, 4 (3): E14-10.1208/ps040314.
Ejendal KF, Hrycyna CA: Multidrug resistance and cancer: the role of the human ABC transporter ABCG2. Curr Protein Pept Sci. 2002, 3 (5): 503-511. 10.2174/1389203023380521.
Dantzig AH, de Alwis DP, Burgess M: Considerations in the design and development of transport inhibitors as adjuncts to drug therapy. Adv Drug Deliv Rev. 2003, 55 (1): 133-150. 10.1016/S0169-409X(02)00175-8.
Smit JW, Duin E, Steen H, Oosting R, Roggeveld J, Meijer DK: Interactions between P-glycoprotein substrates and other cationic drugs at the hepatic excretory level. Br J Pharmacol. 1998, 123 (3): 361-370. 10.1038/sj.bjp.0701606.
Orlowski S, Garrigos M: Multiple recognition of various amphiphilic molecules by the multidrug resistance P-glycoprotein: molecular mechanisms and pharmacological consequences coming from functional interactions between various drugs. Anticancer Res. 1999, 19 (4B): 3109-3123.
Muller M, Mayer R, Hero U, Keppler D: ATP-dependent transport of amphiphilic cations across the hepatocyte canalicular membrane mediated by mdr1 P-glycoprotein. FEBS Lett. 1994, 343 (2): 168-172. 10.1016/0014-5793(94)80312-9.
Wang EJ, Lew K, Casciano CN, Clement RP, Johnson WW: Interaction of common azole antifungals with P glycoprotein. Antimicrob Agents Chemother. 2002, 46 (1): 160-165. 10.1128/AAC.46.1.160-165.2002.
Jedlitschky G, Burchell B, Keppler D: The multidrug resistance protein 5 functions as an ATP-dependent export pump for cyclic nucleotides. J Biol Chem. 2000, 275 (39): 30069-30074. 10.1074/jbc.M005463200.
Borst P, de Wolf C, van de Wetering K: Multidrug resistance-associated proteins 3, 4, and 5. Pflugers Arch. 2007, 453 (5): 661-673. 10.1007/s00424-006-0054-9.
Rosenberg MF, Callaghan R, Ford RC, Higgins CF: Structure of the multidrug resistance P-glycoprotein to 2.5 nm resolution determined by electron microscopy and image analysis. J Biol Chem. 1997, 272 (16): 10685-10694. 10.1074/jbc.272.16.10685.
Rosenberg MF, Mao Q, Holzenburg A, Ford RC, Deeley RG, Cole SP: The structure of the multidrug resistance protein 1 (MRP1/ABCC1). crystallization and single-particle analysis. J Biol Chem. 2001, 276 (19): 16076-16082. 10.1074/jbc.M100176200.
Rosenberg MF, Kamis AB, Callaghan R, Higgins CF, Ford RC: Three-dimensional structures of the mammalian multidrug resistance P-glycoprotein demonstrate major conformational changes in the transmembrane domains upon nucleotide binding. J Biol Chem. 2003, 278 (10): 8294-8299. 10.1074/jbc.M211758200.
Rosenberg MF, Callaghan R, Modok S, Higgins CF, Ford RC: Three-dimensional structure of P-glycoprotein: the transmembrane regions adopt an asymmetric configuration in the nucleotide-bound state. J Biol Chem. 2005, 280 (4): 2857-2862. 10.1074/jbc.M410296200.
Dawson RJ, Locher KP: Structure of a bacterial multidrug ABC transporter. Nature. 2006
Senior AE, Gadsby DC: ATP hydrolysis cycles and mechanism in P-glycoprotein and CFTR. Semin Cancer Biol. 1997, 8 (3): 143-150. 10.1006/scbi.1997.0065.
Borst P, Evers R, Kool M, Wijnholds J: A family of drug transporters: the multidrug resistance-associated proteins. J Natl Cancer Inst. 2000, 92 (16): 1295-1302. 10.1093/jnci/92.16.1295.
Loo TW, Clarke DM: Location of the rhodamine-binding site in the human multidrug resistance P-glycoprotein. J Biol Chem. 2002, 277 (46): 44332-44338. 10.1074/jbc.M208433200.
Loo TW, Bartlett MC, Clarke DM: Methanethiosulfonate derivatives of rhodamine and verapamil activate human P-glycoprotein at different sites. J Biol Chem. 2003, 278 (50): 50136-50141. 10.1074/jbc.M310448200.
Loo TW, Bartlett MC, Clarke DM: Disulfide cross-linking analysis shows that transmembrane segments 5 and 8 of human P-glycoprotein are close together on the cytoplasmic side of the membrane. J Biol Chem. 2004, 279 (9): 7692-7697. 10.1074/jbc.M311825200.
Loo TW, Bartlett MC, Clarke DM: Val133 and Cys137 in transmembrane segment 2 are close to Arg935 and Gly939 in transmembrane segment 11 of human P-glycoprotein. J Biol Chem. 2004, 279 (18): 18232-18238. 10.1074/jbc.M400229200.
Loo TW, Clarke DM: Recent progress in understanding the mechanism of P-glycoprotein-mediated drug efflux. J Membr Biol. 2005, 206 (3): 173-185. 10.1007/s00232-005-0792-1.
Loo TW, Bartlett MC, Clarke DM: Transmembrane segment 1 of human P-glycoprotein contributes to the drug-binding pocket. Biochem J. 2006, 396 (3): 537-545. 10.1042/BJ20060012.
Loo TW, Bartlett MC, Clarke DM: Transmembrane segment 7 of human P-glycoprotein forms part of the drug-binding pocket. Biochem J. 2006
Vandevuer S, Van Bambeke F, Tulkens PM, Prevost M: Predicting the three-dimensional structure of human P-glycoprotein in absence of ATP by computational techniques embodying crosslinking data: Insight into the mechanism of ligand migration and binding sites. Proteins. 2006, 63 (3): 466-478. 10.1002/prot.20892.
Seigneuret M, Garnier-Suillerot A: A structural model for the open conformation of the mdr1 P-glycoprotein based on the MsbA crystal structure. J Biol Chem. 2003, 278 (32): 30115-30124. 10.1074/jbc.M302443200.
Stenham DR, Campbell JD, Sansom MS, Higgins CF, Kerr ID, Linton KJ: An atomic detail model for the human ATP binding cassette transporter P-glycoprotein derived from disulfide cross-linking and homology modeling. Faseb J. 2003, 17 (15): 2287-2289.
Shilling RA, Balakrishnan L, Shahi S, Venter H, van Veen HW: A new dimer interface for an ABC transporter. Int J Antimicrob Agents. 2003, 22 (3): 200-204. 10.1016/S0924-8579(03)00212-7.
Chang G, Roth CB, Reyes CL, Pornillos O, Chen YJ, Chen AP: Retraction. Science. 2006, 314 (5807): 1875-10.1126/science.314.5807.1875b.
Hibert MF, Trumpp-Kallmeyer S, Bruinvels A, Hoflack J: Three-dimensional models of neurotransmitter G-binding protein-coupled receptors. Mol Pharmacol. 1991, 40 (1): 8-15.
Igarashi Y, Aoki KF, Mamitsuka H, Kuma K, Kanehisa M: The evolutionary repertoires of the eukaryotic-type ABC transporters in terms of the phylogeny of ATP-binding domains in eukaryotes and prokaryotes. Mol Biol Evol. 2004, 21 (11): 2149-2160. 10.1093/molbev/msh226.
Borst P, Evers R, Kool M, Wijnholds J: The multidrug resistance protein family. Biochim Biophys Acta. 1999, 1461 (2): 347-357. 10.1016/S0005-2736(99)00167-4.
Beuming T, Shi L, Javitch JA, Weinstein H: A comprehensive structure-based alignment of prokaryotic and eukaryotic neurotransmitter/Na+ symporters (NSS) aids in the use of the LeuT structure to probe NSS structure and function. Mol Pharmacol. 2006
Ravna AW, Jaronczyk M, Sylte I: A homology model of SERT based on the LeuT(Aa) template. Bioorg Med Chem Lett. 2006, 16 (21): 5594-5597. 10.1016/j.bmcl.2006.08.028.
Kaback HR, Wu J: From membrane to molecule to the third amino acid from the left with a membrane transport protein. Q Rev Biophys. 1997, 30 (4): 333-364. 10.1017/S0033583597003387.
Hammes GG: Unifying concept for the coupling between ion pumping and ATP hydrolysis or synthesis. Proc Natl Acad Sci U S A. 1982, 79 (22): 6881-6884. 10.1073/pnas.79.22.6881.
Simons KT, Bonneau R, Ruczinski I, Baker D: Ab initio protein structure prediction of CASP III targets using ROSETTA. Proteins. 1999, Suppl 3: 171-176. 10.1002/(SICI)1097-0134(1999)37:3+<171::AID-PROT21>3.0.CO;2-Z.
Bairoch A, Apweiler R: The SWISS-PROT protein sequence data bank and its supplement TrEMBL in 1999. Nucleic Acids Res. 1999, 27 (1): 49-54. 10.1093/nar/27.1.49.
Abagyan R, Totrov M, Kuznetsov DN: ICM - a new method for protein modeling and design. Applications to docking and structure prediction from the distorted native comformation. J Comp Chem. 1994, 15: 488-506. 10.1002/jcc.540150503.
Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.
T-COFFE version 4.71.http://tcoffee.vital-it.ch/cgi-bin/Tcoffee/tcoffee_cgi/index.cgi
Rost B: PHD: Predicting one-dimensional protein structure by profile-based neural networks. Meth Enzymol. 1996, 266: 525-539.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28 (1): 235-242. 10.1093/nar/28.1.235.
Abagyan R, Totrov M: Biased probability Monte Carlo conformational searches and electrostatic calculations for peptides and proteins. J Mol Biol. 1994, 235 (3): 983-1002. 10.1006/jmbi.1994.1052.
Nemethy G, Gibson KD, Palmer KA, Yoon CN, Paterlini G, Zagari A, Rumsey S, Scheraga HA: Energy parameters in polypeptides. 10. Improved geometrical parameters and nonbonded interactions to use in the ECEPP/3 algorithms, with applications to proline-containing peptides. J Phys Chem. 1992, 96: 6472-6484. 10.1021/j100194a068.
Case DA, Darden TA, Cheatham III TE, Simmerling CL, Wang J, Duke RE, Luo R, Merz KM, Wang B, Pearlman DA, Crowley M, Brozell S, Tsui V, Gohlke H, Mongan J, Hornak V, Cui G, Beroza P, Schafmeister C, Caldwell JW, Ross WS, Kollman PA: AMBER 8. 2004, San Francisco , University of California
Savs Metaserver for structure validation.http://nihserver.mbi.ucla.edu/SAVS/
Laskoswki RA, MacArthur MW, Moss DS, Thorton JM: PROCHECK: a program to check the stereochemical quality of protein structures. J Appl Cryst. 1993, 26: 283-291. 10.1107/S0021889892009944.
Hooft RW, Vriend G, Sander C, Abola EE: Errors in protein structures. Nature. 1996, 381 (6580): 272-10.1038/381272a0.
Colovos C, Yeates TO: Verification of protein structures: patterns of nonbonded atomic interactions. Protein Sci. 1993, 2 (9): 1511-1519.
This study was sponsored by the Norwegian Cancer Society.
The author(s) declare that they have no competing interests.
AWR carried out the molecular modeling studies (homology modeling and model refinement, and quality validation), created sequence alignments and the evolutionary tree, and drafted the manuscript. IS participated in the design of the study, and contributed with bioinformatics advice and critical review of the manuscript. GS conceived of the study and participated in its design, contributed with biological advice and critical review of the manuscript. All authors read and approved the final manuscript.