Multiclass Classifier for P-Glycoprotein Substrates, Inhibitors, and Non-Active Compounds

Natural products can be potential inhibitors of metalloproteinase II from Bacteroides fragilis to intervene colorectal cancer

Bacteroides fragilis, a gram negative and obligate anaerobe bacterium, is a member of normal gut microbiota and facilitates many essential roles being performed in human body in normal circumstances specifically in Gastrointestinal or GI tract. Sometimes, due to genetics, epigenetics, and environmental factors, Bacteroides fragilis and their protein(s) start interacting with intestinal epithelium thus damaging the lining leading to colorectal cancers (CRC). To identify these protein(s), we incorporated a novel subtractive proteomics approach in the study. Metalloproteinase II (MPII), a Bacteroides fragilis toxin (bft), was investigated for its virulence and unique pathways to demonstrate its specificity and uniqueness in pathogenicity followed by molecular docking against a set of small drug-like natural molecules to discover potential inhibitors against the toxin. All these identified inhibitor-like molecules were analyzed for their ADMET calculations and detailed physiochemical properties to predict their druggability, GI absorption, blood brain barrier and skin permeation, and others. Resultantly, a total of ten compounds with the least binding energies were obtained and were subjected to protein-compound interaction analysis. Interaction analysis revealed the most common ligand-interacting residues in MPII are His 345, Glu 346, His 339, Gly 310, Tyr 341, Pro 340, Asp 187, Phe 309, Lys 307, Ile 185, Thr 308, and Pro 184. Therefore, top three compounds complexed with MPII having best binding energies were selected in order to analyze their trajectories. RMSD, RMSF, Rg and MMPBSA analysis revealed that all compounds showed good binding and keeping the complex stable and compact throughout the simulation time in addition to all properties and qualities of being a potential inhibitor against MPII.

Consensus screening for a challenging target: the quest for P-glycoprotein inhibitors

ATP-binding cassette (ABC) transporters are a large family of proteins involved in membrane transport of a wide variety of substrates. Among them, ABCB1, also known as MDR-1 or P-glycoprotein (P-gp), is the most characterized. By exporting xenobiotics out of the cell, P-gp activity can affect the ADME properties of several drugs. Moreover, P-gp has been found to mediate multidrug resistance in cancer cells. Thus, the inhibition of P-gp activity may lead to increased absorption and/or intracellular accumulation of co-administered drugs, enhancing their effectiveness. Using the human-mouse chimeric cryoEM 3D structure of the P-gp in the inhibitor-bound intermediate form (PDBID: 6qee), approximately 200 000 commercially available natural compounds from the ZINC database were virtually screened. To build a model able to discriminate between substrate and inhibitors, two datasets of compounds with known activity, including P-gp inhibitors, substrates, and inactive molecules were also docked. The best docking pose of selected substrates and inhibitors were used to generate 3D common feature pharmacophoric models that were combined with the Autodock Vina binding energy values to prioritize compounds for visual inspection. With this consensus approach, 13 potential candidates were identified and then tested for their ability to inhibit P-gp, using zosuquidar, a third generation P-gp inhibitor, as a reference drug. Eight compounds were found to be active with 6 of them having an IC50 lower than 5 μM in a membrane-based ATPase activity assay. Moreover, the P-gp inhibitory activity was also confirmed by two different cell-based in vitro methods. Both retrospective and prospective results demonstrate the ability of the combined structure-based pharmacophore modeling and docking-based virtual screening approach to predict novel hit compounds with inhibitory activity toward P-gp. The resulting chemical scaffolds could serve as inspiration for the optimization of novel and more potent P-gp inhibitors.

Plasma Concentrations of Rosmarinic Acid in Patients on Antiretroviral Therapy: In Silico Exploration Based on Clinical Data

