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

The ABC transporter MsbA adopts the wide inward-open conformation in E. coli cells

Membrane proteins are currently investigated after detergent extraction from native cellular membranes and reconstitution into artificial liposomes or nanodiscs, thereby removing them from their physiological environment. However, to truly understand the biophysical properties of membrane proteins in a physiological environment, they must be investigated within living cells. Here, we used a spin-labeled nanobody to interrogate the conformational cycle of the ABC transporter MsbA by double electron-electron resonance. Unexpectedly, the wide inward-open conformation of MsbA, commonly considered a nonphysiological state, was found to be prominently populated in Escherichia coli cells. Molecular dynamics simulations revealed that extensive lateral portal opening is essential to provide access of its large natural substrate core lipid A to the binding cavity. Our work paves the way to investigate the conformational landscape of membrane proteins in cells.

Is the emperor wearing shorts? The published structures of ABC transporters

ABC transporters use the energy of ATP binding and hydrolysis to transport substrates across cellular membranes. They comprise two highly conserved nucleotide binding domains and two transmembrane domains that form the transmembrane channel and contain the substrate binding sites. Structural analyses have found a variety of seemingly unrelated folds for the ABC transporter transmembrane domains, and from these, a set of diverse mechanistic models has been inferred. Nevertheless, in spite of the explosion in structure determination of ABC transporters in the last decade, advancement in certainty and clarity as to fundamental aspects of their molecular mechanisms remains elusive. With this in mind, here we put and examine the question: Could current ABC structures differ from the physiologic membrane-embedded forms?

Theoretical insights on helix repacking as the origin of P-glycoprotein promiscuity

P-glycoprotein (P-gp, ABCB1) overexpression is, currently, one of the most important multidrug resistance (MDR) mechanisms in tumor cells. Thus, modulating drug efflux by P-gp has become one of the most promising approaches to overcome MDR in cancer. Yet, more insights on the molecular basis of drug specificity and efflux-related signal transmission mechanism between the transmembrane domains (TMDs) and the nucleotide binding domains (NBDs) are needed to develop molecules with higher selectivity and efficacy. Starting from a murine P-gp crystallographic structure at the inward-facing conformation (PDB ID: 4Q9H), we evaluated the structural quality of the herein generated human P-gp homology model. This initial human P-gp model, in the presence of the “linker” and inserted in a suitable lipid bilayer, was refined through molecular dynamics simulations and thoroughly validated. The best human P-gp model was further used to study the effect of four single-point mutations located at the TMDs, experimentally related with changes in substrate specificity and drug-stimulated ATPase activity. Remarkably, each P-gp mutation is able to induce transmembrane α-helices (TMHs) repacking, affecting the drug-binding pocket volume and the drug-binding sites properties (e.g. volume, shape and polarity) finally compromising drug binding at the substrate binding sites. Furthermore, intracellular coupling helices (ICH) also play an important role since changes in the TMHs rearrangement are shown to have an impact in residue interactions at the ICH-NBD interfaces, suggesting that identified TMHs repacking affect TMD-NBD contacts and interfere with signal transmission from the TMDs to the NBDs.

Replacing the eleven native tryptophans by directed evolution produces an active P-glycoprotein with site-specific, non-conservative substitutions

P-glycoprotein (Pgp) pumps an array of hydrophobic compounds out of cells, and has major roles in drug pharmacokinetics and cancer multidrug resistance. Yet, polyspecific drug binding and ATP hydrolysis-driven drug export in Pgp are poorly understood. Fluorescence spectroscopy using tryptophans (Trp) inserted at strategic positions is an important tool to study ligand binding. In Pgp, this method will require removal of 11 endogenous Trps, including highly conserved Trps that may be important for function, protein-lipid interactions, and/or protein stability. Here, we developed a directed evolutionary approach to first replace all eight transmembrane Trps and select for transport-active mutants in Saccharomyces cerevisiae. Surprisingly, many Trp positions contained non-conservative substitutions that supported in vivo activity, and were preferred over aromatic amino acids. The most active construct, W(3Cyto), served for directed evolution of the three cytoplasmic Trps, where two positions revealed strong functional bias towards tyrosine. W(3Cyto) and Trp-less Pgp retained wild-type-like protein expression, localization and transport function, and purified proteins retained drug stimulation of ATP hydrolysis and drug binding affinities. The data indicate preferred Trp substitutions specific to the local context, often dictated by protein structural requirements and/or membrane lipid interactions, and these new insights will offer guidance for membrane protein engineering.

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.

Comparison of mechanistic transport cycle models of ABC exporters

ABC (ATP binding cassette) transporters, ubiquitous in all kingdoms of life, carry out essential substrate transport reactions across cell membranes. Their transmembrane domains bind and translocate substrates and are connected to a pair of nucleotide binding domains, which bind and hydrolyze ATP to energize import or export of substrates. Over four decades of investigations into ABC transporters have revealed numerous details from atomic-level structural insights to their functional and physiological roles. Despite all these advances, a comprehensive understanding of the mechanistic principles of ABC transporter function remains elusive. The human multidrug resistance transporter ABCB1, also referred to as P-glycoprotein (P-gp), is one of the most intensively studied ABC exporters. Using ABCB1 as the reference point, we aim to compare the dominating mechanistic models of substrate transport and ATP hydrolysis for ABC exporters and to highlight the experimental and computational evidence in their support. In particular, we point out in silico studies that enhance and complement available biochemical data. “This article is part of a Special Issue entitled: Beyond the Structure Function Horizon of Membrane Proteins edited by Ute Hellmich, Rupak Doshi and Benjamin McIlwain.”

Substrate-induced conformational changes in the nucleotide-binding domains of lipid bilayer-associated P-glycoprotein during ATP hydrolysis

P-glycoprotein (Pgp) is an efflux pump important in multidrug resistance of cancer cells and in determining drug pharmacokinetics. Pgp is a prototype ATP-binding cassette transporter with two nucleotide-binding domains (NBDs) that bind and hydrolyze ATP. Conformational changes at the NBDs (the Pgp engines) lead to changes across Pgp transmembrane domains that result in substrate translocation. According to current alternating access models (substrate-binding pocket accessible only to one side of the membrane at a time), binding of ATP promotes NBD dimerization, resulting in external accessibility of the drug-binding site (outward-facing, closed NBD conformation), and ATP hydrolysis leads to dissociation of the NBDs with the subsequent return of the accessibility of the binding site to the cytoplasmic side (inward-facing, open NBD conformation). However, previous work has not investigated these events under near-physiological conditions in a lipid bilayer and in the presence of transport substrate. Here, we used luminescence resonance energy transfer (LRET) to measure the distances between the two Pgp NBDs. Pgp was labeled with LRET probes, reconstituted in lipid nanodiscs, and the distance between the NBDs was measured at 37 °C. In the presence of verapamil, a substrate that activates ATP hydrolysis, the NBDs of Pgp reconstituted in nanodiscs were never far apart during the hydrolysis cycle, and we never observed the NBD-NBD distances of tens of Å that have previously been reported. However, we found two main conformations that coexist in a dynamic equilibrium under all conditions studied. Our observations highlight the importance of performing studies of efflux pumps under near-physiological conditions, in a lipid bilayer, at 37 °C, and during substrate-stimulated hydrolysis.

