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

The Alternating Access Mechanism in Mammalian Multidrug Resistance Transporters and Their Bacterial Homologs

Multidrug resistance (MDR) proteins belonging to the ATP-Binding Cassette (ABC) transporter group play a crucial role in the export of cytotoxic drugs across cell membranes. These proteins are particularly fascinating due to their ability to confer drug resistance, which subsequently leads to the failure of therapeutic interventions and hinders successful treatments. One key mechanism by which multidrug resistance (MDR) proteins carry out their transport function is through alternating access. This mechanism involves intricate conformational changes that enable the binding and transport of substrates across cellular membranes. In this extensive review, we provide an overview of ABC transporters, including their classifications and structural similarities. We focus specifically on well-known mammalian multidrug resistance proteins such as MRP1 and Pgp (MDR1), as well as bacterial counterparts such as Sav1866 and lipid flippase MsbA. By exploring the structural and functional features of these MDR proteins, we shed light on the roles of their nucleotide-binding domains (NBDs) and transmembrane domains (TMDs) in the transport process. Notably, while the structures of NBDs in prokaryotic ABC proteins, such as Sav1866, MsbA, and mammalian Pgp, are identical, MRP1 exhibits distinct characteristics in its NBDs. Our review also emphasizes the importance of two ATP molecules for the formation of an interface between the two binding sites of NBD domains across all these transporters. ATP hydrolysis occurs following substrate transport and is vital for recycling the transporters in subsequent cycles of substrate transportation. Specifically, among the studied transporters, only NBD2 in MRP1 possesses the ability to hydrolyze ATP, while both NBDs of Pgp, Sav1866, and MsbA are capable of carrying out this reaction. Furthermore, we highlight recent advancements in the study of MDR proteins and the alternating access mechanism. We discuss the experimental and computational approaches utilized to investigate the structure and dynamics of MDR proteins, providing valuable insights into their conformational changes and substrate transport. This review not only contributes to an enhanced understanding of multidrug resistance proteins but also holds immense potential for guiding future research and facilitating the development of effective strategies to overcome multidrug resistance, thus improving therapeutic interventions.

Long Range Communication between the Drug-Binding Sites and Nucleotide Binding Domains of the Efflux Transporter ABCB1

The ABC efflux pump P-glycoprotein (P-gp) transports a wide variety of drugs and is inhibited by others. Some drugs stimulate ATP hydrolysis at the nucleotide binding domains (NBDs) and are transported, others uncouple ATP hydrolysis and transport, and others inhibit ATP hydrolysis. The molecular basis for the different behavior of these drugs is not well understood despite the availability of several structural models of P-gp complexes with ligands bound. Hypothetically, ligands differentially alter the conformational dynamics of peptide segments that mediate the coupling between the drug binding sites and the NBDs. Here, we explore by hydrogen-deuterium exchange mass spectrometry the dynamic consequences of a classic substrate and inhibitor, vinblastine and zosuquidar, binding to mouse P-gp (mdr1a) in lipid nanodiscs. The dynamics of P-gp in nucleotide-free, pre-hydrolysis, and post-hydrolysis states in the presence of each drug reveal distinct mechanisms of ATPase stimulation and implications for transport. For both drugs, there are common regions affected in a similar manner, suggesting that particular networks are the key to stimulating ATP hydrolysis. However, drug binding effects diverge in the post-hydrolysis state, particularly in the intracellular helices (ICHs 3 and 4) and neighboring transmembrane helices. The local dynamics and conformational equilibria in this region are critical for the coupling of drug binding and ATP hydrolysis and are differentially modulated in the catalytic cycle.

Study of Allosteric Transitions of Human P-Glycoprotein by Using the Two-State Anisotropic Network Model

Human P-glycoprotein (P-gp) is a kind of ATP-binding cassette (ABC) transporters. Once human P-gp is overexpressed in tumor cells, which can lead to tumor multidrug resistance (MDR). However, the present experimental methods are difficult to obtain the large-scale conformational transition process of human P-gp. In this work, we explored the allosteric pathway of human P-gp from the inward-facing (IF) to the outward-facing (OF) state in the substrate transport process with the two-state anisotropic network model (tANM). These results suggest that the allosteric transitions proceed in a coupled way. The conformational changes of nucleotide-binding domains (NBDs) finally make the transmembrane domains (TMDs) to the OF state via the role of the allosteric propagation of the intracellular helices IH1 and IH2. Additionally, this allosteric pathway is advantageous in energy compared with other methods. This study reveals the conformational transition of P-gp, which contributes to an understanding of the allosteric mechanism of ABC exporters.

