Data-driven homology modelling of P-glycoprotein in the ATP-bound state indicates flexibility of the transmembrane domains

Comprehensive Collection and Prediction of ABC Transmembrane Protein Structures in the AI Era of Structural Biology

The number of unique transmembrane (TM) protein structures doubled in the last four years, which can be attributed to the revolution of cryo-electron microscopy. In addition, AlphaFold2 (AF2) also provided a large number of predicted structures with high quality. However, if a specific protein family is the subject of a study, collecting the structures of the family members is highly challenging in spite of existing general and protein domain-specific databases. Here, we demonstrate this and assess the applicability and usability of automatic collection and presentation of protein structures via the ABC protein superfamily. Our pipeline identifies and classifies transmembrane ABC protein structures using the PFAM search and also aims to determine their conformational states based on special geometric measures, conftors. Since the AlphaFold database contains structure predictions only for single polypeptide chains, we performed AF2-Multimer predictions for human ABC half transporters functioning as dimers. Our AF2 predictions warn of possibly ambiguous interpretation of some biochemical data regarding interaction partners and call for further experiments and experimental structure determination. We made our predicted ABC protein structures available through a web application, and we joined the 3D-Beacons Network to reach the broader scientific community through platforms such as PDBe-KB.

Efflux Pump Mediated Antimicrobial Resistance by Staphylococci in Health-Related Environments: Challenges and the Quest for Inhibition

The increasing emergence of antimicrobial resistance in staphylococcal bacteria is a major health threat worldwide due to significant morbidity and mortality resulting from their associated hospital- or community-acquired infections. Dramatic decrease in the discovery of new antibiotics from the pharmaceutical industry coupled with increased use of sanitisers and disinfectants due to the ongoing COVID-19 pandemic can further aggravate the problem of antimicrobial resistance. Staphylococci utilise multiple mechanisms to circumvent the effects of antimicrobials. One of these resistance mechanisms is the export of antimicrobial agents through the activity of membrane-embedded multidrug efflux pump proteins. The use of efflux pump inhibitors in combination with currently approved antimicrobials is a promising strategy to potentiate their clinical efficacy against resistant strains of staphylococci, and simultaneously reduce the selection of resistant mutants. This review presents an overview of the current knowledge of staphylococcal efflux pumps, discusses their clinical impact, and summarises compounds found in the last decade from plant and synthetic origin that have the potential to be used as adjuvants to antibiotic therapy against multidrug resistant staphylococci. Critically, future high-resolution structures of staphylococcal efflux pumps could aid in design and development of safer, more target-specific and highly potent efflux pump inhibitors to progress into clinical use.

Human ABCB1 with an ABCB11-like degenerate nucleotide binding site maintains transport activity by avoiding nucleotide occlusion

Several ABC exporters carry a degenerate nucleotide binding site (NBS) that is unable to hydrolyze ATP at a rate sufficient for sustaining transport activity. A hallmark of a degenerate NBS is the lack of the catalytic glutamate in the Walker B motif in the nucleotide binding domain (NBD). The multidrug resistance transporter ABCB1 (P-glycoprotein) has two canonical NBSs, and mutation of the catalytic glutamate E556 in NBS1 renders ABCB1 transport-incompetent. In contrast, the closely related bile salt export pump ABCB11 (BSEP), which shares 49% sequence identity with ABCB1, naturally contains a methionine in place of the catalytic glutamate. The NBD-NBD interfaces of ABCB1 and ABCB11 differ only in four residues, all within NBS1. Mutation of the catalytic glutamate in ABCB1 results in the occlusion of ATP in NBS1, leading to the arrest of the transport cycle. Here we show that despite the catalytic glutamate mutation (E556M), ABCB1 regains its ATP-dependent transport activity, when three additional diverging residues are also replaced. Molecular dynamics simulations revealed that the rescue of ATPase activity is due to the modified geometry of NBS1, resulting in a weaker interaction with ATP, which allows the quadruple mutant to evade the conformationally locked pre-hydrolytic state to proceed to ATP-driven transport. In summary, we show that ABCB1 can be transformed into an active transporter with only one functional catalytic site by preventing the formation of the ATP-locked pre-hydrolytic state in the non-canonical site.