Predictive Structure and Topology of Peroxisomal ATP-Binding Cassette (ABC) Transporters

The peroxisomal ATP-binding Cassette (ABC) transporters, which are called ABCD1, ABCD2 and ABCD3, are transmembrane proteins involved in the transport of various lipids that allow their degradation inside the organelle. Defective ABCD1 leads to the accumulation of very long-chain fatty acids and is associated with a complex and severe neurodegenerative disorder called X-linked adrenoleukodystrophy (X-ALD). Although the nucleotide-binding domain is highly conserved and characterized within the ABC transporters family, solid data are missing for the transmembrane domain (TMD) of ABCD proteins. The lack of a clear consensus on the secondary and tertiary structure of the TMDs weakens any structure-function hypothesis based on the very diverse ABCD1 mutations found in X-ALD patients. Therefore, we first reinvestigated thoroughly the structure-function data available and performed refined alignments of ABCD protein sequences. Based on the 2.85  Å resolution crystal structure of the mitochondrial ABC transporter ABCB10, here we propose a structural model of peroxisomal ABCD proteins that specifies the position of the transmembrane and coupling helices, and highlight functional motifs and putative important amino acid residues.

Mechanistic diversity in ATP-binding cassette (ABC) transporters

ABC transporters catalyze transport reactions, such as the high-affinity uptake of micronutrients into bacteria and the export of cytotoxic compounds from mammalian cells. Crystal structures of ABC domains and full transporters have provided a framework for formulating reaction mechanisms of ATP-driven substrate transport, but recent studies have suggested remarkable mechanistic diversity within this protein family. This review evaluates the differing mechanistic proposals and outlines future directions for the exploration of ABC-transporter-catalyzed reactions.

How Intrinsic Dynamics Mediates the Allosteric Mechanism in the ABC Transporter Nucleotide Binding Domain Dimer

A protein's architecture facilitates specific motions-intrinsic dynamic modes-that are employed to effect function. Here we used molecular dynamics (MD) simulations to investigate the dynamics of the MJ0796 ABC transporter nucleotide-binding domain (NBD). ABC transporter NBDs form a rotationally symmetric dimer whereby two equivalent active sites are formed at their interface; in complex with a dimer of transmembrane domains they hydrolyze ATP to energize translocation of substrates across cellular membranes. Our data suggest the ABC NBD's ensemble of functional states can be understood predominately in terms of conformational changes between its major subdomains, occurring along two orthogonal dynamic modes. The data show that ligands and oligomeric interactions modulate the equilibrium conformation of the NBD with respect to these motions, suggesting that allostery is achieved by affecting the energetic profile along these two modes. The observed dynamics and allostery integrate consonantly and logically within a mechanistic framework for the ABC NBD dimer, which is supported by a large body of experimental and theoretical data, providing a higher resolution view of the enzyme's dynamic cycle. Our study shows how valuable mechanistic inferences can be derived from accessible short-time scale MD simulations of an enzyme's substructures.

P-glycoprotein ATPase activity requires lipids to activate a switch at the first transmission interface

P-glycoprotein (P-gp) is an ABC (ATP-Binding Cassette) drug pump. A common feature of ABC proteins is that they are organized into two wings. Each wing contains a transmembrane domain (TMD) and a nucleotide-binding domain (NBD). Drug substrates and ATP bind at the interface between the TMDs and NBDs, respectively. Drug transport involves ATP-dependent conformational changes between inward- (open, NBDs far apart) and outward-facing (closed, NBDs close together) conformations. P-gps crystallized in the presence of detergent show an open structure. Human P-gp is inactive in detergent but basal ATPase activity is restored upon addition of lipids. The lipids might cause closure of the wings to bring the NBDs close together to allow ATP hydrolysis. We show however, that cross-linking the wings together did not activate ATPase activity when lipids were absent suggesting that lipids may induce other structural changes required for ATPase activity. We then tested the effect of lipids on disulfide cross-linking of mutants at the first transmission interface between intracellular loop 4 (TMD2) and NBD1. Mutants L443C/S909C and L443C/R905C but not G471C/S909C and V472C/S909C were cross-linked with oxidant when in membranes. The mutants were then purified and cross-linked with or without lipids. Mutants G471C/S909C and V472C/S909C cross-linked only in the absence of lipids whereas mutants L443C/S909C and L443C/R905C were cross-linked only in the presence of lipids. The results suggest that lipids activate a switch at the first transmission interface and that the structure of P-gp is different in detergents and lipids.

Mapping the Binding Site of the Inhibitor Tariquidar That Stabilizes the First Transmembrane Domain of P-glycoprotein

ABC (ATP-binding cassette) transporters are clinically important because drug pumps like P-glycoprotein (P-gp, ABCB1) confer multidrug resistance and mutant ABC proteins are responsible for many protein-folding diseases such as cystic fibrosis. Identification of the tariquidar-binding site has been the subject of intensive molecular modeling studies because it is the most potent inhibitor and corrector of P-gp. Tariquidar is a unique P-gp inhibitor because it locks the pump in a conformation that blocks drug efflux but activates ATPase activity. In silico docking studies have identified several potential tariquidar-binding sites. Here, we show through cross-linking studies that tariquidar most likely binds to sites within the transmembrane (TM) segments located in one wing or at the interface between the two wings (12 TM segments form 2 divergent wings). We then introduced arginine residues at all positions in the 12 TM segments (223 mutants) of P-gp. The rationale was that a charged residue in the drug-binding pocket would disrupt hydrophobic interaction with tariquidar and inhibit its ability to rescue processing mutants or stimulate ATPase activity. Arginines introduced at 30 positions significantly inhibited tariquidar rescue of a processing mutant and activation of ATPase activity. The results suggest that tariquidar binds to a site within the drug-binding pocket at the interface between the TM segments of both structural wings. Tariquidar differed from other drug substrates, however, as it stabilized the first TM domain. Stabilization of the first TM domain appears to be a key mechanism for high efficiency rescue of ABC processing mutants that cause disease.

Structure, mechanism and cooperation of bacterial multidrug transporters

Cells from all domains of life encode energy-dependent trans-membrane transporters that can expel harmful substances including clinically applied therapeutic agents. As a collective body, these transporters perform as a super-system that confers tolerance to an enormous range of harmful compounds and consequently aid survival in hazardous environments. In the Gram-negative bacteria, some of these transporters serve as energy-transducing components of tripartite assemblies that actively efflux drugs and other harmful compounds, as well as deliver virulence agents across the entire cell envelope. We draw together recent structural and functional data to present the current models for the transport mechanisms for the main classes of multi-drug transporters and their higher-order assemblies.

In silico identified targeted inhibitors of P-glycoprotein overcome multidrug resistance in human cancer cells in culture