Rosmarinic acid (RA) is a phenolic compound with antiviral properties, often encountered in dietary supplements and herbal drugs. Data on the pharmacokinetics of RA are lacking in cases of the chronic use of supplements containing this compound, and only limited data on the metabolism and distribution of RA are available. The aim of the study was to investigate the plasma levels of RA after 12 weeks of use and determine potential interactions of RA and selected antiretroviral drugs. Patients infected with human immunodeficiency virus took a supplement containing RA for 12 weeks, after which the RA concentrations in the plasma samples were analyzed. A detailed in silico analysis was conducted in order to elucidate the potential interactions between RA and the drugs efavirenz, darunavir and raltegravir. It was found that RA can be detected in patients' plasma samples, mainly in the form of sulphoglucuronide. The potential interactions are suggested on the level of liver metabolizing enzymes and efflux P-glycoprotein, with RA competing with antiretroviral drugs as a substrate in metabolism and distribution systems. The present study suggests that the simultaneous use of RA and antiretroviral therapy (containing efavirenz, darunavir or raltegravir) may affect the plasma levels of RA after prolonged supplementation.

Antiparasitic Activity of Plumbago auriculata Extracts and Its Naphthoquinone Plumbagin against Trypanosoma cruzi

Chagas disease (CD) caused by the protozoan Trypanosoma cruzi affects more than six million people worldwide. Treatment is restricted to benznidazole (Bz) and nifurtimox (Nf) that display low activity in the later chronic stage besides triggering toxic events that result in treatment abandonment. Therefore, new therapeutic options are necessary. In this scenario, natural products emerge as promising alternatives to treat CD. In the family Plumbaginaceae, Plumbago sp. exhibits a broad spectrum of biological and pharmacological activities. Thus, our main objective was to evaluate, in vitro and in silico, the biological effect of crude extracts of root and of aerial parts of P. auriculata, as well as its naphthoquinone Plumbagin (Pb) against T. cruzi. The phenotypic assays revealed potent activity of the root extract against different forms (trypomastigote and intracellular forms) and strains (Y and Tulahuen), with a compound concentration that reduced 50% of the number of the parasite (EC₅₀) values ranging from 1.9 to 3.9 µg/mL. In silico analysis showed that Pb is predicted to have good oral absorption and permeability in Caco2 cells, besides excellent probability of absorption by human intestinal cells, without toxic or mutagenic potential effects, not being predicted as a substrate or inhibitor of P-glycoprotein. Pb was as potent as Bz against intracellular forms and displayed a superior trypanosomicidal effect (about 10-fold) in bloodstream forms (EC₅₀ = 0.8 µM) as compared to the reference drug (8.5 µM). The cellular targets of Pb on T. cruzi were evaluated using electron microscopy assays and the findings on bloodstream trypomastigotes showed several cellular insults related to the autophagic process. Regarding toxicity in mammalian cells, the root extracts and the naphthoquinone present a moderate toxic profile on fibroblasts and cardiac cell lines. Then, aiming to reduce host toxicity, the root extract and Pb were tested in combination with Bz, and the data showed additive profiles with the sum of the fractional inhibitory concentration indexes (ΣFICIs) being 1.45 and 0.87, respectively. Thus, our work reveals the promising antiparasitic activity of Plumbago auriculata crude extracts and its purified naphthoquinone Plumbagin against different forms and strains of Trypanosoma cruzi in vitro.

Assessment of Malondialdehyde Belonging to Modulators and Substrates of the P-Glycoprotein Transporter Protein

In the study on cells of the Caco-2 line, the affiliation of malondialdehyde (MDA) to modulators and substrates of P-glycoprotein (Pgp) was assessed, and the biological role of Pgp in conditions of oxidative stress (OS) was studied. MDA was used at concentrations of 10, 50, 100, and 150 μM; OS was simulated by incubation with hydrogen peroxide (H₂O₂) at concentrations of 0.1-100 μM for 24 h. The relative amount of Pgp was evaluated by the Western blot hybridization, and the activity was estimated by the transport of its substrate fexofenadine (HPLC with UV detection, HPLC MS/MS). In this study, it was shown that MDA at concentrations of 10 and 50 μM and exposure duration of 24 h increases the relative amount and activity of Pgp by acting through CAR and PXR, and MDA can be transported by Pgp. The induction of Pgp under the action of MDA during the development of OS can have a protective significance, ensuring the removal of the peroxidation product from cells into the extracellular space and thereby increasing the viability of cells.

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.