ABC transporters are one of the largest membrane protein superfamilies, present in all organisms from archaea to humans. They transport a wide range of molecules including amino acids, sugars, vitamins, nucleotides, peptides, lipids, metabolites, antibiotics, and xenobiotics. ABC transporters energize substrate transport by hydrolyzing ATP in two symmetrically arranged nucleotide binding sites (NBSs). The human multidrug resistance transporter ABCB1 has two active NBSs, and it is generally believed that integrity and cooperation of both sites are needed for transport. Several human ABC transporters, such as the bile salt transporter ABCB11, have one degenerate NBS, which has significantly reduced ATPase activity. Interestingly, unilateral mutations affecting one of the two NBSs completely abolish the function of symmetrical ABC transporters. Here we engineered an ABCB1 variant with a degenerate, ABCB11-like NBS1, which can nevertheless transport substrates. Our results indicate that ABCB1 can mediate active transport with a single active site, questioning the validity of models assuming strictly alternating catalysis.

Cross-linking, DEER-spectroscopy and molecular dynamics confirm the inward facing state of P-glycoprotein in a lipid membrane

The drug efflux pump P-glycoprotein (P-gp) displays a complex transport mechanism involving multiple drug binding sites and two centres for nucleotide hydrolysis. Elucidating the molecular mechanism of transport remains elusive and the availability of P-gp structures in distinct natural and ligand trapped conformations will accelerate our understanding. The present investigation sought to provide biochemical data to validate specific features of these structures; with particular focus on the transmembrane domain that provides the transport conduit. Hence our focus was on transmembrane helices six and twelve (TM6/TM12), which are believed to participate in drug binding, as they line the central transport conduit and provide a direct link to the catalytic centres. A series of P-gp mutants were generated with a single cysteine in both TM6 and TM12 to facilitate measurement of inter-helical distances using cross-linking and DEER strategies. Experimental results were compared to published structures per se and those refined by MD simulations. This analysis revealed that the refined inward-facing murine structure (4M1M) of P-gp provides a good representation of the proximity, topography and relative motions of TM6 and TM12 in reconstituted human P-gp.

Conversion of chemical to mechanical energy by the nucleotide binding domains of ABCB1

P-glycoprotein (ABCB1) is an important component of barrier tissues that extrudes a wide range of chemically unrelated compounds. ABCB1 consists of two transmembrane domains forming the substrate binding and translocation domain, and of two cytoplasmic nucleotide binding domains (NBDs) that provide the energy by binding and hydrolyzing ATP. We analyzed the mechanistic and energetic properties of the NBD dimer via molecular dynamics simulations. We find that MgATP stabilizes the NBD dimer through strong attractive forces by serving as an interaction hub. The irreversible ATP hydrolysis step converts the chemical energy stored in the phosphate bonds of ATP into potential energy. Following ATP hydrolysis, interactions between the NBDs and the ATP hydrolysis products MgADP + Pi remain strong, mainly because Mg2+ forms stabilizing interactions with ADP and Pi. Despite these stabilizing interactions MgADP + Pi are unable to hold the dimer together, which becomes separated by avid interactions of MgADP + Pi with water. ATP binding to the open NBDs and ATP hydrolysis in the closed NBD dimer represent two steps of energy input, each leading to the formation of a high energy state. Relaxation from these high energy states occurs through conformational changes that push ABCB1 through the transport cycle.

Quantitative comparison of ABC membrane protein type I exporter structures in a standardized way

