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

Study Models of Drug-Drug Interactions Involving P-Glycoprotein: The Potential Benefit of P-Glycoprotein Modulation at the Kidney and Intestinal Levels

P-glycoprotein (P-gp) is a crucial membrane transporter situated on the cell's apical surface, being responsible for eliminating xenobiotics and endobiotics. P-gp modulators are compounds that can directly or indirectly affect this protein, leading to changes in its expression and function. These modulators can act as inhibitors, inducers, or activators, potentially causing drug-drug interactions (DDIs). This comprehensive review explores diverse models and techniques used to assess drug-induced P-gp modulation. We cover several approaches, including in silico, in vitro, ex vivo, and in vivo methods, with their respective strengths and limitations. Additionally, we explore the therapeutic implications of DDIs involving P-gp, with a special focus on the renal and intestinal elimination of P-gp substrates. This involves enhancing the removal of toxic substances from proximal tubular epithelial cells into the urine or increasing the transport of compounds from enterocytes into the intestinal lumen, thereby facilitating their excretion in the feces. A better understanding of these interactions, and of the distinct techniques applied for their study, will be of utmost importance for optimizing drug therapy, consequently minimizing drug-induced adverse and toxic effects.

Cross-Linking Mass Spectrometry on P-Glycoprotein

The ABC transporter P-glycoprotein (Pgp) has been found to be involved in multidrug resistance in tumor cells. Lipids and cholesterol have a pivotal role in Pgp's conformations; however, it is often difficult to investigate it with conventional structural biology techniques. Here, we applied robust approaches coupled with cross-linking mass spectrometry (XL-MS), where the natural lipid environment remains quasi-intact. Two experimental approaches were carried out using different cross-linkers (i) on living cells, followed by membrane preparation and immunoprecipitation enrichment of Pgp, and (ii) on-bead, subsequent to membrane preparation and immunoprecipitation. Pgp-containing complexes were enriched employing extracellular monoclonal anti-Pgp antibodies on magnetic beads, followed by on-bead enzymatic digestion. The LC-MS/MS results revealed mono-links on Pgp's solvent-accessible residues, while intraprotein cross-links confirmed a complex interplay between extracellular, transmembrane, and intracellular segments of the protein, of which several have been reported to be connected to cholesterol. Harnessing the MS results and those of molecular docking, we suggest an epitope for the 15D3 cholesterol-dependent mouse monoclonal antibody. Additionally, enriched neighbors of Pgp prove the strong connection of Pgp to the cytoskeleton and other cholesterol-regulated proteins. These findings suggest that XL-MS may be utilized for protein structure and network analyses in such convoluted systems as membrane proteins.

Probing the Allosteric Modulation of P-Glycoprotein: A Medicinal Chemistry Approach Toward the Identification of Noncompetitive P-Gp Inhibitors

A medicinal chemistry approach combining in silico and in vitro methodologies was performed aiming at identifying and characterizing putative allosteric drug-binding sites (aDBSs) at the interface of the transmembrane- and nucleotide-binding domains (TMD-NBD) of P-glycoprotein. Two aDBSs were identified, one in TMD1/NBD1 and another one in TMD2/NBD2, by means of in silico fragment-based molecular dynamics and characterized in terms of size, polarity, and lining residues. From a small library of thioxanthone and flavanone derivatives, experimentally described to bind at the TMD-NBD interfaces, several compounds were identified to be able to decrease the verapamil-stimulated ATPase activity. An IC₅₀ of 81 ± 6.6 μM is reported for a flavanone derivative in the ATPase assays, providing evidence for an allosteric efflux modulation in P-glycoprotein. Molecular docking and molecular dynamics gave additional insights on the binding mode on how flavanone derivatives may act as allosteric inhibitors.

Long-range communication between transmembrane- and nucleotide-binding domains does not depend on drug binding to mutant P-glycoprotein

In this study, the impact of four P-gp mutations (G185V, G830V, F978A and ΔF335) on drug-binding and efflux-related signal-transmission mechanism was comprehensively evaluated in the presence of ligands within the drug-binding pocket (DBP), experimentally related with changes in their drug efflux profiles. The severe repacking of the transmembrane helices (TMH), induced by mutations and exacerbated by the presence of ligands, indicates that P-gp is sensitive to perturbations in the transmembrane region. Alterations on drug-binding were also observed as a consequence of the TMH repacking, but were not always correlated with alterations on ligands binding mode and/or binding affinity. Finally, and although all P-gp variants holo systems showed considerable changes in the intracellular coupling helices/nucleotide-binding domain (ICH-NBD) interactions, they seem to be primarily induced by the mutation itself rather than by the presence of ligands within the DBP. The data further suggest that the changes in drug efflux experimentally reported are mostly related with changes on drug specificity rather than effects on signal-transmission mechanism. We also hypothesize that an increase in the drug-binding affinity may also be related with the decreased drug efflux, while minor changes in binding affinities are possibly related with the increased drug efflux observed in transfected cells.Communicated by Ramaswamy H. Sarma.