Failure of cancer chemotherapies is often linked to the over expression of ABC efflux transporters like the multidrug resistance P-glycoprotein (P-gp). P-gp expression in cells leads to the elimination of a variety of chemically unrelated, mostly cytotoxic compounds. Administration of chemotherapeutics during therapy frequently selects for cells that over express P-gp and are therefore capable of robustly exporting diverse compounds, including chemotherapeutics, from the cells. P-gp thus confers multidrug resistance to a majority of drugs currently available for the treatment of cancers and diseases like HIV/AIDS. The search for P-gp inhibitors for use as co-therapeutics to combat multidrug resistances has had little success to date. In a previous study (Brewer et al., Mol Pharmacol 86: 716–726, 2014), we described how ultrahigh throughput computational searches led to the identification of four drug-like molecules that specifically interfere with the energy harvesting steps of substrate transport and inhibit P-gp catalyzed ATP hydrolysis in vitro. In the present study, we demonstrate that three of these compounds reversed P-gp-mediated multidrug resistance of cultured prostate cancer cells to restore sensitivity comparable to naïve prostate cancer cells to the chemotherapeutic drug, paclitaxel. Potentiation concentrations of the inhibitors were <3 μmol/L. The inhibitors did not exhibit significant toxicity to noncancerous cells at concentrations where they reversed multidrug resistance in cancerous cells. Our results indicate that these compounds with novel mechanisms of P-gp inhibition are excellent leads for the development of co-therapeutics for the treatment of multidrug resistances.

The Nucleotide-Free State of the Multidrug Resistance ABC Transporter LmrA: Sulfhydryl Cross-Linking Supports a Constant Contact, Head-to-Tail Configuration of the Nucleotide-Binding Domains

ABC transporters are integral membrane pumps that are responsible for the import or export of a diverse range of molecules across cell membranes. ABC transporters have been implicated in many phenomena of medical importance, including cystic fibrosis and multidrug resistance in humans. The molecular architecture of ABC transporters comprises two transmembrane domains and two ATP-binding cassettes, or nucleotide-binding domains (NBDs), which are highly conserved and contain motifs that are crucial to ATP binding and hydrolysis. Despite the improved clarity of recent structural, biophysical, and biochemical data, the seemingly simple process of ATP binding and hydrolysis remains controversial, with a major unresolved issue being whether the NBD protomers separate during the catalytic cycle. Here chemical cross-linking data is presented for the bacterial ABC multidrug resistance (MDR) transporter LmrA. These indicate that in the absence of nucleotide or substrate, the NBDs come into contact to a significant extent, even at 4°C, where ATPase activity is abrogated. The data are clearly not in accord with an inward-closed conformation akin to that observed in a crystal structure of V. cholerae MsbA. Rather, they suggest a head-to-tail configuration ‘sandwich’ dimer similar to that observed in crystal structures of nucleotide-bound ABC NBDs. We argue the data are more readily reconciled with the notion that the NBDs are in proximity while undergoing intra-domain motions, than with an NBD ‘Switch’ mechanism in which the NBD monomers separate in between ATP hydrolysis cycles.

The Transmission Interfaces Contribute Asymmetrically to the Assembly and Activity of Human P-glycoprotein

P-glycoprotein (P-gp; ABCB1) is an ABC drug pump that protects us from toxic compounds. It is clinically important because it confers multidrug resistance. The homologous halves of P-gp each contain a transmembrane (TM) domain (TMD) with 6 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. Each NBD is connected to the TMDs by a transmission interface involving a pair of intracellular loops (ICLs) that form ball-and-socket joints. P-gp is different from CFTR (ABCC7) in that deleting NBD2 causes misprocessing of only P-gp. Therefore, NBD2 might be critical for stabilizing ICLs 2 and 3 that form a tetrahelix bundle at the NBD2 interface. Here we report that the NBD1 and NBD2 transmission interfaces in P-gp are asymmetric. Point mutations to 25 of 60 ICL2/ICL3 residues at the NBD2 transmission interface severely reduced P-gp assembly while changes to the equivalent residues in ICL1/ICL4 at the NBD1 interface had little effect. The hydrophobic nature at the transmission interfaces was also different. Mutation of Phe-1086 or Tyr-1087 to arginine at the NBD2 socket blocked activity or assembly while the equivalent mutations at the NBD1 socket had only modest effects. The results suggest that the NBD transmission interfaces are asymmetric. In contrast to the ICL2/3-NBD2 interface, the ICL1/4-NBD1 transmission interface is more hydrophilic and insensitive to mutations. Therefore the ICL2/3-NBD2 transmission interface forms a precise hydrophobic connection that acts as a linchpin for assembly and trafficking of P-gp.

Distinct conformational spectrum of homologous multidrug ABC transporters

ATP-binding cassette (ABC) exporters are ubiquitously found in all kingdoms of life and their members play significant roles in mediating drug pharmacokinetics and multidrug resistance in the clinic. Significant questions and controversies remain regarding the relevance of their conformations observed in X-ray structures, their structural dynamics, and mechanism of transport. Here, we used single particle electron microscopy (EM) to delineate the entire conformational spectrum of two homologous ABC exporters (bacterial MsbA and mammalian P-glycoprotein) and the influence of nucleotide and substrate binding. Newly developed amphiphiles in complex with lipids that support high protein stability and activity enabled EM visualization of individual complexes in a membrane-mimicking environment. The data provide a comprehensive view of the conformational flexibility of these ABC exporters under various states and demonstrate not only similarities but striking differences between their mechanistic and energetic regulation of conformational changes.

Equilibrated atomic models of outward-facing P-glycoprotein and effect of ATP binding on structural dynamics

P-glycoprotein (Pgp) is an ATP-binding cassette (ABC) transporter that alternates between inward- and outward-facing conformations to capture and force substrates out of cells like a peristaltic pump. The high degree of similarity in outward-facing structures across evolution of ABC transporters allowed construction of a high-confidence outward-facing Pgp atomic model based on crystal structures of outward-facing Sav1866 and inward-facing Pgp. The model adhered to previous experimentally determined secondary- and tertiary- configurations during all-atom molecular dynamics simulations in the presence or absence of MgATP. Three long lasting (>100 ns) meta-stable states were apparent in the presence of MgATP revealing new insights into alternating access. The two ATP-binding pockets are highly asymmetric resulting in differential control of overall structural dynamics and allosteric regulation of the drug-binding pocket. Equilibrated Pgp has a considerably different electrostatic profile compared to Sav1866 that implicates significant kinetic and thermodynamic differences in transport mechanisms.

Molecular basis of the polyspecificity of P-glycoprotein (ABCB1): recent biochemical and structural studies