An increasing number of ABC membrane protein structures are determined by cryo-electron microscopy and X-ray crystallography, consequently identifying differences between their conformations has become an arising issue. Therefore, we propose to define standardized measures for ABC Type I exporter structure characterization. We set conformational vectors, conftors, which describe the relative orientation of domains and can highlight structural differences. In addition, continuum electrostatics calculations were performed to characterize the energetics of membrane insertion illuminating functionally crucial regions. In summary, the proposed metrics contribute to deeper understanding of ABC membrane proteins' structural features, structure validation, and analysis of movements observed in a molecular dynamics trajectory. Moreover, the concept of standardized metrics can be applied not only to ABC membrane protein structures (http://conftors.hegelab.org).

Dissecting the Forces that Dominate Dimerization of the Nucleotide Binding Domains of ABCB1

P-glycoprotein, also known as multidrug resistance protein 1 or ABCB1, can export a wide range of chemically unrelated compounds, including chemotherapeutic drugs. ABCB1 consists of two transmembrane domains that form the substrate binding and translocation domain, and of two cytoplasmic nucleotide binding domains (NBDs) that energize substrate transport by ATP binding and hydrolysis. ATP binding triggers dimerization of the NBDs, which switches the transporter from an inward facing to an outward facing transmembrane domain conformation. We performed MD simulations to study the dynamic behavior of the NBD dimer in the presence or absence of nucleotides. In the apo configuration, the NBDs were overall attractive to each other as shown in the potential of mean force profile, but the energy well was shallow and broad. In contrast, a sharp and deep energy minimum (∼-42 kJ/mol) was found in the presence of ATP, leading to a well-defined conformation. Motif interaction network analyses revealed that ATP stabilizes the NBD dimer by serving as the central hub for interdomain connections. Simulations showed that forces promoting dimerization are multilayered, dominated by electrostatic interactions between the nucleotide and conserved amino acids of the signature sequence and the Walker A motif. In addition, direct and water-bridged hydrogen bonds between NBDs provided conformation-defining interactions. Importantly, we characterized a largely unrecognized but essential contribution from hydrophobic interactions between the adenine moiety of the nucleotides and a hydrophobic surface of the X-loop to the stabilization of the nucleotide-bound NBD dimer. These hydrophobic interactions lead to a sharp energy minimum, thereby conformationally restricting the nucleotide-bound state.

Molecular Mechanism of Taurocholate Transport by the Bile Salt Export Pump, an ABC Transporter Associated with Intrahepatic Cholestasis

The bile salt export pump (BSEP/ABCB11) transports bile salts from hepatocytes into bile canaliculi. Its malfunction is associated with severe liver disease. One reason for functional impairment of BSEP is systemic administration of drugs, which as a side effect inhibit the transporter. Therefore, drug candidates are routinely screened for potential interaction with this transporter. Hence, understanding the functional biology of BSEP is of key importance. In this study, we engineered the transporter to dissect interdomain communication paths. We introduced mutations in noncanonical and in conserved residues of either of the two nucleotide binding domains and determined the effect on BSEP basal and substrate-stimulated ATPase activity as well as on taurocholate transport. Replacement of the noncanonical methionine residue M584 (Walker B sequence of nucleotide binding site 1) by glutamate imparted hydrolysis competency to this site. Importantly, this mutation was able to sustain 15% of wild-type transport activity, when the catalytic glutamate of the canonical nucleotide binding site 2 was mutated to glutamine. Kinetic modeling of experimental results for the ensuing M584E/E1244Q mutant suggests that a transfer of hydrolytic capacity from the canonical to the noncanonical nucleotide binding site results in loss of active and adoption of facilitative characteristics. This facilitative transport is ATP-gated. To the best of our knowledge, this result is unprecedented in ATP-binding cassette proteins with one noncanonical nucleotide binding site. Our study promotes an understanding of the domain interplay in BSEP as a basis for exploration of drug interactions with this transporter.

Folding correction of ABC-transporter ABCB1 by pharmacological chaperones: a mechanistic concept