A Structure-Based View on ABC-Transporter Linked to Multidrug Resistance

The discovery of the first ATP-binding cassette (ABC) transporter, whose overexpression in cancer cells is responsible for exporting anticancer drugs out of tumor cells, initiated enormous efforts to overcome tumor cell multidrug resistance (MDR) by inhibition of ABC-transporter. Because of its many physiological functions, diverse studies have been conducted on the mechanism, function and regulation of this important group of transmembrane transport proteins. In this review, we will focus on the structural aspects of this transporter superfamily. Since the resolution revolution of electron microscope, experimentally solved structures increased rapidly. A summary of the structures available and an overview of recent structure-based studies are provided. More specifically, the artificial intelligence (AI)-based predictions from AlphaFold-2 will be discussed.

Interaction of a Homologous Series of Amphiphiles with P-glycoprotein in a Membrane Environment-Contributions of Polar and Non-Polar Interactions

The transport of drugs by efflux transporters in biomembranes limits their bioavailability and is a major determinant of drug resistance development by cancer cells and pathogens. A large number of chemically dissimilar drugs are transported, and despite extensive studies, the molecular determinants of substrate specificity are still not well understood. In this work, we explore the role of polar and non-polar interactions on the interaction of a homologous series of fluorescent amphiphiles with the efflux transporter P-glycoprotein. The interaction of the amphiphiles with P-glycoprotein is evaluated through effects on ATPase activity, efficiency in inhibition of [¹²⁵I]-IAAP binding, and partition to the whole native membranes containing the transporter. The results were complemented with partition to model membranes with a representative lipid composition, and details on the interactions established were obtained from MD simulations. We show that when the total concentration of amphiphile is considered, the binding parameters obtained are apparent and do not reflect the affinity for P-gp. A new formalism is proposed that includes sequestration of the amphiphiles in the lipid bilayer and the possible binding of several molecules in P-gp's substrate-binding pocket. The intrinsic binding affinity thus obtained is essentially independent of amphiphile hydrophobicity, highlighting the importance of polar interactions. An increase in the lipophilicity and amphiphilicity led to a more efficient association with the lipid bilayer, which maintains the non-polar groups of the amphiphiles in the bilayer, while the polar groups interact with P-gp's binding pocket. The presence of several amphiphiles in this orientation is proposed as a mechanism for inhibition of P-pg function.

Recent Advances on P-Glycoprotein (ABCB1) Transporter Modelling with In Silico Methods

ABC transporters play a critical role in both drug bioavailability and toxicity, and with the discovery of the P-glycoprotein (P-gp), this became even more evident, as it plays an important role in preventing intracellular accumulation of toxic compounds. Over the past 30 years, intensive studies have been conducted to find new therapeutic molecules to reverse the phenomenon of multidrug resistance (MDR) ), that research has found is often associated with overexpression of P-gp, the most extensively studied drug efflux transporter; in MDR, therapeutic drugs are prevented from reaching their targets due to active efflux from the cell. The development of P-gp inhibitors is recognized as a good way to reverse this type of MDR, which has been the subject of extensive studies over the past few decades. Despite the progress made, no effective P-gp inhibitors to reverse multidrug resistance are yet on the market, mainly because of their toxic effects. Computational studies can accelerate this process, and in silico models such as QSAR models that predict the activity of compounds associated with P-gp (or analogous transporters) are of great value in the early stages of drug development, along with molecular modelling methods, which provide a way to explain how these molecules interact with the ABC transporter. This review highlights recent advances in computational P-gp research, spanning the last five years to 2022. Particular attention is given to the use of machine-learning approaches, drug-transporter interactions, and recent discoveries of potential P-gp inhibitors that could act as modulators of multidrug resistance.

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.

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.

Structure activity relationships and the binding mode of quinolinone-pyrimidine hybrids as reversal agents of multidrug resistance mediated by P-gp