ABCB1 (P-glycoprotein/P-gp) is an ATP-binding cassette transporter well known for its association with multidrug resistance in cancer cells. Powered by the hydrolysis of ATP, it effluxes structurally diverse compounds. In this chapter, we discuss current views on the molecular basis of the substrate polyspecificity of P-gp. One of the features that accounts for this property is the structural flexibility observed in P-gp. Several X-ray crystal structures of mouse P-gp have been published recently in the absence of nucleotide, with and without bound inhibitors. All the structures are in an inward-facing conformation exhibiting different degrees of domain separation, thus revealing a highly flexible protein. Biochemical and biophysical studies also demonstrate this flexibility in mouse as well as human P-gp. Site-directed mutagenesis has revealed the existence of multiple transport-active binding sites in P-gp for a single substrate. Thus, drugs can bind at either primary or secondary sites. Biochemical, molecular modeling, and structure-activity relationship studies suggest a large, common drug-binding pocket with overlapping sites for different substrates. We propose that in addition to the structural flexibility, the molecular or chemical flexibility also contributes to the binding of substrates to multiple sites forming the basis of polyspecificity.

Tariquidar inhibits P-glycoprotein drug efflux but activates ATPase activity by blocking transition to an open conformation

P-glycoprotein (P-gp, ABCB1) is a drug pump that confers multidrug resistance. Inhibition of P-gp would improve chemotherapy. Tariquidar is a potent P-gp inhibitor but its mechanism is unknown. Here, we tested our prediction that tariquidar inhibits P-gp cycling between the open and closed states during the catalytic cycle. Transition of P-gp to an open state can be monitored in intact cells using reporter cysteines introduced into extracellular loops 1 (A80C) and 4 (R741C). Residues A80C/R741C come close enough (<7Å) to spontaneously cross-link in the open conformation (<7Å) but are widely separated (>30Å) in the closed conformation. Cross-linking of A80C/R741C can be readily detected because it causes the mutant protein to migrate slower on SDS-PAGE gels. We tested whether drug substrates or inhibitors could inhibit cross-linking of the mutant. It was found that only tariquidar blocked A80C/R741C cross-linking. Tariquidar was also a more potent pharmacological chaperone than other P-gp substrates/modulators such as cyclosporine A. Only tariquidar promoted maturation of misprocessed mutant F804D to yield mature P-gp. Tariquidar interacted with the transmembrane domains because it could rescue a misprocessed truncation mutant lacking the nucleotide-binding domains. These results show that tariquidar is a potent pharmacological chaperone and inhibits P-gp drug efflux by blocking transition to the open state during the catalytic cycle.

Cysteines introduced into extracellular loops 1 and 4 of human P-glycoprotein that are close only in the open conformation spontaneously form a disulfide bond that inhibits drug efflux and ATPase activity

P-glycoprotein (P-gp) is an ATP-binding cassette drug pump that protects us from toxic compounds and confers multidrug resistance. The protein is organized into two halves. The halves contain a transmembrane domain (TMD) with six transmembrane segments and a nucleotide-binding domain (NBD). The drug- and ATP-binding sites reside at the TMD1/TMD2 and NBD1/NBD2 interfaces, respectively. ATP-dependent drug efflux involves changes between the open inward-facing (NBDs apart, extracellular loops (ECLs) close together) and the closed outward-facing (NBDs close together, ECLs apart) conformations. It is controversial, however, whether the open conformation only exists transiently in intact cells because of the presence of high levels of ATP. To test for the presence of an open conformation in intact cells, reporter cysteines were placed in extracellular loops 1 (A80C, N half) and 4 (R741C, C half). The rationale was that cysteines A80C/R741C would only come close enough to form a disulfide bond in an open conformation (6.9 Å apart) because they are separated widely (30.4 Å apart) in the closed conformation. It was observed that the mutant A80C/R741C cross-linked spontaneously (>90%) when expressed in cells. In contrast to previous reports showing that trapping P-gp in a closed conformation highly activated ATPase activity, here we show that A80C/R741C cross-linking inhibited ATPase activity and drug efflux. Both activities were restored when the cross-linked mutant was treated with a thiol-reducing agent. The results show that an open conformation can be readily detected in cells and that cross-linking of cysteines placed in ECLs 1 and 4 inhibits activity.

Structural characterization of two metastable ATP-bound states of P-glycoprotein

ATP Binding Cassette (ABC) transporters couple the binding and hydrolysis of ATP to the transport of substrate molecules across the membrane. The mechanism by which ATP binding and/or hydrolysis drives the conformational changes associated with substrate transport has not yet been characterized fully. Here, changes in the conformation of the ABC export protein P-glycoprotein on ATP binding are examined in a series of molecular dynamics simulations. When one molecule of ATP is placed at the ATP binding site associated with each of the two nucleotide binding domains (NBDs), the membrane-embedded P-glycoprotein crystal structure adopts two distinct metastable conformations. In one, each ATP molecule interacts primarily with the Walker A motif of the corresponding NBD. In the other, the ATP molecules interacts with both Walker A motif of one NBD and the Signature motif of the opposite NBD inducing the partial dimerization of the NBDs. This interaction is more extensive in one of the two ATP binding site, leading to an asymmetric structure. The overall conformation of the transmembrane domains is not altered in either of these metastable states, indicating that the conformational changes associated with ATP binding observed in the simulations in the absence of substrate do not lead to the outward-facing conformation and thus would be insufficient in themselves to drive transport. Nevertheless, the metastable intermediate ATP-bound conformations observed are compatible with a wide range of experimental cross-linking data demonstrating the simulations do capture physiologically important conformations. Analysis of the interaction between ATP and its cofactor Mg(2+) with each NBD indicates that the coordination of ATP and Mg(2+) differs between the two NBDs. The role structural asymmetry may play in ATP binding and hydrolysis is discussed. Furthermore, we demonstrate that our results are not heavily influenced by the crystal structure chosen for initiation of the simulations.

Identification of the distance between the homologous halves of P-glycoprotein that triggers the high/low ATPase activity switch

P-glycoprotein (P-gp, ABCB1) is an ATP-binding cassette drug pump that protects us from toxic compounds and confers multidrug resistance. Each homologous half contains a transmembrane domain with six transmembrane segments followed by a nucleotide-binding domain (NBD). The drug- and ATP-binding sites reside at the interface between the transmembrane domain and NBDs, respectively. Drug binding activates ATPase activity by an unknown mechanism. There is no high resolution structure of human P-gp, but homology models based on the crystal structures of bacterial, mouse, and Caenorhabditis elegans ATP-binding cassette drug pumps yield both open (NBDs apart) and closed (NBDs together) conformations. Molecular dynamics simulations predict that the NBDs can be separated over a range of distances (over 20 Å). To determine the distance that show high or low ATPase activity, we cross-linked reporter cysteines L175C (N-half) and N820C (C-half) with cross-linkers of various lengths that separated the halves between 6 and 30 Å (α-carbons). We observed that ATPase activity increased over 10-fold when the cysteines were cross-linked at distances between 6 and 19 Å, although cross-linking at distances greater than 20 Å yielded basal levels of activity. The results suggest that the ATPase activation switch appears to be turned on or off when L175C/N820 are clamped at distances less than or greater than 20 Å, respectively. We predict that the high/low ATPase activity switch may occur at a distance where the NBDs are predicted in molecular dynamic simulations to undergo pronounced twisting as they approach each other (Wise, J. G. (2012) Biochemistry 51, 5125-5141).