Versatile in silico modelling of microplastics adsorption capacity in aqueous environment based on molecular descriptor and machine learning

To comprehensively evaluate the hazards of microplastics and their coexisting organic pollutants, the sorption capacity of microplastics is a major issue that is quantified through the microplastic-aqueous sorption coefficient (K_d). Almost all quantitative structure-property relationship (QSPR) models that describe K_d apply only to narrow, relatively homogeneous groups of reactants. Herein, non-hybrid QSPR-based models were developed to predict PE-water (K_PE-w), PE-seawater (K_PE-sw), PVC-water (K_PVC-w) and PP-seawater (K_PP-sw) sorption coefficients at different temperatures, with eight machine learning algorithms. Moreover, novel hybrid intelligent models for predicting K_d more accurately were innovatively developed by applying GA, PSO and AdaBoost algorithms to optimize MLP and ELM models. The results indicated that all three optimization algorithms could improve the robustness and predictability of the standalone MLP and ELM models. In all models trained with K_PE-w, K_PE-sw, K_PVC-w and K_PP-sw data sets, GBDT-1 and XGBoost-1 models, MLP-GA-2 and MLP-PSO-2 models, MLR-3 and MLR-4 models performed better in terms of goodness of fit (R_adj²: 0.907-0.999), robustness (Q_BOOT²: 0.900-0.937) and predictability (R_ext²: 0.889-0.970), respectively. Analyzing the descriptors revealed that temperature, lipophilicity, ionization potential and molecular size were correlated closely with the adsorption capacity of microplastics to organic pollutants. The proposed QSPR models may assist in initial environmental exposure assessments without imposing heavy costs in the early experimental phase.

Probabilistic modelling of developmental neurotoxicity based on a simplified adverse outcome pathway network

•

A developmental neurotoxicity Adverse Outcome Pathway network was simplified.

•

Common key events were chosen based on topology analysis and expert judgement.

•

Quantification of causal relationships was informed by key event relationships.

•

Various types of information were integrated for probability prediction.

•

Bayesian hierarchical modelling was applied for hazard identification.

In a century where toxicology and chemical risk assessment are embracing alternative methods to animal testing, there is an opportunity to understand the causal factors of neurodevelopmental disorders such as learning and memory disabilities in children, as a foundation to predict adverse effects. New testing paradigms, along with the advances in probabilistic modelling, can help with the formulation of mechanistically-driven hypotheses on how exposure to environmental chemicals could potentially lead to developmental neurotoxicity (DNT). This investigation aimed to develop a Bayesian hierarchical model of a simplified AOP network for DNT. The model predicted the probability that a compound induces each of three selected common key events (CKEs) of the simplified AOP network and the adverse outcome (AO) of DNT, taking into account correlations and causal relations informed by the key event relationships (KERs). A dataset of 88 compounds representing pharmaceuticals, industrial chemicals and pesticides was compiled including physicochemical properties as well as in silico and in vitro information. The Bayesian model was able to predict DNT potential with an accuracy of 76%, classifying the compounds into low, medium or high probability classes. The modelling workflow achieved three further goals: it dealt with missing values; accommodated unbalanced and correlated data; and followed the structure of a directed acyclic graph (DAG) to simulate the simplified AOP network. Overall, the model demonstrated the utility of Bayesian hierarchical modelling for the development of quantitative AOP (qAOP) models and for informing the use of new approach methodologies (NAMs) in chemical risk assessment.

Approaches to minimize the effects of P-glycoprotein in drug transport: A review