Point mutations of ATP‐binding cassette (ABC) proteins are a common cause of human diseases. Available crystal structures indicate a similarity in the architecture of several members of this protein family. Their molecular architecture makes these proteins vulnerable to mutation, when critical structural elements are affected. The latter preferentially involve the two transmembrane domain (TMD)/nucleotide‐binding domain (NBD) interfaces (transmission interfaces), formation of which requires engagement of coupling helices of intracellular loops with NBDs. Both, formation of the active sites and engagement of the coupling helices, are contingent on correct positioning of ICLs 2 and 4 and thus an important prerequisite for proper folding. Here, we show that active site compounds are capable of rescuing P‐glycoprotein (P‐gp) mutants ∆Y490 and ∆Y1133 in a concentration‐dependent manner. These trafficking deficient mutations are located at the transmission interface in pseudosymmetric position to each other. In addition, the ability of propafenone analogs to correct folding correlates with their ability to inhibit transport of model substrates. This finding indicates that folding correction and transport inhibition by propafenone analogs are brought about by binding to the active sites. Furthermore, this study demonstrates an asymmetry in folding correction with cis‐flupentixol, which reflects the asymmetric binding properties of this modulator to P‐gp. Our results suggest a mechanistic model for corrector action in a model ABC transporter based on insights into the molecular architecture of these transporters.

Unidirectional Transport Mechanism in an ATP Dependent Exporter

ATP-binding cassette (ABC) transporters use the energy of ATP binding and hydrolysis to move a large variety of compounds across biological membranes. P-glycoprotein, involved in multidrug resistance, is the most investigated eukaryotic family member. Although a large number of biochemical and structural approaches have provided important information, the conformational dynamics underlying the coupling between ATP binding/hydrolysis and allocrite transport remains elusive. To tackle this issue, we performed molecular dynamic simulations for different nucleotide occupancy states of Sav1866, a prokaryotic P-glycoprotein homologue. The simulations reveal an outward-closed conformation of the transmembrane domain that is stabilized by the binding of two ATP molecules. The hydrolysis of a single ATP leads the X-loop, a key motif of the ATP binding cassette, to interfere with the transmembrane domain and favor its outward-open conformation. Our findings provide a structural basis for the unidirectionality of transport in ABC exporters and suggest a ratio of one ATP hydrolyzed per transport cycle.

Interactions and cooperativity between P-glycoprotein structural domains determined by thermal unfolding provides insights into its solution structure and function

Structural changes in mouse P-glycoprotein (Pgp) induced by thermal unfolding were studied by differential scanning calorimetry (DSC), circular dichroism and fluorescence spectroscopy to gain insight into the solution conformation(s) of this ABC transporter that may not be apparent from current crystal structures. DSC of reconstituted Pgp showed two thermal unfolding transitions in the absence of MgATP, suggesting that each transition involved the cooperative unfolding of two or more interacting structural domains. A low calorimetric unfolding enthalpy and minimal structural changes were observed, which are hallmarks of the thermal unfolding of α-helical membrane proteins, because generally only the extramembranous regions undergo significant unfolding. Nucleotide binding increased the unfolding temperature of both transitions to the same extent, suggesting that one nucleotide binding domain (NBD) unfolds with each transition. Combined with the results from the two isolated NBDs, we propose that each DSC transition represents the cooperative unfolding of one NBD and the two contacting intracellular loops. Further, the presence of two transitions in both apo and MgATP bound wild-type Pgp suggests the NBD-dimeric conformation is transient, and that Pgp resides predominantly in the crystallographically observed inward-facing conformation with NBDs separated, even under conditions supporting continuous MgATP hydrolysis. In contrast, DSC of the vanadate-trapped MgADP·Pgp complex and the MgATP-bound catalytically inactive mutant, E552A/E1197A, show an additional transition at much higher temperature, corresponding to the unfolding of the nucleotide-trapped NBD-dimeric outward-facing conformation. The collective results indicate a strong preference for an NBD dissociated, inward-facing conformation of Pgp.