Effect of the Force Field on Molecular Dynamics Simulations of the Multidrug Efflux Protein P-Glycoprotein

Molecular dynamics (MD) simulations have been used extensively to study P-glycoprotein (P-gp), a flexible multidrug transporter that is a key player in the development of multidrug resistance to chemotherapeutics. A substantial body of literature has grown from simulation studies that have employed various simulation conditions and parameters, including AMBER, CHARMM, OPLS, GROMOS, and coarse-grained force fields, drawing conclusions from simulations spanning hundreds of nanoseconds. Each force field is typically parametrized and validated on different data and observables, usually of small molecules and peptides; there have been few comparisons of force field performance on large protein-membrane systems. Here we compare the conformational ensembles of P-gp embedded in a POPC/cholesterol bilayer generated over 500 ns of replicate simulation with five force fields from popular biomolecular families: AMBER 99SB-ILDN, CHARMM 36, OPLS-AA/L, GROMOS 54A7, and MARTINI. We find considerable differences among the ensembles with little conformational overlap, although they correspond to similar extents to structural data obtained from electron paramagnetic resonance and cross-linking studies. Moreover, each trajectory was still sampling new conformations at a high rate after 500 ns of simulation, suggesting the need for more sampling. This work highlights the need to consider known limitations of the force field used (e.g., biases toward certain secondary structures) and the simulation itself (e.g., whether sufficient sampling has been achieved) when interpreting accumulated results of simulation studies of P-gp and other transport proteins.

Drug Resistance in Osteosarcoma: Emerging Biomarkers, Therapeutic Targets and Treatment Strategies

Simple Summary

Despite the adoption of aggressive, multimodal treatment schedules, the cure rate of high-grade osteosarcoma (HGOS) has not significantly improved in the last 30 years. The most relevant problem preventing improvement in HGOS prognosis is drug resistance. Therefore, validated novel biomarkers that help to identify those patients who could benefit from innovative treatment options and the development of drugs enabling personalized therapeutic protocols are necessary. The aim of this review was to give an overview on the most relevant emerging drug resistance-related biomarkers, therapeutic targets and new agents or novel candidate treatment strategies, which have been highlighted and suggested for HGOS to improve the success rate of clinical trials.

Abstract

High-grade osteosarcoma (HGOS), the most common primary malignant tumor of bone, is a highly aggressive neoplasm with a cure rate of approximately 40–50% in unselected patient populations. The major clinical problems opposing the cure of HGOS are the presence of inherent or acquired drug resistance and the development of metastasis. Since the drugs used in first-line chemotherapy protocols for HGOS and clinical outcome have not significantly evolved in the past three decades, there is an urgent need for new therapeutic biomarkers and targeted treatment strategies, which may increase the currently available spectrum of cure modalities. Unresponsive or chemoresistant (refractory) HGOS patients usually encounter a dismal prognosis, mostly because therapeutic options and drugs effective for rescue treatments are scarce. Tailored treatments for different subgroups of HGOS patients stratified according to drug resistance-related biomarkers thus appear as an option that may improve this situation. This review explores drug resistance-related biomarkers, therapeutic targets and new candidate treatment strategies, which have emerged in HGOS. In addition to consolidated biomarkers, specific attention has been paid to the role of non-coding RNAs, tumor-derived extracellular vesicles, and cancer stem cells as contributors to drug resistance in HGOS, in order to highlight new candidate markers and therapeutic targets. The possible use of new non-conventional drugs to overcome the main mechanisms of drug resistance in HGOS are finally discussed.

Physiological Functions of Bacterial "Multidrug" Efflux Pumps