Locking intracellular helices 2 and 3 together inactivates human P-glycoprotein

The P-glycoprotein (P-gp) drug pump (ABCB1) has two transmembrane domains and two nucleotide-binding domains (NBDs). Coupling of the drug-binding sites in the transmembrane domains to the NBDs occurs through interaction of the intracellular helices (IHs) with residues in the NBDs (IH1/IH4/NBD1 and IH2/IH3/NBD2). We showed previously that cross-linking of cysteines in IH3 and IH1 with a short cross-linker mimicked drug binding as it activated P-gp ATPase activity. Here we show that residue A259C(IH2) could be directly cross-linked to W803C(IH3). Cross-linking was inhibited by the presence of ATP and adenosine 5'-(β,γ-imino)triphosphate but not by ADP. Cross-linking of mutant A259C/W803C inhibited its verapamil-stimulated ATPase activity mutant, but activity was restored after addition of dithiothreitol. Because these residues are close to the ball-and-socket joint A266C(IH2)/Phe(1086)(NBD2), we mutated the adjacent Tyr(1087)(NBD2) close to IH3. Mutants Y1087A and Y1087L, but not Y1087F, were misprocessed, and all inhibited ATPase activity. Mutation of hydrophobic residues (F793A, L797A, L814A, and L818A) flanking IH3 also inhibited maturation. The results suggest that these residues, together with Trp(803) and Phe(804), form a large hydrophobic pocket. The results show that there is an important hydrophobic network at the IH2/IH3/NBD2 transmission interface that is critical for folding and activity of P-gp.

Refined structures of mouse P-glycoprotein

The recently determined C. elegans P-glycoprotein (Pgp) structure revealed significant deviations compared to the original mouse Pgp structure, which suggested possible misinterpretations in the latter model. To address this concern, we generated an experimental electron density map from single-wavelength anomalous dispersion phasing of an original mouse Pgp dataset to 3.8 Å resolution. The map exhibited significantly more detail compared to the original MAD map and revealed several regions of the structure that required de novo model building. The improved drug-free structure was refined to 3.8 Å resolution with a 9.4 and 8.1% decrease in R_work and R_free, respectively, (R_work = 21.2%, R_free = 26.6%) and a significant improvement in protein geometry. The improved mouse Pgp model contains ∼95% of residues in the favorable Ramachandran region compared to only 57% for the original model. The registry of six transmembrane helices was corrected, revealing amino acid residues involved in drug binding that were previously unrecognized. Registry shifts (rotations and translations) for three transmembrane (TM)4 and TM5 and the addition of three N-terminal residues were necessary, and were validated with new mercury labeling and anomalous Fourier density. The corrected position of TM4, which forms the frame of a portal for drug entry, had backbone atoms shifted >6 Å from their original positions. The drug translocation pathway of mouse Pgp is 96% identical to human Pgp and is enriched in aromatic residues that likely play a collective role in allowing a high degree of polyspecific substrate recognition.

Microparticle drug sequestration provides a parallel pathway in the acquisition of cancer drug resistance

Expanding on our previous findings demonstrating that microparticles (MPs) spread cancer multidrug resistance, we now show that MPs sequester drugs, reducing the free drug concentration available to cells. MPs were isolated from drug-sensitive and drug-resistant sub-clones of a human breast adenocarcinoma cell line and from human acute lymphoblastic leukemia cells. MPs were assessed for size, mitochondria, RNA and phospholipid content, P-glycoprotein (P-gp) expression and orientation and ATPase activity relative to drug sequestration capacity. Of the drug classes examined, MPs sequestered the anthracycline class to a significant degree. The degree of sequestration was likely due to the size of MPs and thus the amount of cargo they contain, to which the anthracyclines bind. Moreover, a proportion of the P-gp present on MPs was inside-out in orientation, enabling it to influx drugs rather than its typical efflux function. This was confirmed by surface immunofluorescence and by assessment of drug-stimulated ATPase activity following MP permeabilization. Thus we determined that breast cancer MPs carried a proportion of their P-gp oriented inside-out, providing active sequestration within the microvesicular compartment. These results demonstrate a capacity for MPs to sequester chemotherapeutic drugs, which has a predominantly active sequestration component for MPs derived from drug-resistant cells and a predominantly passive component for MPs derived from drug-sensitive cells. This reduction in available drug concentration has potential to contribute to a parallel pathway and complements that of the intercellular transfer of P-gp. These findings lend further support to the role of MPs in limiting the successful management of cancer.

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.

On the origin of large flexibility of P-glycoprotein in the inward-facing state

P-glycoprotein (Pgp) is one of the most biomedically relevant transporters in the ATP binding cassette (ABC) superfamily due to its involvement in developing multidrug resistance in cancer cells. Employing molecular dynamics simulations and double electron-electron resonance spectroscopy, we have investigated the structural dynamics of membrane-bound Pgp in the inward-facing state and found that Pgp adopts an unexpectedly wide range of conformations, highlighted by the degree of separation between the two nucleotide-binding domains (NBDs). The distance between the two NBDs in the equilibrium simulations covers a range of at least 20 Å, including, both, more open and more closed NBD configurations than the crystal structure. The double electron-electron resonance measurements on spin-labeled Pgp mutants also show wide distributions covering both longer and shorter distances than those observed in the crystal structure. Based on structural and sequence analyses, we propose that the transmembrane domains of Pgp might be more flexible than other structurally known ABC exporters. The structural flexibility of Pgp demonstrated here is not only in close agreement with, but also helps rationalize, the reported high NBD fluctuations in several ABC exporters and possibly represents a fundamental difference in the transport mechanism between ABC exporters and ABC importers. In addition, during the simulations we have captured partial entrance of a lipid molecule from the bilayer into the lumen of Pgp, reaching the putative drug binding site. The location of the protruding lipid suggests a putative pathway for direct drug recruitment from the membrane.

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

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

An asymmetric post-hydrolysis state of the ABC transporter ATPase dimer

ABC transporters are a superfamily of enzyme pumps that hydrolyse ATP in exchange for translocation of substrates across cellular membranes. Architecturally, ABC transporters are a dimer of transmembrane domains coupled to a dimer of nucleotide binding domains (NBDs): the NBD dimer contains two ATP-binding sites at the intersubunit interface. A current controversy is whether the protomers of the NBD dimer separate during ATP hydrolysis cycling, or remain in constant contact. In order to investigate the ABC ATPase catalytic mechanism, MD simulations using the recent structure of the ADP+Pi-bound MJ0796 isolated NBD dimer were performed. In three independent simulations of the ADP+Pi/apo state, comprising a total of >0.5 µs, significant opening of the apo (empty) active site was observed; occurring by way of intrasubunit rotations between the core and helical subdomains within both NBD monomers. In contrast, in three equivalent simulations of the ATP/apo state, the NBD dimer remained close to the crystal structure, and no opening of either active site occurred. The results thus showed allosteric coupling between the active sites, mediated by intrasubunit conformational changes. Opening of the apo site is exquisitely tuned to the nature of the ligand, and thus to the stage of the reaction cycle, in the opposite site. In addition to this, in also showing how one active site can open, sufficient to bind nucleotide, while the opposite site remains occluded and bound to the hydrolysis products ADP+Pi, the results are consistent with a Constant Contact Model. Conversely, they show how there may be no requirement for the NBD protomers to separate to complete the catalytic cycle.