Molecular insight into conformational transmission of human P-glycoprotein

P-glycoprotein (P-gp), a kind of ATP-binding cassette transporter, can export candidates through a channel at the two transmembrane domains (TMDs) across the cell membranes using the energy released from ATP hydrolysis at the two nucleotide-binding domains (NBDs). Considerable evidence has indicated that human P-gp undergoes large-scale conformational changes to export a wide variety of anti-cancer drugs out of the cancer cells. However, molecular mechanism of the conformational transmission of human P-gp from the NBDs to the TMDs is still unclear. Herein, targeted molecular dynamics simulations were performed to explore the atomic detail of the conformational transmission of human P-gp. It is confirmed that the conformational transition from the inward- to outward-facing is initiated by the movement of the NBDs. It is found that the two NBDs move both on the two directions (x and y). The movement on the x direction leads to the closure of the NBDs, while the movement on the y direction adjusts the conformations of the NBDs to form the correct ATP binding pockets. Six key segments (KSs) protruding from the TMDs to interact with the NBDs are identified. The relative movement of the KSs along the y axis driven by the NBDs can be transmitted through α-helices to the rest of the TMDs, rendering the TMDs to open towards periplasm in the outward-facing conformation. Twenty eight key residue pairs are identified to participate in the interaction network that contributes to the conformational transmission from the NBDs to the TMDs of human P-gp. In addition, 9 key residues in each NBD are also identified. The studies have thus provided clear insight into the conformational transmission from the NBDs to the TMDs in human P-gp.

The Q loops of the human multidrug resistance transporter ABCB1 are necessary to couple drug binding to the ATP catalytic cycle

For a primary active pump, such as the human ATP-binding-cassette (ABC) transporter ABCB1, coupling of drug-binding by the two transmembrane domains (TMDs) to the ATP catalytic cycle of the two nucleotide-binding domains (NBDs) is fundamental to the transport mechanism, but is poorly understood at the biochemical level. Structure data suggest that signals are transduced through intracellular loops of the TMDs that slot into grooves on the NBDs. At the base of these grooves is the Q loop. We therefore mutated the eponymous glutamine in one or both NBD Q loops and measured the effect on conformation and function by using a conformation-sensitive antibody (UIC2) and a fluorescent drug (Bodipy-verapamil), respectively. We showed that the double mutant is trapped in the inward-open state, which binds the drug, but cannot couple to the ATPase cycle. Our data also describe marked redundancy within the transport mechanism, because single-Q-loop mutants are functional for Bodipy-verapamil transport. This result allowed us to elucidate transduction pathways from twin drug-binding cavities to the Q loops using point mutations to favor one cavity over the other. Together, the data show that the Q loop is the central flexion point where the aspect of the drug-binding cavities is coupled to the ATP catalytic cycle.

Pore-exposed tyrosine residues of P-glycoprotein are important hydrogen-bonding partners for drugs

The multispecific efflux transporter, P-glycoprotein, plays an important role in drug disposition. Substrate translocation occurs along the interface of its transmembrane domains. The rotational C2 symmetry of ATP-binding cassette transporters implies the existence of two symmetry-related sets of substrate-interacting amino acids. These sets are identical in homodimeric transporters, and remain evolutionary related in full transporters, such as P-glycoprotein, in which substrates bind preferentially, but nonexclusively, to one of two binding sites. We explored the role of pore-exposed tyrosines for hydrogen-bonding interactions with propafenone type ligands in their preferred binding site 2. Tyrosine 953 is shown to form hydrogen bonds not only with propafenone analogs, but also with the preferred site 1 substrate rhodamine123. Furthermore, an accessory role of tyrosine 950 for binding of selected propafenone analogs is demonstrated. The present study demonstrates the importance of domain interface tyrosine residues for interaction of small molecules with P-glycoprotein.

Overview of experimental models of the blood-brain barrier in CNS drug discovery