P-gp-associated multidrug resistance is a major impediment to the success of chemotherapy. With the aim of finding non-toxic and effective P-gp inhibitors, we investigated a panel of quinolin-2-one-pyrimidine hybrids. Among the active compounds, two of them significantly increased intracellular doxorubicin and rhodamine 123 accumulation by inhibiting the efflux mediated by P-gp and restored doxorubicin toxicity at nanomolar range. Structure-activity relationships showed that the number of methoxy groups, an optimal length of the molecule in its extended conformation, and at least one flexible methylene group bridging the quinolinone to the moiety bearing the pyrimidine favored the inhibitory potency of P-gp. The best compounds showed a similar binding pattern and interactions to those of doxorubicin and tariquidar, as revealed by MD and hybrid QM/MM simulations performed with the recent experimental structure of P-gp co-crystallized with paclitaxel. Analysis of the molecular interactions stabilizing the different molecular complexes determined by MD and QTAIM showed that binding to key residues from TMH 4-7 and 12 is required for inhibition.

Screening of Natural Compounds as P-Glycoprotein Inhibitors against Multidrug Resistance

Multidrug resistance (MDR) is a common problem when fighting cancer with chemotherapy. P-glycoprotein (P-gp, or MDR1) is an active pump responsible for the efflux of xenobiotics out of the cell, including anti-cancer drugs. It is a validated target against MDR. No crystal structure of the human P-gp is available to date, and only recently several cryo-EM structures have been solved. In this paper, we present a comprehensive computational approach that includes constructing the full-length three-dimensional structure of the human P-gp and its refinement using molecular dynamics. We assessed its flexibility and conformational diversity, compiling a dynamical ensemble that was used to dock a set of lignan compounds, previously reported as active P-gp inhibitors, and disclose their binding modes. Based on the statistical analysis of the docking results, we selected a system for performing the structure-based virtual screening of new potential P-gp inhibitors. We tested the method on a library of 87 natural flavonoids described in the literature, and 10 of those were experimentally assayed. The results reproduced the theoretical predictions only partially due to various possible factors. However, at least two of the predicted natural flavonoids were demonstrated to be effective P-gp inhibitors. They were able to increase the accumulation of doxorubicin inside the human promyelocytic leukemia HL60/MDR cells overexpressing P-gp and potentiate the antiproliferative activity of this anti-cancer drug.

Lysine 268 adjacent to transmembrane helix 5 of hamster P-glycoprotein is the major photobinding site of iodomycin in CHO B30 cells

P‐glycoprotein (Pgp) detoxifies cells by exporting hundreds of chemically dissimilar hydrophobic and amphipathic compounds and is implicated in multidrug resistance (MDR) in the treatment of cancers. Photoaffinity labeling of plasma membrane vesicles of MDR CHO B30 cells with the anthracycline [¹²⁵I]‐iodomycin, subsequent sequential cleavage with BNPS‐skatol and endoproteinase Lys‐C, and the Edman sequencing of the purified photoaffinity‐labeled peptide identified the lysine residue at position 268 in the hamster Pgp primary sequence as the major photobinding site of iodomycin in CHO B30 cells. Lysine 268 is located adjacent to the cytosolic terminus of transmembrane 5. According to thermodynamic and kinetic analyses, this location should present the equilibrium binding site of ATP‐free Pgp for daunomycin and iodomycin in B30 cells.

Cholesterol Asymmetrically Modulates the Conformational Ensemble of the Nucleotide-Binding Domains of P-Glycoprotein in Lipid Nanodiscs

P-Glycoprotein (P-gp) is an ATP-dependent efflux pump that clears a wide variety of drugs and toxins from cells. P-gp undergoes large-scale structural changes and demonstrates conformational heterogeneity even within a single catalytic or drug-bound state, although the role of heterogeneity remains unclear. P-gp is found in a variety of cell types that vary in lipid composition, which modulates its activity. An understanding of structural or dynamic changes due to the lipid environment is lacking. We aimed to determine the effects of cholesterol in a membrane on the conformational behavior of P-gp in lipid nanodiscs. The presence of cholesterol stimulates ATP hydrolysis and alters lipid order and fluidity. Hydrogen/deuterium exchange mass spectrometry demonstrates that cholesterol in the membrane induces asymmetric, long-range changes in the distributions and exchange kinetics of conformations of the nucleotide-binding domains, correlating the effects of lipid composition on activity with specific changes in the P-gp conformational landscape.