Bacterial multidrug efflux pumps have come to prominence in human and veterinary pathogenesis because they help bacteria protect themselves against the antimicrobials used to overcome their infections. However, it is increasingly realized that many, probably most, such pumps have physiological roles that are distinct from protection of bacteria against antimicrobials administered by humans. Here we undertake a broad survey of the proteins involved, allied to detailed examples of their evolution, energetics, structures, chemical recognition, and molecular mechanisms, together with the experimental strategies that enable rapid and economical progress in understanding their true physiological roles. Once these roles are established, the knowledge can be harnessed to design more effective drugs, improve existing microbial production of drugs for clinical practice and of feedstocks for commercial exploitation, and even develop more sustainable biological processes that avoid, for example, utilization of petroleum.

Structural dynamics of ABC transporters: molecular simulation studies

The biological activities of living organisms involve various inputs and outputs. The ATP-driven substances (biomolecules) responsible for these kinds of activities through membrane (i.e. uptake and efflux of substrates) include ATP-binding cassette (ABC) transporters, some of which play important roles in multidrug resistance. The basic architecture of ABC transporters comprises transmembrane domains (TMDs) and nucleotide-binding domains (NBDs). The functional dynamics (substrate transport) of ABC transporters are realized by concerted motions, such as NBD dimerization, mechanical transmission via coupling helices (CHs), and the translocation of substrates through TMDs, which are induced by the binding and/or hydrolysis of ATP molecules and substrates. In this mini-review, we briefly discuss recent progresses in the structural dynamics as revealed by molecular simulation studies at all-atom (AA), coarse-grained (CG), and quantum mechanics/molecular mechanics (QM/MM) levels.

Identification and characterisation of putative drug binding sites in human ATP-binding cassette B5 (ABCB5) transporter

The human ATP-binding cassette B5 (ABCB5) transporter, a member of the ABC transporter superfamily, is linked to chemoresistance in tumour cells by drug effluxion. However, little is known about its structure and drug-binding sites. In this study, we generated an atomistic model of the full-length human ABCB5 transporter with the highest quality using the X-ray crystal structure of mouse ABCB1 (Pgp1), a close homologue of ABCB5 and a well-studied member of the ABC family. Molecular dynamics simulations were used to validate the atomistic model of ABCB5 and characterise its structural properties in model cell membranes. Molecular docking simulations of known ABCB5 substrates such as taxanes, anthracyclines, camptothecin and etoposide were then used to identify at least three putative binding sites for chemotherapeutic drugs transported by ABCB5. The location of these three binding sites is predicted to overlap with the corresponding binding sites in Pgp1. These findings will serve as the basis for future in vitro studies to validate the nature of the identified substrate-binding sites in the full-length ABCB5 transporter.

The role of the degenerate nucleotide binding site in type I ABC exporters

ATP‐binding cassette (ABC) transporters are fascinating molecular machines that are capable of transporting a large variety of chemically diverse compounds. The energy required for translocation is derived from binding and hydrolysis of ATP. All ABC transporters share a basic architecture and are composed of two transmembrane domains and two nucleotide binding domains (NBDs). The latter harbor all conserved sequence motifs that hallmark the ABC transporter superfamily. The NBDs form the nucleotide binding sites (NBSs) in their interface. Transporters with two active NBSs are called canonical transporters, while ABC exporters from eukaryotic organisms, including humans, frequently have a degenerate NBS1 containing noncanonical residues that strongly impair ATP hydrolysis. Here, we summarize current knowledge on degenerate ABC transporters. By integrating structural information with biophysical and biochemical evidence of asymmetric function, we develop a model for the transport cycle of degenerate ABC transporters. We will elaborate on the unclear functional advantages of a degenerate NBS.

The transport pathway in the ABCG2 protein and its regulation revealed by molecular dynamics simulations

Atomic-level structural insight on the human ABCG2 membrane protein, a pharmacologically important transporter, has been recently revealed by several key papers. In spite of the wealth of structural data, the pathway of transmembrane movement for the large variety of structurally different ABCG2 substrates and the physiological lipid regulation of the transporter has not been elucidated. The complex molecular dynamics simulations presented here may provide a breakthrough in understanding the steps of the substrate transport process and its regulation by cholesterol. Our analysis revealed drug binding cavities other than the central binding site and delineated a putative dynamic transport pathway for substrates with variable structures. We found that membrane cholesterol accelerated drug transport by promoting the closure of cytoplasmic protein regions. Since ABCG2 is present in all major biological barriers and drug-metabolizing organs, influences the pharmacokinetics of numerous clinically applied drugs, and plays a key role in uric acid extrusion, this information may significantly promote a reliable prediction of clinically important substrate characteristics and drug-drug interactions.