Molecular disruption of the power stroke in the ATP-binding cassette transport protein MsbA

ATP-binding cassette transporters affect drug pharmacokinetics and are associated with inherited human diseases and impaired chemotherapeutic treatment of cancers and microbial infections. Current alternating access models for ATP-binding cassette exporter activity suggest that ATP binding at the two cytosolic nucleotide-binding domains provides a power stroke for the conformational switch of the two membrane domains from the inward-facing conformation to the outward-facing conformation. In outward-facing crystal structures of the bacterial homodimeric ATP-binding cassette transporters MsbA from gram-negative bacteria and Sav1866 from Staphylococcus aureus, two transmembrane helices (3 and 4) in the membrane domains have their cytoplasmic extensions in close proximity, forming a tetrahelix bundle interface. In biochemical experiments on MsbA from Escherichia coli, we show for the first time that a robust network of inter-monomer interactions in the tetrahelix bundle is crucial for the transmission of nucleotide-dependent conformational changes to the extracellular side of the membrane domains. Our observations are the first to suggest that modulation of tetrahelix bundle interactions in ATP-binding cassette exporters might offer a potent strategy to alter their transport activity.

Mechanism of the ABC transporter ATPase domains: catalytic models and the biochemical and biophysical record

ABC transporters comprise a large, diverse, and ubiquitous superfamily of membrane active transporters. Their core architecture is a dimer of dimers, comprising two transmembrane domains that bind substrate and form the channel, and two ATP-binding cassettes, which bind and hydrolyze ATP to energize the translocase function. The prevailing paradigm for the ABC transport mechanism is the Switch Model, in which the nucleotide binding domains are proposed to dimerise upon binding of two ATP molecules, and thence dissociate upon sequential hydrolysis of the ATP. This idea appears consistent with crystal structures of both isolated subunits and whole transporters, as well as with a significant body of biochemical data. Nonetheless, an alternative Constant Contact Model has been proposed, in which the nucleotide binding domains do not fully dissociate, and ATP hydrolysis occurs alternately at each of the two active sites. Here, we review the biochemical and biophysical data relating to the ABC catalytic mechanism, to show how they may be construed as consistent with a Constant Contact Model, and to assess to what extent they support the Switch Model.

Mammalian P4-ATPases and ABC transporters and their role in phospholipid transport

Transport of phospholipids across cell membranes plays a key role in a wide variety of biological processes. These include membrane biosynthesis, generation and maintenance of membrane asymmetry, cell and organelle shape determination, phagocytosis, vesicle trafficking, blood coagulation, lipid homeostasis, regulation of membrane protein function, apoptosis, etc. P(4)-ATPases and ATP binding cassette (ABC) transporters are the two principal classes of membrane proteins that actively transport phospholipids across cellular membranes. P(4)-ATPases utilize the energy from ATP hydrolysis to flip aminophospholipids from the exocytoplasmic (extracellular/lumen) to the cytoplasmic leaflet of cell membranes generating membrane lipid asymmetry and lipid imbalance which can induce membrane curvature. Many ABC transporters play crucial roles in lipid homeostasis by actively transporting phospholipids from the cytoplasmic to the exocytoplasmic leaflet of cell membranes or exporting phospholipids to protein acceptors or micelles. Recent studies indicate that some ABC proteins can also transport phospholipids in the opposite direction. The importance of P(4)-ATPases and ABC transporters is evident from the findings that mutations in many of these transporters are responsible for severe human genetic diseases linked to defective phospholipid transport. This article is part of a Special Issue entitled Phospholipids and Phospholipid Metabolism.

P-glycoprotein trafficking at the blood-brain barrier altered by peripheral inflammatory hyperalgesia

P-glycoprotein (ABCB1/MDR1, EC 3.6.3.44), the major efflux transporter at the blood-brain barrier (BBB), is a formidable obstacle to CNS pharmacotherapy. Understanding the mechanism(s) for increased P-glycoprotein activity at the BBB during peripheral inflammatory pain is critical in the development of novel strategies to overcome the significant decreases in CNS analgesic drug delivery. In this study, we employed the λ-carrageenan pain model (using female Sprague-Dawley rats), combined with confocal microscopy and subcellular fractionation of cerebral microvessels, to determine if increased P-glycoprotein function, following the onset of peripheral inflammatory pain, is associated with a change in P-glycoprotein trafficking which leads to pain-induced effects on analgesic drug delivery. Injection of λ-carrageenan into the rat hind paw induced a localized, inflammatory pain (hyperalgesia) and simultaneously, at the BBB, a rapid change in colocalization of P-glycoprotein with caveolin-1, a key scaffolding/trafficking protein. Subcellular fractionation of isolated cerebral microvessels revealed that the bulk of P-glycoprotein constitutively traffics to membrane domains containing high molecular weight, disulfide-bonded P-glycoprotein-containing structures that cofractionate with membrane domains enriched with monomeric and high molecular weight, disulfide-bonded, caveolin-1-containing structures. Peripheral inflammatory pain promoted a dynamic redistribution between membrane domains of P-glycoprotein and caveolin-1. Disassembly of high molecular weight P-glycoprotein-containing structures within microvascular endothelial luminal membrane domains was accompanied by an increase in ATPase activity, suggesting a potential for functionally active P-glycoprotein. These results are the first observation that peripheral inflammatory pain leads to specific structural changes in P-glycoprotein responsible for controlling analgesic drug delivery to the CNS.

Perspectives on the structure-function of ABC transporters: the Switch and Constant Contact models

ABC transporters constitute one of the largest protein families across the kingdoms of archaea, eubacteria and eukarya. They couple ATP hydrolysis to vectorial translocation of diverse substrates across membranes. The ABC transporter architecture comprises two transmembrane domains and two cytosolic ATP-binding cassettes. During 2002-2012, nine prokaryotic ABC transporter structures and two eukaryotic structures have been solved to medium resolution. Despite a wealth of biochemical, biophysical, and structural data, fundamental questions remain regarding the coupling of ATP hydrolysis to unidirectional substrate translocation, and the mechanistic suite of steps involved. The mechanics of the ATP cassette dimer is defined most popularly by the 'Switch Model', which proposes that hydrolysis in each protomer is sequential, and that as the sites are freed of nucleotide, the protomers lose contact across a large solvent-filled gap of 20-30 Å; as captured in several X-ray solved structures. Our 'Constant Contact' model for the operational mechanics of ATP binding and hydrolysis in the ATP-binding cassettes is derived from the 'alternating sites' model, proposed in 1995, and which requires an intrinsic asymmetry in the ATP sites, but does not require the partner protomers to lose contact. Thus one of the most debated issues regarding the function of ABC transporters is whether the cooperative mechanics of ATP hydrolysis requires the ATP cassettes to separate or remain in constant contact and this dilemma is discussed at length in this review.