The blood-brain barrier (BBB) is a physical and metabolic entity that isolates the brain from the systemic circulation. The barrier consists of tight junctions between endothelial cells that contain egress transporters and catabolic enzymes. To cross the BBB, a drug must possess the appropriate physicochemical properties to achieve a sufficient time-concentration profile in brain interstitial fluid (ISF). In this overview, we review techniques to measure BBB permeation, which is evidenced by the free concentration of compound in brain ISF over time. We consider a number of measurement techniques, including in vivo microdialysis and brain receptor occupancy following perfusion. Consideration is also given to the endothelial and nonendothelial cell systems used to assess both the BBB permeation of a test compound and its interactions with egress transporters, and computer models employed for predicting passive permeation and the probability of interactions with BBB transporters.

Substrate versus inhibitor dynamics of P-glycoprotein

By far the most studied multidrug resistance protein is P-glycoprotein. Despite recent structural data, key questions about its function remain. P-glycoprotein (P-gp) is flexible and undergoes large conformational changes as part of its function and in this respect, details not only of the export cycle, but also the recognition stage are currently lacking. Given the flexibility, molecular dynamics (MD) simulations provide an ideal tool to examine this aspect in detail. We have performed MD simulations to examine the behaviour of P-gp. In agreement with previous reports, we found that P-gp undergoes large conformational changes which tended to result in the nucleotide-binding domains coming closer together. In all simulations, the approach of the NBDs was asymmetrical in agreement with previous observations for other ABC transporter proteins. To validate the simulations, we make extensive comparison to previous cross-linking data. Our results show very good agreement with the available data. We then went on to compare the influence of inhibitor compounds bound with simulations of a substrate (daunorubicin) bound. Our results suggest that inhibitors may work by keeping the NBDs apart, thus preventing ATP-hydrolysis. On the other hand, repeat simulations of daunorubicin (substrate) in one particular binding pose suggest that the approach of the NBDs is not impaired and that the structure would be still be competent to perform ATP hydrolysis, thus providing a model for inhibition or substrate transport. Finally we compare the latter to earlier QSAR data to provide a model for the first part of substrate transport within P-gp.

Prediction of promiscuous p-glycoprotein inhibition using a novel machine learning scheme

BACKGROUND

P-glycoprotein (P-gp) is an ATP-dependent membrane transporter that plays a pivotal role in eliminating xenobiotics by active extrusion of xenobiotics from the cell. Multidrug resistance (MDR) is highly associated with the over-expression of P-gp by cells, resulting in increased efflux of chemotherapeutical agents and reduction of intracellular drug accumulation. It is of clinical importance to develop a P-gp inhibition predictive model in the process of drug discovery and development.

METHODOLOGY/PRINCIPAL FINDINGS

An in silico model was derived to predict the inhibition of P-gp using the newly invented pharmacophore ensemble/support vector machine (PhE/SVM) scheme based on the data compiled from the literature. The predictions by the PhE/SVM model were found to be in good agreement with the observed values for those structurally diverse molecules in the training set (n = 31, r(2) = 0.89, q(2) = 0.86, RMSE = 0.40, s = 0.28), the test set (n = 88, r(2) = 0.87, RMSE = 0.39, s = 0.25) and the outlier set (n = 11, r(2) = 0.96, RMSE = 0.10, s = 0.05). The generated PhE/SVM model also showed high accuracy when subjected to those validation criteria generally adopted to gauge the predictivity of a theoretical model.

CONCLUSIONS/SIGNIFICANCE

This accurate, fast and robust PhE/SVM model that can take into account the promiscuous nature of P-gp can be applied to predict the P-gp inhibition of structurally diverse compounds that otherwise cannot be done by any other methods in a high-throughput fashion to facilitate drug discovery and development by designing drug candidates with better metabolism profile.

ATP hydrolysis at one of the two sites in ABC transporters initiates transport related conformational transitions