The critical size of gold nanoparticles for overcoming P-gp mediated multidrug resistance

Multidrug resistance (MDR) remains a huge obstacle during cancer treatment. One of the most studied MDR mechanisms is P-glycoprotein (P-gp) mediated drug efflux. Based on the three-dimensional structural characteristics of P-gp, gold nanoparticles (AuNPs) with average sizes of 4.1 nm and 5.4 nm were designed for the construction of nanodrug delivery systems (NanoDDSs), with the anticancer molecules 2-(9-anthracenylmethylene)-hydrazinecarbothioamide (ANS) and 6-mercaptopurine (6-MP) modified on the AuNP surfaces through the thiol group. In vitro cytotoxicity results suggested that the larger sized AuNPs can effectively decrease the drug resistance index of MCF-7/ADR cells to ∼2. Verapamil and P-gp antibody competitive experiments, combined with the cellular uptake of AuNPs, indicated that larger NanoDDSs were more conducive to intracellular drug accumulation and thus had improved anticancer activities, due to a size mismatch between the nanoparticles and the active site of P-gp, and, therefore, reduced drug efflux was seen. Measurements of ATPase activity and intracellular ATP levels indicated that the larger nanoparticles do not bind well to P-gp, thus avoiding effective recognition by P-gp. This was further evidenced by the observation that 4.1 nm and 5.4 nm NanoDDS-treated MCF-7/ADR cells showed remarkable differences in energy-related metabolic pathways. Therefore, the critical size of AuNPs for overcoming MDR was identified to be between 4.1 nm and 5.4 nm. This provides a more accurate description of the composite dimension requirements for NanoDDSs that are designed to overcome MDR.

Optimizing Targeted Inhibitors of P-Glycoprotein Using Computational and Structure-Guided Approaches

Overexpression of ABC transporters like P-glycoprotein (P-gp) has been correlated with resistances in cancer chemotherapy. Intensive efforts to identify P-gp inhibitors for use in combination therapy have not led to clinically approved inhibitors to date. Here, we describe computational approaches combined with structure-based design to improve the characteristics of a P-gp inhibitor previously identified by us. This hit compound represents a novel class of P-gp inhibitors that specifically targets and inhibits P-gp ATP hydrolysis while not being transported by the pump. We describe here a new program for virtual chemical synthesis and computational assessment, ChemGen, to produce hit compound variants with improved binding characteristics. The chemical syntheses of several variants, efficacy in reversing multidrug resistance in cell culture, and biochemical assessment of the inhibition mechanism are described. The usefulness of the computational predictions of binding characteristics of the inhibitor variants is discussed and compared to more traditional structure-based approaches.

ATP-binding Cassette Exporters: Structure and Mechanism with a Focus on P-glycoprotein and MRP1

BACKGROUND

Proteins that belong to the ATP-binding cassette superfamily include transporters that mediate the efflux of substrates from cells. Among these exporters, P-glycoprotein and MRP1 are involved in cancer multidrug resistance, protection from endo and xenobiotics, determination of drug pharmacokinetics, and the pathophysiology of a variety of disorders.

OBJECTIVE

To review the information available on ATP-binding cassette exporters, with a focus on Pglycoprotein, MRP1 and related proteins. We describe tissue localization and function of these transporters in health and disease, and discuss the mechanisms of substrate transport. We also correlate recent structural information with the function of the exporters, and discuss details of their molecular mechanism with a focus on the nucleotide-binding domains.

METHODS

Evaluation of selected publications on the structure and function of ATP-binding cassette proteins.

CONCLUSIONS

Conformational changes on the nucleotide-binding domains side of the exporters switch the accessibility of the substrate-binding pocket between the inside and outside, which is coupled to substrate efflux. However, there is no agreement on the magnitude and nature of the changes at the nucleotide- binding domains side that drive the alternate-accessibility. Comparison of the structures of Pglycoprotein and MRP1 helps explain differences in substrate selectivity and the bases for polyspecificity. P-glycoprotein substrates are hydrophobic and/or weak bases, and polyspecificity is explained by a flexible hydrophobic multi-binding site that has a few acidic patches. MRP1 substrates are mostly organic acids, and its polyspecificity is due to a single bipartite binding site that is flexible and displays positive charge.