P-glycoprotein (P-gp) is a transporter protein that is come under the ATP binding cassette family of proteins. It is situated on the surface of the intestine epithelium, where P-gp substrate binds to the transporter and is pumped into the intestine lumen by the ATP-driven energy-dependent process. In this review, we summarize the role of the P-gp efflux transporter situated on the intestine, the clinical importance of P-gp related drug interactions, and approaches to minimize the effect of P-gp in drug transport. This review also focuses on the impact of P-gp on the bioavailability of the orally administered drug. Many drug's oral bioavailabilities can improve by concomitant use of P-gp inhibitors. Multidrug resistance are reduced by using some naturally occurring compounds obtained from plants and several synthetic P-gp inhibitors. Formulation strategies, one of the most important approaches to mimic the P-gp transporter's action, finally enhancing the oral bioavailability of the drug by inhibiting its P-gp efflux. Vitamin E TPGS, Gelucire 44/14 and other pharmaceutical/formulation excipients inhibit the P-gp efflux. A prodrug approach might be a useful strategy to overcome drug resistance. Prodrug helps to enhance the solubility or alter the pharmacokinetic properties but does not diminish the pharmacological action.

Structure-Function Relationships in the Human P-Glycoprotein (ABCB1): Insights from Molecular Dynamics Simulations

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.

ABCB1 Does Not Require the Side-Chain Hydrogen-Bond Donors Gln347, Gln725, Gln990 to Confer Cellular Resistance to the Anticancer Drug Taxol

The multidrug efflux transporter ABCB1 is clinically important for drug absorption and distribution and can be a determinant of chemotherapy failure. Recent structure data shows that three glutamines donate hydrogen bonds to coordinate taxol in the drug binding pocket. This is consistent with earlier drug structure-activity relationships that implicated the importance of hydrogen bonds in drug recognition by ABCB1. By replacing the glutamines with alanines we have tested whether any, or all, of Gln³⁴⁷, Gln⁷²⁵, and Gln⁹⁹⁰ are important for the transport of three different drug classes. Flow cytometric transport assays show that Q347A and Q990A act synergistically to reduce transport of Calcein-AM, BODIPY-verapamil, and OREGON GREEN-taxol bisacetate but the magnitude of the effect was dependent on the test drug and no combination of mutations completely abrogated function. Surprisingly, Q725A mutants generally improved transport of Calcein-AM and BODIPY-verapamil, suggesting that engagement of the wild-type Gln⁷²⁵ in a hydrogen bond is inhibitory for the transport mechanism. To test transport of unmodified taxol, stable expression of Q347/725A and the triple mutant was engineered and shown to confer equivalent resistance to the drug as the wild-type transporter, further indicating that none of these potential hydrogen bonds between transporter and transport substrate are critical for the function of ABCB1. The implications of the data for plasticity of the drug binding pocket are discussed.

Synthesis and Characterization of Some New Quinoxalin-2(1H)one and 2-Methyl-3H-quinazolin-4-one Derivatives Targeting the Onset and Progression of CRC with SRA, Molecular Docking, and ADMET Analyses

The pathogenesis of colorectal cancer is a multifactorial process. Dysbiosis and the overexpression of COX-2 and LDHA are important effectors in the initiation and development of the disease through chromosomal instability, PGE2 biosynthesis, and induction of the Warburg effect, respectively. Herein, we report the in vitro testing of some new quinoxalinone and quinazolinone Schiff’s bases as: antibacterial, COX-2 and LDHA inhibitors, and anticolorectal agents on HCT-116 and LoVo cells. Moreover, molecular docking and SAR analyses were performed to identify the structural features contributing to the biological activities. Among the synthesized molecules, the most active cytotoxic agent, (6d) was also a COX-2 inhibitor. In silico ADMET studies predicted that (6d) would have high Caco-2 permeability, and %HIA (99.58%), with low BBB permeability, zero hepatotoxicity, and zero risk of sudden cardiac arrest, or mutagenicity. Further, (6d) is not a potential P-gp substrate, instead, it is a possible P-gpI and II inhibitor, therefore, it can prevent or reverse the multidrug resistance of the anticancer drugs. Collectively, (6d) can be considered as a promising lead suitable for further optimization to develop anti-CRC agents or glycoproteins inhibitors.

Homology Modeling of the Human P-glycoprotein (ABCB1) and Insights into Ligand Binding through Molecular Docking Studies

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.