ABC transporters play important roles in all types of organisms by participating in physiological and pathological processes. In order to modulate the function of ABC transporters, detailed knowledge regarding their structure and dynamics is necessary. Available structures of ABC proteins indicate three major conformations, a nucleotide-bound "bottom-closed" state with the two nucleotide binding domains (NBDs) tightly closed, and two nucleotide-free conformations, the "bottom-closed" and the "bottom-open", which differ in the extent of separation of the NBDs. However, it remains a question how the widely open conformation should be interpreted, and whether hydrolysis at one of the sites can drive conformational transitions while the NBDs remain in contact. To extend our knowledge, we have investigated the dynamic properties of the Sav1866 transporter using molecular dynamics (MD) simulations. We demonstrate that the replacement of one ATP by ADP alters the correlated motion patterns of the NBDs and the transmembrane domains (TMD). The results suggest that the hydrolysis of a single nucleotide could lead to extracellular closure, driving the transport cycle. Essential dynamics analysis of simulations suggests that single nucleotide hydrolysis can drive the system toward a "bottom-closed" apo conformation similar to that observed in the structure of the MsbA transporter. We also found significant structural instability of the "bottom-open" form of the transporters in simulations. Our results suggest that ATP hydrolysis at one of the sites promotes transport related conformational changes leading to the "bottom-closed" apo conformation, which could thus be physiologically more relevant for describing the structure of the apo state.

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.

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.

Molecular dissection of dual pseudosymmetric solute translocation pathways in human P-glycoprotein

The human multispecific drug efflux transporter P-glycoprotein (P-gp) causes drug resistance and modulates the pharmacological profile of systemically administered medicines. It has arisen from a homodimeric ancestor by gene duplication. Crystal structures of mouse MDR1A indicate that P-gp shares the overall architecture with two homodimeric bacterial exporters, Sav1866 and MsbA, which have complete rotational symmetry. For ATP-binding cassette transporters, nucleotide binding occurs in two symmetric positions in the motor domains. Based on the homology with entirely symmetric half-transporters, the present study addressed the key question: can biochemical evidence for the existence of dual drug translocation pathways in the transmembrane domains of P-gp be found? P-gp was photolabeled with propafenone analogs, purified, and digested proteolytically, and peptide fragments were identified by high-resolution mass spectrometry. Labeling was assigned to two regions in the protein by projecting data into homology models. Subsequently, symmetric residue pairs in the putative translocation pathways were identified and replaced by site-directed mutagenesis. Transport assays corroborated the existence of two pseudosymmetric translocation pathways. Although rhodamine123 has a preference to take one path, verapamil, propafenones, and vinblastine preferentially use the other. Two major findings ensued from this study: the existence of two solute translocation pathways in P-gp as a reflection of evolutionary origin from a homodimeric ancestor and selective but not exclusive use of one of these pathways by different P-gp solutes. The pseudosymmetric behavior reconciles earlier kinetic and thermodynamic data, suggesting an alternative concept of drug transport by P-gp that will aid in understanding the off-target quantitative structure activity relationships of P-gp interacting drugs.

Human P-glycoprotein is active when the two halves are clamped together in the closed conformation

The P-glycoprotein (P-gp, ABCB1) drug pump protects us from toxic compounds and confers multidrug resistance. Each of the two homologous halves of P-gp is composed of a transmembrane domain (TMD) with six TM segments followed by a nucleotide-binding domain (NBD). The drug- and ATP-binding sites reside at the interface between the TMDs and NBDs, respectively. Crystal structures show drug pumps in the open and closed conformations, where the drug-binding pocket and NBDs are open or closed at the cytoplasmic side, respectively. Although it has been postulated that drug substrates enter the drug-binding pocket in the open conformation, it is unknown if they can enter in the closed conformation. To determine this, we introduced cysteines into regions of TM3 (residues 175-178) and TM9 (residues 820-822) that extend into the cytoplasm and are 4 A and 20 A apart in the closed and open conformations, respectively. The 12 double cysteine mutants were then cross-linked with a short cross-linker, M1M (4 A) at 0 degrees C to reduce thermal motion in the protein. Only mutant L175C/N820C was cross-linked. Cross-linking was not increased in the presence of ATP or drug substrates. Cross-linking increased its basal ATPase activity about 3-fold. Activity could be increased further by drug substrates such as verapamil and rhodamine B. These results suggest that P-gp in the membrane is in the closed conformation that has a high affinity for drug substrates.