Endoplasmic reticulum-targeting doxorubicin: a new tool effective against doxorubicin-resistant osteosarcoma

Doxorubicin is one of the most effective drugs for the first-line treatment of high-grade osteosarcoma. Several studies have demonstrated that the major cause for doxorubicin resistance in osteosarcoma is the increased expression of the drug efflux transporter ABCB1/P-glycoprotein (Pgp). We recently identified a library of H2S-releasing doxorubicins (Sdox) that were more effective than doxorubicin against resistant osteosarcoma cells. Here we investigated the molecular mechanisms of the higher efficacy of Sdox in human osteosarcoma cells with increasing resistance to doxorubicin. Differently from doxorubicin, Sdox preferentially accumulated within the endoplasmic reticulum (ER), and its accumulation was only modestly reduced in Pgp-expressing osteosarcoma cells. The increase in doxorubicin resistance was paralleled by the progressive down-regulation of genes of ER-associated protein degradation/ER-quality control (ERAD/ERQC), two processes that remove misfolded proteins and protect cell from ER stress-triggered apoptosis. Sdox, that sulfhydrated ER-associated proteins and promoted their subsequent ubiquitination, up-regulated ERAD/ERQC genes. This up-regulation, however, was insufficient to protect cells, since Sdox activated ER stress-dependent apoptotic pathways, e.g., the C/EBP-β LIP/CHOP/PUMA/caspases 12-7-3 axis. Sdox also promoted the sulfhydration of Pgp that was subsequently ubiquitinated: this process further enhanced Sdox retention and toxicity in resistant cells. Our work suggests that Sdox overcomes doxorubicin resistance in osteosarcoma cells by at least two mechanisms: it induces the degradation of Pgp following its sulfhydration and produces a huge misfolding of ER-associated proteins, triggering ER-dependent apoptosis. Sdox may represent the prototype of innovative anthracyclines, effective against doxorubicin-resistant/Pgp-expressing osteosarcoma cells by perturbing the ER functions.

Different structures of berberine and five other protoberberine alkaloids that affect P-glycoprotein-mediated efflux capacity

Berberine, berberrubine, thalifendine, demethyleneberberine, jatrorrhizine, and columbamine are six natural protoberberine alkaloid (PA) compounds that display extensive pharmacological properties and share the same protoberberine molecular skeleton with only slight substitution differences. The oral delivery of most PAs is hindered by their poor bioavailability, which is largely caused by P-glycoprotein (P-gp)-mediated drug efflux. Meanwhile, P-gp undergoes large-scale conformational changes (from an inward-facing to an outward-facing state) when transporting substrates, and these changes might strongly affect the P-gp-binding specificity. To confirm whether these six compounds are substrates of P-gp, to investigate the differences in efflux capacity caused by their trivial structural differences and to reveal the key to increasing their binding affinity to P-gp, we conducted a series of in vivo, in vitro, and in silico assays. Here, we first confirmed that all six compounds were substrates of P-gp by comparing the drug concentrations in wild-type and P-gp-knockout mice in vivo. The efflux capacity (net efflux) ranked as berberrubine > berberine > columbamine ~ jatrorrhizine > thalifendine > demethyleneberberine based on in vitro transport studies in Caco-2 monolayers. Using molecular dynamics simulation and molecular docking techniques, we determined the transport pathways of the six compounds and their binding affinities to P-gp. The results suggested that at the early binding stage, different hydrophobic and electrostatic interactions collectively differentiate the binding affinities of the compounds to P-gp, whereas electrostatic interactions are the main determinant at the late release stage. In addition to hydrophobic interactions, hydrogen bonds play an important role in discriminating the binding affinities.

Allosteric Role of Substrate Occupancy Toward the Alignment of P-glycoprotein Nucleotide Binding Domains

