The ABCB1 transporter also known as P-glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette super-family of transporters; it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping the compounds out of cells. P-gp contributes to a decrease of toxicity and possesses broad substrate specificity. It is involved in the failure of numerous anticancer and antiviral chemotherapies due to the multidrug resistance (MDR) phenomenon, where it removes the chemotherapeutics out of the targeted cells. Understanding the details of the ligand–P-gp interaction is therefore crucial for the development of drugs that might overcome the MRD phenomenon and for obtaining a more effective prediction of the toxicity of certain compounds. In this work, an in silico modeling was performed using homology modeling and molecular docking methods with the aim of better understanding the ligand–P-gp interactions. Based on different mouse P-gp structural templates from the PDB repository, a 3D model of the human P-gp (hP-gp) was constructed by means of protein homology modeling. The homology model was then used to perform molecular docking calculations on a set of thirteen compounds, including some well-known compounds that interact with P-gp as substrates, inhibitors, or both. The sum of ranking differences (SRD) was employed for the comparison of the different scoring functions used in the docking calculations. A consensus-ranking scheme was employed for the selection of the top-ranked pose for each docked ligand. The docking results showed that a high number of π interactions, mainly π–sigma, π–alkyl, and π–π type of interactions, together with the simultaneous presence of hydrogen bond interactions contribute to the stability of the ligand–protein complex in the binding site. It was also observed that some interacting residues in hP-gp are the same when compared to those observed in a co-crystallized ligand (PBDE-100) with mouse P-gp (PDB ID: 4XWK). Our in silico approach is consistent with available experimental results regarding P-gp efflux transport assay; therefore it could be useful in the prediction of the role of new compounds in systemic toxicity.

P-glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette superfamily of transporters, and it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping compounds out of cells. P-gp contributes to a reduction in toxicity, and has broad substrate specificity. It is involved in the failure of many cancer and antiviral chemotherapies due to the phenomenon of multidrug resistance (MDR), in which the membrane transporter removes chemotherapeutic drugs from target cells. Understanding the details of the ligand–P-gp interaction is therefore critical for the development of drugs that can overcome the MDR phenomenon, for the early identification of P-gp substrates that will help us to obtain a more effective prediction of toxicity, and for the subsequent outdesign of substrate properties if needed. In this work, a series of molecular dynamics (MD) simulations of human P-gp (hP-gp) in an explicit membrane-and-water environment were performed to investigate the effects of binding different compounds on the conformational dynamics of P-gp. The results revealed significant differences in the behaviour of P-gp in the presence of active and non-active compounds within the binding pocket, as different patterns of movement were identified that could be correlated with conformational changes leading to the activation of the translocation mechanism. The predicted ligand–P-gp interactions are in good agreement with the available experimental data, as well as the estimation of the binding-free energies of the studied complexes, demonstrating the validity of the results derived from the MD simulations.

P-glycoprotein (P-gp) is a transmembrane protein that actively transports a wide variety of chemically diverse compounds out of the cell. It is highly associated with the ADMET (absorption, distribution, metabolism, excretion and toxicity) properties of drugs/drug candidates and contributes to decreasing toxicity by eliminating compounds from cells, thereby preventing intracellular accumulation. Therefore, in the drug discovery and toxicological assessment process it is advisable to pay attention to whether a compound under development could be transported by P-gp or not. In this study, an in silico multiclass classification model capable of predicting the probability of a compound to interact with P-gp was developed using a counter-propagation artificial neural network (CP ANN) based on a set of 2D molecular descriptors, as well as an extensive dataset of 2512 compounds (1178 P-gp inhibitors, 477 P-gp substrates and 857 P-gp non-active compounds). The model provided a good classification performance, producing non error rate (NER) values of 0.93 for the training set and 0.85 for the test set, while the average precision (AvPr) was 0.93 for the training set and 0.87 for the test set. An external validation set of 385 compounds was used to challenge the model’s performance. On the external validation set the NER and AvPr values were 0.70 for both indices. We believe that this in silico classifier could be effectively used as a reliable virtual screening tool for identifying potential P-gp ligands.