Impact of the Recent Mouse P-Glycoprotein Structure for Structure-Based Ligand Design

P-Glycoprotein (P-gp), a transmembrane, ATP-dependent drug efflux transporter, has attracted considerable interest both with respect to its role in tumour cell multidrug resistance and in absorption-distribution and elimination of drugs. Although known since more than 30 years, the understanding of the molecular basis of drug/transporter interaction is still limited, which is mainly due to the lack of structural information available. However, within the past decade X-ray structures of several bacterial homologues as well as very recently also of mouse P-gp have become available. Within this review we give an overview on the current status of structural information available and on its impact for structure-based drug design.

Predicting ligand interactions with ABC transporters in ADME

ABC-type drug efflux pumps, e.g., ABCB1 (=P-glycoprotein, =MDR1), ABCC1 (=MRP1), and ABCG2 (=MXR, =BCRP), confer a multi-drug resistance (MDR) phenotype to cancer cells. Furthermore, the important contribution of ABC transporters for bioavailability, distribution, elimination, and blood-brain barrier permeation of drug candidates is increasingly recognized. This review presents an overview on the different computational methods and models pursued to predict ABC transporter substrate properties of drug-like compounds. They encompass ligand-based approaches ranging from 'simple rule'-based efforts to sophisticated machine learning methods. Many of these models show excellent performance for the data sets used. However, due to the complex nature of the applied methods, useful interpretation of the models that can be directly translated into chemical structures by the medicinal chemist is rather difficult. Additionally, very recent and promising attempts in the field of structure-based modeling of ABC transporters, which embody homology modeling as well as recently published X-ray structures of murine ABCB1, will be discussed.

Identification of residues in the drug translocation pathway of the human multidrug resistance P-glycoprotein by arginine mutagenesis

P-glycoprotein (P-gp, ATP-binding cassette B1) is a drug pump that extracts toxic drug substrates from the plasma membrane and catalyzes their ATP-dependent efflux. To map the residues in the drug translocation pathway, we performed arginine-scanning mutagenesis on all transmembrane (TM) segments (total = 237 residues) of a P-gp processing mutant (G251V) defective in folding (15% maturation efficiency) (glycosylation state used to monitor folding). The rationale was that arginines introduced into the drug-binding sites would mimic drug rescue and enhance maturation of wild-type or processing mutants of P-gp. It was found that 38 of the 89 mutants that matured had enhanced maturation. Enhancer mutations were found in 11 of the 12 TM segments with the largest number found in TMs 6 and 12 (seven in each), TMs that are critical for P-gp-drug substrate interactions. Modeling of the TM segments showed that the enhancer arginines were found on the hydrophilic face, whereas inhibitory arginines were located on a hydrophobic face that may be in contact with the lipid bilayer. It was found that many of the enhancer arginines caused large alterations in P-gp-drug interactions in ATPase assays. For example, mutants A302R (TM5), L339R (TM6), G872R (TM10), F942R (TM11), Q946R (TM11), V982R (TM12), and S993R (TM12) reduced the apparent affinity for verapamil by approximately 10-fold, whereas the F336R (TM6) and M986R (TM12) mutations caused at least a 10-fold increase in apparent affinity for rhodamine B. The results suggest that P-gp contains a large aqueous-filled drug translocation pathway with multiple drug-binding sites that can accommodate the bulky arginine side chains to promote folding of the protein.