The ATPase activity of the P-glycoprotein drug pump is highly activated when the N-terminal and central regions of the nucleotide-binding domains are linked closely together

The P-glycoprotein (P-gp, ABCB1) drug pump protects us from toxic compounds and confers multidrug resistance. Each of the homologous halves of P-gp is composed of a transmembrane domain (TMD) with 6 TM segments followed by a nucleotide-binding domain (NBD). The predicted drug- and ATP-binding sites reside at the interface between the TMDs and NBDs, respectively. Crystal structures and EM projection images suggest that the two halves of P-gp are separated by a central cavity that closes upon binding of nucleotide. Binding of drug substrates may induce further structural rearrangements because they stimulate ATPase activity. Here, we used disulfide cross-linking with short (8 Å) or long (22 Å) cross-linkers to identify domain-domain interactions that activate ATPase activity. It was found that cross-linking of cysteines that lie close to the LSGGQ (P517C) and Walker A (I1050C) sites of NBD1 and NBD2, respectively, as well as the cytoplasmic extensions of TM segments 3 (D177C or L175C) and 9 (N820C) with a short cross-linker activated ATPase activity over 10-fold. A pyrylium compound that inhibits ATPase activity blocked cross-linking at these sites. Cross-linking between the NBDs was not inhibited by tariquidar, a drug transport inhibitor that stimulates P-gp ATPase activity but is not transported. Cross-linking between extracellular cysteines (T333C/L975C) predicted to lock P-gp into a conformation that prevents close NBD association inhibited ATPase activity. The results suggest that trapping P-gp in a conformation in which the NBDs are closely associated likely mimics the structural rearrangements caused by binding of drug substrates that stimulate ATPase activity.

Catalytic transitions in the human MDR1 P-glycoprotein drug binding sites

Multidrug resistance proteins that belong to the ATP-binding cassette family like the human P-glycoprotein (ABCB1 or Pgp) are responsible for many failed cancer and antiviral chemotherapies because these membrane transporters remove the chemotherapeutics from the targeted cells. Understanding the details of the catalytic mechanism of Pgp is therefore critical to the development of inhibitors that might overcome these resistances. In this work, targeted molecular dynamics techniques were used to elucidate catalytically relevant structures of Pgp. Crystal structures of homologues in four different conformations were used as intermediate targets in the dynamics simulations. Transitions from conformations that were wide open to the cytoplasm to transition state conformations that were wide open to the extracellular space were studied. Twenty-six nonredundant transitional protein structures were identified from these targeted molecular dynamics simulations using evolutionary structure analyses. Coupled movement of nucleotide binding domains (NBDs) and transmembrane domains (TMDs) that form the drug binding cavities were observed. Pronounced twisting of the NBDs as they approached each other as well as the quantification of a dramatic opening of the TMDs to the extracellular space as the ATP hydrolysis transition state was reached were observed. Docking interactions of 21 known transport ligands or inhibitors were analyzed with each of the 26 transitional structures. Many of the docking results obtained here were validated by previously published biochemical determinations. As the ATP hydrolysis transition state was approached, drug docking in the extracellular half of the transmembrane domains seemed to be destabilized as transport ligand exit gates opened to the extracellular space.

The Effect of Environment on the Structure of a Membrane Protein: P-Glycoprotein under Physiological Conditions

The stability of the crystal structure of the multidrug transporter P-glycoprotein proposed by Aller et al. (PDBid 3G5U ) has been examined under different environmental conditions using molecular dynamics. We show that in the presence of the detergent cholate, the structure of P-glycoprotein solved at pH 7.5 is stable. However, when incorporated into a cholesterol-enriched POPC membrane in the presence of 150 mM NaCl, the structure rapidly deforms. Only when the simulation conditions closely matched the experimental conditions under which P-glycoprotein is transport active was a stable conformation obtained. Specifically, the presence of Mg(2+), which bound to distinct sites in the nucleotide binding domains (NBDs), and the double protonation of the catalytic histidines (His583 and His1228) and His149 were required. While the structure obtained in a membrane environment under these conditions is very similar to the crystal structure of Aller et al., there are several key differences. The NBDs are in direct contact, reminiscent of the open state of MalK. The angle between the transmembrane domains is also increased, resulting in an outward motion of the intracellular loops. Notably, the structures obtained from the simulations provide a better match to a range of experimental cross-linking data than does the original 3G5U-a crystal structure. This work highlights the effect small changes in environmental conditions can have of the conformation of a membrane protein and the importance of representing the experimental conditions appropriately in modeling studies.

Novel cGMP efflux inhibitors identified by virtual ligand screening (VLS) and confirmed by experimental studies

Elevated intracellular levels of cyclic guanosine monophosphate (cGMP) may induce apoptosis, and at least some cancer cells seem to escape this effect by increased efflux of cGMP, as clinical studies have shown that extracellular cGMP levels are elevated in various types of cancer. The human ATP binding cassette (ABC) transporter ABCC5 transports cGMP out of cells, and inhibition of ABCC5 may have cytotoxic effects. Sildenafil inhibits cGMP efflux by binding to ABCC5, and in order to search for potential novel ABCC5 inhibitors, we have identified sildenafil derivates using structural and computational guidance and tested them for the cGMP efflux effect. Eleven compounds from virtual ligand screening (VLS) were tested in vitro, using inside-out vesicles (IOV), for inhibition of cGMP efflux. Seven of 11 compounds predicted by VLS to bind to ABCC5 were more potent than sildenafil, and the two most potent showed K(i) of 50-100 nM.

Dynamic ligand-induced conformational rearrangements in P-glycoprotein as probed by fluorescence resonance energy transfer spectroscopy

P-glycoprotein (Pgp), a member of the ATP-binding cassette transporter family, functions as an ATP hydrolysis-driven efflux pump to rid the cell of toxic organic compounds, including a variety of drugs used in anticancer chemotherapy. Here, we used fluorescence resonance energy transfer (FRET) spectroscopy to delineate the structural rearrangements the two nucleotide binding domains (NBDs) are undergoing during the catalytic cycle. Pairs of cysteines were introduced into equivalent regions in the N- and C-terminal NBDs for labeling with fluorescent dyes for ensemble and single-molecule FRET spectroscopy. In the ensemble FRET, a decrease of the donor to acceptor (D/A) ratio was observed upon addition of drug and ATP. Vanadate trapping further decreased the D/A ratio, indicating close association of the two NBDs. One of the cysteine mutants was further analyzed using confocal single-molecule FRET spectroscopy. Single Pgp molecules showed fast fluctuations of the FRET efficiencies, indicating movements of the NBDs on a time scale of 10-100 ms. Populations of low, medium, and high FRET efficiencies were observed during drug-stimulated MgATP hydrolysis, suggesting the presence of at least three major conformations of the NBDs during catalysis. Under conditions of vanadate trapping, most molecules displayed high FRET efficiency states, whereas with cyclosporin, more molecules showed low FRET efficiency. Different dwell times of the FRET states were found for the distinct biochemical conditions, with the fastest movements during active turnover. The FRET spectroscopy observations are discussed in context of a model of the catalytic mechanism of Pgp.