P-glycoprotein (Pgp) is an ATP-binding cassette transporter that eliminates toxins from the cell but causes multidrug resistance in chemotherapies. The crystal structures of Pgp revealed drug-like compounds bound to an inward-facing conformation in which the energy-harnessing nucleotide binding domains (NBDs) were widely separated with no interfacial interaction. Following drug binding, inward-facing Pgp must transition to an NBD dimer conformation to achieve ATP binding and hydrolysis at canonical sites defined by both halves of the interface. However, given the high degree of flexibility shown for this transporter, it is difficult to envision how NBDs overcome entropic considerations for achieving proper alignment in order to form the canonical ATP binding site. We explored the hypothesis that substrate occupancy of the polyspecific drug-binding cavity plays a role in the proper alignment of NBDs using computational approaches. We conducted twelve atomistic molecular dynamics (MD) simulations (100-300 ns) on inward-facing Pgp in a lipid bilayer with and without small molecule substrates to ascertain effects of drug occupancy on NBD dimerization. Both apo- and drug-occupied simulations showed NBDs approaching each other compared to the crystal structures. Apo-Pgp reached a pseudo-dimerization in which NBD signature motifs for ATP binding exhibited a significant misalignment during closure. In contrast, occupancy of three established substrates positioned by molecular docking achieved NBD alignment that was much more compatible with a canonical NBD dimerization trajectory. Additionally, aromatic amino acids, known to confer the polyspecific drug-binding characteristic of the internal pocket, may also govern polyspecific drug access to the cavity. The enrichment of aromatics comprising the TM4-TM6 portal suggested a preferential pathway over the aromatic-poor TM10-TM12 for lateral drug entry from the lipid bilayer. Our study also suggested that drug polyspecificity is enhanced due to a synergism between multiple drug-domain interactions involving 36 residues identified in TM1, 5, 6, 7, 11 and 12.

Evidence for the critical role of transmembrane helices 1 and 7 in substrate transport by human P-glycoprotein (ABCB1)

P-glycoprotein (P-gp) is an ABC transporter that exports many amphipathic or hydrophobic compounds, including chemically and functionally dissimilar anticancer drugs, from cells. To understand the role of transmembrane helices (TMH) 1 and 7 in drug-binding and transport, we selected six residues from both TMH1 (V53, I59, I60, L65, M68 and F72) and TMH7 (V713, I719, I720, Q725, F728 and F732); and substituted them with alanine by gene synthesis to generate a variant termed "TMH1,7 mutant P-gp". The expression and function of TMH1,7 mutant P-gp with twelve mutations was characterized using the BacMam baculovirus-HeLa cell expression system. The expression and conformation of TMH1,7 mutant P-gp was not altered by the introduction of the twelve mutations, as confirmed by using the human P-gp-specific antibodies UIC2, MRK16 and 4E3. We tested 25 fluorescently-labeled substrates and found that only three substrates, NBD-cyclosporine A, Rhod-2-AM and X-Rhod-1-AM were transported by the TMH1,7 mutant. The basal ATPase activity of TMH1,7 mutant P-gp was lower (40-50%) compared to wild-type (WT) P-gp, despite similar level of expression. Although most of the substrates modulate ATPase activity of P-gp, the activity of TMH1,7 mutant transporter was not significantly modulated by any of the tested substrates. Docking of selected substrates in homology models showed comparable docking scores for the TMH1,7 mutant and WT P-gp, although the binding conformations were different. Both the ATPase assay and in silico docking analyses suggest that the interactions with residues in the drug-binding pocket are altered as a consequence of the mutations. We demonstrate that it is possible to generate a variant of P-gp with a loss of broad substrate specificity and propose that TMH1 and TMH7 play a critical role in the drug efflux function of this multidrug transporter.

Exploring movement and energy in human P-glycoprotein conformational rearrangement

Human P-glycoprotein (P-gp), a kind of ATP-Binding Cassette transporter, can export a diverse variety of anti-cancer drugs out of the tumor cell. Its overexpression is one of the main reasons for the multidrug resistance (MDR) of tumor cells. It has been confirmed that during the substrate transport process, P-gp experiences a large-scale structural rearrangement from the inward- to outward-facing states. However, the mechanism of how the nucleotide-binding domains (NBDs) control the transmembrane domains (TMDs) to open towards the periplasm in the outward-facing state has not yet been fully characterized. Herein, targeted molecular dynamics simulations were performed to explore the conformational rearrangement of human P-gp. The results show that the allosteric process proceeds in a coupled way, and first the transition is driven by the NBDs, and then transmitted to the cytoplasmic parts of TMDs, finally to the periplasmic parts. The trajectories show that besides the translational motions, the NBDs undergo a rotation movement, which mainly occurs in xy plane and ensures the formation of the correct ATP-binding pockets. The analyses on the interaction energies between the six structure segments (cICLs) from the TMDs and NBDs reveal that their subtle energy differences play an important role in causing the periplasmic parts of the transmembrane helices to separate from each other in the established directions and in appropriate amplitudes. This conclusion can explain the two experimental phenomena about human P-gp in some extent. These studies have provided a detailed exploration into human P-gp rearrangement process and given an energy insight into the TMD reorientation during P-gp transition.