Predicting P-glycoprotein-mediated drug transport based on support vector machine and three-dimensional crystal structure of P-glycoprotein

Human P-glycoprotein (P-gp) is an ATP-binding cassette multidrug transporter that confers resistance to a wide range of chemotherapeutic agents in cancer cells by active efflux of the drugs from cells. P-gp also plays a key role in limiting oral absorption and brain penetration and in facilitating biliary and renal elimination of structurally diverse drugs. Thus, identification of drugs or new molecular entities to be P-gp substrates is of vital importance for predicting the pharmacokinetics, efficacy, safety, or tissue levels of drugs or drug candidates. At present, publicly available, reliable in silico models predicting P-gp substrates are scarce. In this study, a support vector machine (SVM) method was developed to predict P-gp substrates and P-gp-substrate interactions, based on a training data set of 197 known P-gp substrates and non-substrates collected from the literature. We showed that the SVM method had a prediction accuracy of approximately 80% on an independent external validation data set of 32 compounds. A homology model of human P-gp based on the X-ray structure of mouse P-gp as a template has been constructed. We showed that molecular docking to the P-gp structures successfully predicted the geometry of P-gp-ligand complexes. Our SVM prediction and the molecular docking methods have been integrated into a free web server (http://pgp.althotas.com), which allows the users to predict whether a given compound is a P-gp substrate and how it binds to and interacts with P-gp. Utilization of such a web server may prove valuable for both rational drug design and screening.

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.

Central neurotoxicity in cancer chemotherapy: pharmacogenetic insights

Central neurotoxicity of chemotherapy is likely to be multifactorial. There are two hypotheses regarding endogenous mechanisms that may be involved, namely the target and the blood-brain barrier transporter hypotheses. Here, we will review candidate genetic determinants for the risk of chemotherapy-induced neurotoxicity, such as polymorphisms involved in the target mechanism. These include polymorphisms in folate metabolizing enzymes and apolipoprotein E, as well as those in blood-brain barrier transporter genes. Currently, the exact role of pharmacogenetics in mechanisms that lead to central neurotoxicity of chemotherapy has not been fully unraveled. Larger, prospective, longitudinal and more uniform studies are needed, with prechemotherapy and follow-up measurements of neuropsychological performance, MRI, PET, genetic profiles and biomarkers relevant for the proposed target and transporter mechanisms.

Sitamaquine overcomes ABC-mediated resistance to miltefosine and antimony in Leishmania

Although oral miltefosine represented an important therapeutic advance in the treatment of leishmaniasis, the appearance of resistance remains a serious threat. LMDR1/LABCB4, a P-glycoprotein-like transporter included in the Leishmania ABC (ATP-binding cassette) family, was the first molecule shown to be involved in experimental miltefosine resistance. LMDR1 pumps drugs out of the parasite, thereby decreasing their intracellular accumulation. Sitamaquine, another promising oral drug for leishmaniasis, is currently in phase 2b clinical trials. The physicochemical features of this drug suggested to us that it could be considered for use as an LMDR1 inhibitor. Indeed, we report herein that nonleishmanicidal concentrations of sitamaquine reverse miltefosine resistance in a multidrug resistance Leishmania tropica line that overexpresses LMDR1. This reversal effect is due to modulation of the LMDR1-mediated efflux of miltefosine. In addition, sitamaquine is not a substrate of LMDR1, as this transporter does not affect sitamaquine accumulation or sensitivity in the parasite. Likewise, we show that ketoconazole, another oral leishmanicidal drug known to interact with ABC transporters, is also able to reverse LMDR1-mediated miltefosine resistance, although with a lower efficiency than sitamaquine. Molecular docking on a three-dimensional homology model of LMDR1 showed different preferential binding sites for each substrate-inhibitor pair, thus explaining this different behavior. Finally, we show that sitamaquine is also able to modulate the antimony resistance mediated by MRPA/LABCC3, another ABC transporter involved in experimental and clinical antimony resistance in this parasite. Taken together, these data suggest that the combination of sitamaquine with miltefosine or antimony could avoid the appearance of resistance mediated by these membrane transporters in Leishmania.

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.

Conformational changes induced by ATP-hydrolysis in an ABC transporter: a molecular dynamics study of the Sav1866 exporter

ATP-Binding Cassette (ABC) transporters are ubiquitous membrane proteins that use energy from ATP binding or/and hydrolysis to actively transport allocrites across membranes. In this study, we identify ATP-hydrolysis induced conformational changes in a complete ABC exporter (Sav1866) from Staphylococcus aureaus, using molecular dynamics (MD) simulations. By performing MD simulations on the ATP and ADP+IP bound states, we identify the conformational consequences of hydrolysis, showing that the major rearrangements are not restricted to the NBDs, but extend to the transmembrane domains (TMDs) external regions. For the first time, to our knowledge, we see, within the context of a complete transporter, NBD dimer opening in the ADP+IP state in contrast with all ATP-bound states. This opening results from the dissociation of the ABC signature motif from the nucleotide. In addition, in both states, we observe the opening of a gate entrance in the intracellular loop region leading to the exposure of the TMDs internal cavity to the cytoplasm. To see if this opening was large enough to allow allocrite transport, the adiabatic energy profile for doxorubicin passage was determined. For both states, this profile, although an approximation, is overall downhill from the cytoplasmatic to the extracellular side, and the local energy barriers along the TMDs are relatively small, evidencing the exporter nature of Sav1866. The major difference between states is an energy barrier located in the cytoplasmic gate region, which becomes reduced upon hydrolysis, suggesting that allocrite passage is facilitated, and evidencing a possible molecular mechanism for the active transport in these proteins.

Homology modeling and binding site assessment of the human P-glycoprotein

Background: Due to its impact on multidrug resistance and pharmacokinetics P-glycoprotein (P-gp) has been identified as an important anti-target in pharmaceutical research. Recent publication of the mouse P-gp structure prompted us to build a new model for human P-gp and investigate its binding-site characteristics. Results: We developed and validated the human P-gp model that was used for induced-fit docking of experimentally characterized P-gp substrates. Residues located in the binding pocket are in good correlation with the results of side-directed mutagenesis studies. However, enrichment studies aimed at discriminating inhibitors and substrates from decoys resulted in only limited enrichments. Conclusion: A mouse P-gp-based homology model might be useful when analyzing protein–ligand interactions of known human P-gp substrates if induced-fit effects are considered.