Conformational dynamics of P-glycoprotein in lipid nanodiscs and detergent micelles reveal complex motions on a wide time scale

P-glycoprotein (P-gp) is a highly substrate-promiscuous efflux transporter that plays a critical role in drug disposition. P-gp utilizes ATP hydrolysis by nucleotide-binding domains (NBDs) to drive transitions between inward-facing (IF) conformations that bind drugs and outward-facing (OF) conformations that release them to the extracellular solution. However, the details of the protein dynamics within either macroscopic IF or OF conformation remain uncharacterized, and the functional role of local dynamics has not been determined. In this work we measured the local dynamics of the IF state of P-gp in lipid nanodiscs and in detergent solution by hydrogen-deuterium (H/D) exchange MS. We observed "EX1 exchange kinetics," or bimodal kinetics, for several peptides distributed in both NBDs, particularly for P-gp in the lipid nanodiscs. Remarkably, the EX1 kinetics occurred on several time scales, ranging from seconds to hours, suggesting highly complex, and correlated, motions. The results indicate at least three distinct conformational states in the ligand-free P-gp and suggest a rough conformational landscape. Addition of excess ATP and vanadate, to favor the OF conformations, caused a generalized, but modest, decrease in H/D exchange throughout the NBDs and slowed the EX1 kinetic transitions of several peptides. The functional implications of the results are consistent with the possibility that conformational selection provides a source of substrate promiscuity.

The reliability of molecular dynamics simulations of the multidrug transporter P-glycoprotein in a membrane environment

Despite decades of research, the mechanism of action of the ABC multidrug transporter P-glycoprotein (P-gp) remains elusive. Due to experimental limitations, many researchers have turned to molecular dynamics simulation studies in order to investigate different aspects of P-gp function. However, such studies are challenging and caution is required when interpreting the results. P-gp is highly flexible and the time scale on which it can be simulated is limited. There is also uncertainty regarding the accuracy of the various crystal structures available, let alone the structure of the protein in a physiologically relevant environment. In this study, three alternative structural models of mouse P-gp (3G5U, 4KSB, 4M1M), all resolved to 3.8 Å, were used to initiate sets of simulations of P-gp in a membrane environment in order to determine: a) the sensitivity of the results to differences in the starting configuration; and b) the extent to which converged results could be expected on the times scales commonly simulated for this system. The simulations suggest that the arrangement of the nucleotide binding domains (NBDs) observed in the crystal structures is not stable in a membrane environment. In all simulations, the NBDs rapidly associated (within 10 ns) and changes within the transmembrane helices were observed. The secondary structure within the transmembrane domain was best preserved in the 4M1M model under the simulation conditions used. However, the extent to which replicate simulations diverged on a 100 to 200 ns timescale meant that it was not possible to draw definitive conclusions as to which structure overall was most stable, or to obtain converged and reliable results for any of the properties examined. The work brings into question the reliability of conclusions made in regard to the nature of specific interactions inferred from previous simulation studies on this system involving similar sampling times. It also highlights the need to demonstrate the statistical significance of any results obtained in simulations of large flexible proteins, especially where the initial structure is uncertain.

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.

Mechanisms of TPGS and its derivatives inhibiting P-glycoprotein efflux pump and application for reversing multidrug resistance in hepatocellular carcinoma

We have developed a TPGS–GA conjugate and TPGS–LA conjugate which possess more effective P-gp inhibition compared to TPGS because of the enhancement of hydrophilicity and negative charge.

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.

Luminescence resonance energy transfer spectroscopy of ATP-binding cassette proteins

The ATP-binding cassette (ABC) superfamily includes regulatory and transport proteins. Most human ABC exporters pump substrates out of cells using energy from ATP hydrolysis. Although major advances have been made toward understanding the molecular mechanism of ABC exporters, there are still many issues unresolved. During the last few years, luminescence resonance energy transfer has been used to detect conformational changes in real time, with atomic resolution, in isolated ABC nucleotide binding domains (NBDs) and full-length ABC exporters. NBDs are particularly interesting because they provide the power stroke for substrate transport. Luminescence resonance energy transfer (LRET) is a spectroscopic technique that can provide dynamic information with atomic-resolution of protein conformational changes under physiological conditions. Using LRET, it has been shown that NBD dimerization, a critical step in ABC proteins catalytic cycle, requires binding of ATP to two nucleotide binding sites. However, hydrolysis at just one of the sites can drive dissociation of the NBD dimer. It was also found that the NBDs of the bacterial ABC exporter MsbA reconstituted in a lipid bilayer membrane and studied at 37°C never separate as much as suggested by crystal structures. This observation stresses the importance of performing structural/functional studies of ABC exporters under physiologic conditions. 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.

Strategies of Drug Transporter Quantitation by LC-MS: Importance of Peptide Selection and Digestion Efficiency

Huge variation of drug transporter abundance was seen in the literature, making PBPK prediction difficult when transporters play a major role. Among multiple factors such as membrane fraction, digestion, and peptide selection that contributed to such variation, peptide selection is the least discussed. Herein, a strategy was established by using a small amount of purified protein standard to select a peptide with near 100% digestion efficiency for quantitation of a transporter protein MDR1. The impact of native membrane protein's tertiary structure on the digestion efficiency of surrogate peptides of MDR1 was investigated. Peptides in more solvent accessible regions are found to be digested much more efficiently than those in large stretches of helical structures. The concentration of peptide EALDESIPPVSFWR(EAL) in the most solvent accessible linker region of MDR1 was found closest to the true protein concentration. When using EAL for MDR1 quantitation, the abundance is over 10 times higher than previously reported, indicating the importance of peptide selection for transporter quantitation. In addition, this study also proposes a screening strategy to select peptides appropriate for relative quantitation for in vitro-in vivo extrapolation in the absence of any protein standard.

Structural Dynamics of the Heterodimeric ABC Transporter TM287/288 Induced by ATP and Substrate Binding

TM287/288 is a heterodimeric ATP-binding cassette (ABC) transporter, which harnesses the energy of ATP binding and hydrolysis at the nucleotide-binding domains (NBDs) to transport a wide variety of molecules through the transmembrane domains (TMDs) by alternating inward- and outward-facing conformations. Here, we conducted multiple 100 ns molecular dynamics simulations of TM287/288 in different ATP- and substrate-bound states to elucidate the effects of ATP and substrate binding. As a result, the binding of two ATP molecules to the NBDs induced the formation of the consensus ATP-binding pocket (ABP2) or the NBD dimerization, whereas these processes did not occur in the presence of a single ATP molecule or when the protein was in its apo state. Moreover, binding of the substrate to the TMDs enhanced the formation of ABP2 through allosteric TMD-NBD communication. Furthermore, in the apo state, α-helical subdomains of the NBDs approached each other, acquiring a conformation with core half-pockets exposed to the solvent, appropriate for ATP binding. We propose a "core-exposed" model for this novel conformation found in the apo state of ABC transporters. These findings provide important insights into the structural dynamics of ABC transporters.

Homology modelling of human P-glycoprotein

P-glycoprotein (P-gp) is an ATP-binding cassette transporter that exports a huge range of compounds out of cells and is thus one of the key proteins in conferring multi-drug resistance in cancer. Understanding how it achieves such a broad specificity and the series of conformational changes that allow export to occur form major, on-going, research objectives around the world. Much of our knowledge to date has been derived from mutagenesis and assay data. However, in recent years, there has also been great progress in structural biology and although the structure of human P-gp has not yet been solved, there are now a handful of related structures on which homology models can be built to aid in the interpretation of the vast amount of experimental data that currently exists. Many models for P-gp have been built with this aim, but the situation is complicated by the apparent flexibility of the system and by the fact that although many potential templates exist, there is large variation in the conformational state in which they have been crystallized. In this review, we summarize how homology modelling has been used in the past, how models are typically selected and finally illustrate how MD simulations can be used as a means to give more confidence about models that have been generated via this approach.