In Silico Targeting Human Multidrug Transporter ABCG2 in Breast Cancer: Database Screening, Molecular Docking, and Molecular Dynamics Study

Foretinib, a c-MET receptor tyrosine kinase inhibitor, tackles multidrug resistance in cancer cells by inhibiting ABCB1 and ABCG2 transporters

BACKGROUND

ABC transporter-mediated multidrug resistance (MDR) remains a major obstacle for cancer pharmacological treatment. Some tyrosine kinase inhibitors (TKIs) have been shown to reverse MDR. The present study was designed to evaluate for the first time whether foretinib, a multitargeted TKI, can circumvent ABCB1 and ABCG2-mediated MDR in treatment-resistant cancer models.

METHODS

Accumulation of fluorescent substrates of ABCB1 and ABCG2 in ABCB1-overexpressing MES-SA/DX5 and ABCG2-overexpressing MCF-7/MX and their parenteral cells was evaluated by flow cytometry. The growth inhibitory activity of single and combination therapy of foretinib and chemotherapeutic drugs on MDR cells was examined by MTT assay. Analysis of combined interaction effects was performed using CalcuSyn software.

RESULTS

It was firstly proved that foretinib increased the intracellular accumulation of rhodamine 123 and mitoxantrone in MES-SA/DX5 and MCF-7/MX cancer cells, with accumulation ratios of 12 and 2.2 at 25 μM concentration, respectively. However, it did not affect the accumulation of fluorescent substrates in the parental cells. Moreover, foretinib synergistically improved the cytotoxic effects of doxorubicin and mitoxantrone. The means of combination index (CI) values at fraction affected (Fa) values of 0.5, 0.75, and 0.9 were 0.64 ± 0.08 and 0.47 ± 0.09, in MES-SA/DX5 and MCF-7/MX cancer cells, respectively. In silico analysis also suggested that the drug-binding domain of ABCB1 and ABCG2 transporters could be considered as potential target for foretinib.

CONCLUSION

Overall, our results suggest that foretinib can target MDR-linked ABCB1 and ABCG2 transporters in clinical cancer therapy.

Computational and artificial intelligence-based approaches for drug metabolism and transport prediction

Drug metabolism and transport, orchestrated by drug-metabolizing enzymes (DMEs) and drug transporters (DTs), are implicated in drug-drug interactions (DDIs) and adverse drug reactions (ADRs). Reliable and precise predictions of DDIs and ADRs are critical in the early stages of drug development to reduce the rate of drug candidate failure. A variety of experimental and computational technologies have been developed to predict DDIs and ADRs. Recent artificial intelligence (AI) approaches offer new opportunities for better predicting and understanding the complex processes related to drug metabolism and transport. We summarize the role of major DMEs and DTs, and provide an overview of current progress in computational approaches for the prediction of drug metabolism, transport, and DDIs, with an emphasis on AI including machine learning (ML) and deep learning (DL) modeling.

Unraveling the therapeutic potential of Senna singueana phytochemicals to attenuate pancreatic cancer using protein-protein interactions, molecular docking, and MD simulation

Pancreatic cancer (PC) presents challenges due to limited treatment options and patients seek complementary therapies alongside conventional treatments to improve well-being. This study uses computational drug discovery approaches to find potential phytochemicals from S. singueana for PC treatment. Among the 38 phytochemicals screened from S. singueana, specific inhibitors against PC were selected. Protein-protein interaction (PPI) network analysis highlighted key targets with high degrees, including PTEN (8) and PTK2 (7) genes, along with their respective proteins 5BZX and 3BZ3, which were employed for molecular docking studies. 1-methylchrysene and 3-methyl-1,8,9-anthracenetriol showed strong binding affinities of - 9.2 and - 8.1 Kcal/mol, respectively. Molecular dynamics simulations lasting 300 ns assessed structural stability and interaction energy of compound-target dockings: 1-methylchrysene-PTEN and 3-methyl-1,8,9-anthracenetriol-PTK2. In molecular dynamics simulations, the 3-methyl-1,8,9-anthracenetriol-PTK2 complex showed lower RMSD, RMSF, radius of gyration, solvent-accessible surface area, and more hydrogen bonds than the 1-methylchrysene-PTEN complex. The 3-methyl-1,8,9-anthracenetriol-PTK2 complex exhibited significantly stronger binding with a binding free energy (ΔG) of - 21.92 kcal/mol compared to the less favourable ΔG of - 10.65 kcal/mol for the 1-methylchrysene-PTEN complex. The consistent and stable binding interaction observed in the 3-methyl-1,8,9-anthracenetriol-PTK2 complex highlights its potential as a potent inhibitor of Focal Adhesion Kinase 1. Consequently, it emerges as a promising lead compound for the development of pancreatic cancer therapeutics. Conversely, the fluctuations observed in the 1-methylchrysene-PTEN complex indicate a less stable binding interaction. This indicates the potential of 3-methyl-1,8,9-anthracenetriol as a primary candidate for pancreatic cancer treatment. These findings improve our grasp of S. singueana's multi-target effects and its promise in addressing pancreatic cancer. Nevertheless, additional in-vivo and in-vitro studies are required to validate their effectiveness and therapeutic potential.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40203-023-00179-9.

Identification of novel pyrazole containing ɑ-glucosidase inhibitors: insight into pharmacophore, 3D-QSAR, virtual screening, and molecular dynamics study

Pharmacophore modelling, 3 D QSAR modelling, virtual screening, and molecular dynamics study, all-in-one combination were employed successfully design and develop an alpha-glucosidase inhibitor. To explain the structural prerequisites of biologically active components, 3 D-QSAR models were generated using the selected best hypothesis (AARRR) for compounds 55 included in the model C. The selection of 3 D-QSAR models showed that the Gaussian steric characteristic is crucial to alpha glucosidase's inhibitory potential. The alpha-glucosidase inhibitory potency of the compound is enhanced by other components, including Gaussian hydrophobic groups, Gaussian hydrogen bond acceptor or donor groups, Gaussian electrostatic characteristics, and a Gaussian steric feature. An identification of structure-activity relationships can be obtained from the developed 3 D-QSAR, C model, with R2 = 0.77 and SD = 0.02 for training set, and Q2 = 0.66, RMSE 0.02, and Pearson R = 0.81 for testing set, corresponding to elevated predictive ability. Additionally, docking and MM/GBSA experiments on 1146023 showed that it interacts with critical amino acids in the binding site when coupled with acarbose. Further, five compounds that display a high affinity for alpha-glucosidase were found, and these compounds may serve as potent leads for alpha-glucosidase inhibitor development. Biological activity will be tested for these compounds in the future.Communicated by Ramaswamy H. Sarma.

In silico drug repurposing and lipid bilayer molecular dynamics puzzled out potential breast cancer resistance protein (BCRP/ABCG2) inhibitors

Multidrug resistance (MDR) is a fundamental reason for the fiasco of carcinoma chemotherapy. A wide variety of anticarcinoma drugs are expelled from neoplasm cells through the ATP-binding cassette (ABC) transporter superfamily, rendering the neoplasm cells resistant to treatment. The ATP-binding cassette transporter G2 (ABCG2, gene symbol BCRP) is an ABC efflux transporter that plays a key function in MDR to antineoplastic therapies. For these reasons, the identification of medicaments as BCRP inhibitors could assist in discovering better curative approaches for breast cancer therapy. Because of the deficiency of prospective BCRP inhibitors, the SuperDRUG2 database was virtually screened for inhibitor activity towards the BCRP transporter using molecular docking computations. The most potent drug candidates were then characterized utilizing molecular dynamics (MD) simulations. Furthermore, molecular mechanics-generalized Born surface area (MM-GBSA) binding affinities of the most potent drug candidates were estimated. Based on the MM-GBSA binding affinities throughout 150 ns MD simulations, three drugs-namely zotarolimus (SD002595), temsirolimus (SD003393), and glecaprevir (SD006009)-revealed greater binding affinities towards BCRP transporter compared to the co-crystallized BWQ ligand with ΔGbinding values of -86.6 ± 5.6, -79.5 ± 8.0, -75.8 ± 4.6 and -59.5 ± 4.1 kcal/mol, respectively. The steadiness of these promising drugs bound with BCRP transporter was examined utilizing their structural and energetical analyses throughout a 150 ns MD simulation. To imitate the physiological environment, 150 ns MD simulations for the identified drugs bound with BCRP transporter were conducted in the 1-palmitoyl-2-oleoyl-phosphatidylcholine lipid bilayer. These findings identify zotarolimus, temsirolimus and glecaprevir as auspicious anti-MDR drug leads that warrant further experimental assays.Communicated by Ramaswamy H. Sarma.

SuperNatural inhibitors to reverse multidrug resistance emerged by ABCB1 transporter: Database mining, lipid-mediated molecular dynamics, and pharmacokinetics study

An effective approach to reverse multidrug resistance (MDR) is P-glycoprotein (P-gp, ABCB1) transport inhibition. To identify such molecular regulators, the SuperNatural II database, which comprises > 326,000 compounds, was virtually screened for ABCB1 transporter inhibitors. The Lipinski rule was utilized to initially screen the SuperNatural II database, identifying 128,126 compounds. Those natural compounds were docked against the ABCB1 transporter, and those with docking scores less than zosuquidar (ZQU) inhibitor were subjected to molecular dynamics (MD) simulations. Based on MM-GBA binding energy (ΔGbinding) estimations, UMHSN00009999 and UMHSN00097206 demonstrated ΔGbinding values of -68.3 and -64.1 kcal/mol, respectively, compared to ZQU with a ΔGbinding value of -49.8 kcal/mol. For an investigation of stability, structural and energetic analyses for UMHSN00009999- and UMHSN00097206-ABCB1 complexes were performed and proved the high steadiness of these complexes throughout 100 ns MD simulations. Pharmacokinetic properties of the identified compounds were also predicted. To mimic the physiological conditions, MD simulations in POPC membrane surroundings were applied to the UMHSN00009999- and UMHSN00097206-ABCB1 complexes. These results demonstrated that UMHSN00009999 and UMHSN00097206 are promising ABCB1 inhibitors for reversing MDR in cancer and warrant additional in-vitro/in-vivo studies.

In-Silico Mining of the Toxins Database (T3DB) towards Hunting Prospective Candidates as ABCB1 Inhibitors: Integrated Molecular Docking and Lipid Bilayer-Enhanced Molecular Dynamics Study

Multidrug resistance (MDR) is one of the most problematic issues in chemotherapeutic carcinoma therapy. The ABCB1 transporter, a drug efflux pump overexpressed in cancer cells, has been thoroughly investigated for its association with MDR. Thus, discovering ABCB1 inhibitors can reverse the MDR in cancer cells. In the current work, a molecular docking technique was utilized for hunting the most prospective ABCB1 inhibitors from the Toxin and Toxin-Target Database (T3DB). Based on the docking computations, the most promising T3DB compounds complexed with the ABCB1 transporter were subjected to molecular dynamics (MD) simulations over 100 ns. Utilizing the MM-GBSA approach, the corresponding binding affinities were computed. Compared to ZQU (calc. -49.8 kcal/mol), Emamectin B1a (T3D1043), Emamectin B1b (T3D1044), Vincristine (T3D4016), Vinblastine (T3D4017), and Vindesine (T3D2479) complexed with ABCB1 transporter demonstrated outstanding binding affinities with ΔG_binding values of -93.0, -92.6, -93.8, -92.2, and -90.8 kcal/mol, respectively. The structural and energetic investigations confirmed the constancy of the identified T3DB compounds complexed with the ABCB1 transporter during the 100 ns MD course. To mimic the physiological conditions, MD simulations were conducted for those identified inhibitors complexed with ABCB1 transporter in the presence of a POPC membrane. These findings revealed that Emamectin B1a, Emamectin B1b, Vincristine, Vinblastine, and Vindesine are promising ABCB1 inhibitors that can reverse the MDR. Therefore, subjecting those compounds to further in-vitro and in-vivo investigations is worthwhile.

Exploring Probenecid Derived 1,3,4-Oxadiazole-Phthalimide Hybrid as α-Amylase Inhibitor: Synthesis, Structural Investigation, and Molecular Modeling

1,3,4-Oxadiazole moiety is a crucial pharmacophore in many biologically active compounds. In a typical synthesis, probenecid was subjected to a sequence of reactions to obtain a 1,3,4-oxadiazole–phthalimide hybrid (PESMP) in high yields. The NMR (1H and 13C) spectroscopic analysis initially confirmed the structure of PESMP. Further spectral aspects were validated based on a single-crystal XRD analysis. Experimental findings were confirmed afterwards by executing a Hirshfeld surface (HS) analysis and quantum mechanical computations. The HS analysis showed the role of the π⋯π stacking interactions in PESMP. PESMP was found to have a high stability and lower reactivity in terms of global reactivity parameters. α-Amylase inhibition studies revealed that the PESMP was a good inhibitor of α-amylase with an s value of 10.60 ± 0.16 μg/mL compared with that of standard acarbose (IC50 = 8.80 ± 0.21 μg/mL). Molecular docking was also utilized to reveal the binding pose and features of PESMP against the α-amylase enzyme. Via docking computations, the high potency of PESMP and acarbose towards the α-amylase enzyme was unveiled and confirmed by docking scores of −7.4 and −9.4 kcal/mol, respectively. These findings shine a new light on the potential of PESMP compounds as α-amylase inhibitors.

Potential drug candidates as P-glycoprotein inhibitors to reverse multidrug resistance in cancer: an in silico drug discovery study

The failure of chemotherapy in the treatment of carcinoma is mainly due to the development of multidrug resistance (MDR), which is largely caused by the overexpression of P-glycoprotein (P-gp/ABCB1/MDR1). Until recently, the 3D structure of the P-gp transporter has not been experimentally resolved, which restricted the discovery of prospective P-gp inhibitors utilizing in silico techniques. In this study, the binding energies of 512 drug candidates in clinical or investigational stages were assessed as potential P-gp inhibitors employing in silico methods. On the basis of the available experimental data, the performance of the AutoDock4.2.6 software to predict the drug-P-gp binding mode was initially validated. Molecular docking and molecular dynamics (MD) simulations combined with molecular mechanics-generalized Born surface area (MM-GBSA) binding energy computations were subsequently conducted to screen the investigated drug candidates. Based on the current results, five promising drug candidates, namely valspodar, dactinomycin, elbasvir, temsirolimus, and sirolimus, showed promising binding energies against P-gp transporter with ΔGbinding values of -126.7, -112.1, -111.9, -102.9, and -101.4 kcal/mol, respectively. The post-MD analyses revealed the energetical and structural stabilities of the identified drug candidates in complex with the P-gp transporter. Furthermore, in order to mimic the physiological conditions, the potent drugs complexed with the P-gp were subjected to 100 ns MD simulations in an explicit membrane-water environment. The pharmacokinetic properties of the identified drugs were predicted and demonstrated good ADMET characteristics. Overall, these results indicated that valspodar, dactinomycin, elbasvir, temsirolimus, and sirolimus hold promise as prospective P-gp inhibitors and warrant further invitro/invivo investigations.

Computer-based identification of olive oil components as a potential inhibitor of neirisaral adhesion a regulatory protein

In silico methods such as molecular docking and molecular dynamic (MD) simulations have significant interest due to their ability to identify the protein-ligand interactions at the atomic level. In this work, different computational methods were used to elucidate the ability of some olive oil components to act as Neisseria adhesion A Regulatory protein (NadR) inhibitors. The frontier molecular orbitals (FMOs) and the global properties such as global hardness, electronegativity, and global softness of ten olive oil components (α-Tocopherol, Erythrodiol, Hydroxytyrosol, Linoleic acid, Apigenin, Luteolin, Oleic acid, Oleocanthal, Palmitic acid, and Tyrosol) were reported using Density Functional Theory (DFT) methods. Among all investigated compounds, Erythrodiol, Apigenin, and Luteolin demonstrated the highest binding affinities (-8.72, -7.12, and -8.24 kcal/mol, respectively) against NadR, compared to -8.21 kcal/mol of the native ligand based on molecular docking calculations. ADMET properties and physicochemical features showed that Erythrodiol, Apigenin, and Luteolin have good physicochemical features and can act as drugs candidate. Molecular dynamics (MD) simulations demonstrated that Erythrodiol, Apigenin, and Luteolin show stable binding affinity and molecular interaction with NadR. Further Molecular Mechanics Poisson-Boltzmann Surface Area (MM-PBSA) analyses using the MD trajectories also demonstrated the higher binding affinity of Erythrodiol, Apigenin and Luteolin inside NadR protein. The overall study provides a rationale to use Erythrodiol, Apigenin, and Luteolin in the drug development as anti-adhesive drugs lead. Communicated by Ramaswamy H. Sarma.

Venetoclax analogs as promising anticancer therapeutics via targeting Bcl-2 protein: in-silico drug discovery study

B-cell lymphoma 2 (Bcl-2) protein plays a vital role in enhancing malignant cell survival by alleviating programmed cell death. Therefore, Bcl-2 protein has been identified as a charming druggable target for cancer treatment. Venetoclax has enticed considerable attention as a potential Bcl-2 inhibitor. Herein, in-silico computations were executed to search for new venetoclax analogs against the Bcl-2 protein. A library involving 4112 was collected, prepared, and virtually screened against Bcl-2 protein using AutoDock Vina1.1.2 software. Promising analogs in complex with Bcl-2 protein were further submitted to molecular dynamics (MD) simulations, pursued by binding energy computations using the MM-GBSA approach. Compared to venetoclax (ΔGbinding = -51.2 kcal/mol), PubChem-873-158-83 and PubChem-148-422-478 demonstrated greater binding affinities with Bcl-2 protein throughout 100 ns MD simulations with ΔGbinding values of -69.1 and -62.4 kcal/mol, respectively. Structural and energetical analyses unveiled good stabilization of the identified analogs complexed with Bcl-2 protein over the MD course. The pharmacokinetic features of the two identified analogs were anticipated and unveiled the oral bioavailability of these compounds. Further in-vitro/in-vivo biological evaluations around these compounds could assist in identifying anticancer leads towards Bcl-2 protein.Communicated by Ramaswamy H. Sarma.

Click-designed vanilloid-triazole conjugates as dual inhibitors of AChE and Aβ aggregation

Based on their reported neuroprotective properties, vanilloids provide a good starting point for the synthesis of anti-Alzheimer's disease (AD) agents. In this context, nine new 1,2,3-triazole conjugates of vanilloids were synthesized via click chemistry. The compounds were tested for their effect on acetylcholine esterase (AChE) and amyloid-beta peptide (Aβ) aggregation. The triazole esters (E)-(1-(4-hydroxy-3-methoxybenzyl)-1H-1,2,3-triazol-4-yl)methyl 3-(4-hydroxy-3 methoxyphenyl)acrylate 9 and (1-(4-hydroxy-3-methoxybenzyl)-1H-1,2,3-triazol-4-yl)methyl-4-hydroxy-3-methoxybenzoate 8 displayed dual inhibitory activity for AChE and Aβ aggregation with IC50 values of 0.47/0.31 μM and 1.2/0.95 μM, respectively, as compared to donepezil (0.27 μM) and tacrine (0.41 μM), respectively. The results showed that the triazole ester moiety is more favorable for the activity than the triazole ether moiety. This could be attributed to the longer length of the spacer between the two vanillyl moieties in the triazole esters. Furthermore, the binding affinities and modes of the triazole esters 9 and 8 were examined against AChE and Aβ utilizing a combination of docking predictions and molecular dynamics (MD) simulations. Docking computations revealed promising binding affinity of triazole esters 9 and 8 as potential AChE, Aβ40, and Aβ42 inhibitors with docking scores of -10.4 and -9.4 kcal mol-1, -5.8 and -4.7 kcal mol-1, and -3.3 and -2.9 kcal mol-1, respectively. The stability and binding energies of triazole esters 9 and 8 complexed with AChE, Aβ40, and Aβ42 were measured and compared to donepezil and tacrine over 100 ns MD simulations. According to the estimated binding energies, compounds 9 and 8 displayed good binding affinities with AChE, Aβ42, and Aβ40 with average ΔG binding values of -32.9 and -31.8 kcal mol-1, -12.0 and -10.5 kcal mol-1, and -20.4 and -16.6 kcal mol-1, respectively. Post-MD analyses demonstrated high steadiness for compounds 9 and 8 with AChE and Aβ during the 100 ns MD course. This work suggests the triazole conjugate of vanilloids as a promising skeleton for developing multi-target potential AD therapeutics.

Chetomin, a SARS-CoV-2 3C-like Protease (3CLpro) Inhibitor: In Silico Screening, Enzyme Docking, Molecular Dynamics and Pharmacokinetics Analysis

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to over 6 million deaths. The 3C-like protease (3CLpro) enzyme of the SARS-CoV-2 virus is an attractive druggable target for exploring therapeutic drug candidates to combat COVID-19 due to its key function in viral replication. Marine natural products (MNPs) have attracted considerable attention as alternative sources of antiviral drug candidates. In looking for potential 3CLpro inhibitors, the MNP database (>14,000 molecules) was virtually screened against 3CLpro with the assistance of molecular docking computations. The performance of AutoDock and OEDocking software in anticipating the ligand-3CLpro binding mode was first validated according to the available experimental data. Based on the docking scores, the most potent MNPs were further subjected to molecular dynamics (MD) simulations, and the binding affinities of those molecules were computed using the MM-GBSA approach. According to MM-GBSA//200 ns MD simulations, chetomin (UMHMNP1403367) exhibited a higher binding affinity against 3CLpro than XF7, with ΔGbinding values of −55.5 and −43.7 kcal/mol, respectively. The steadiness and tightness of chetomin with 3CLpro were evaluated, revealing the high stabilization of chetomin (UMHMNP1403367) inside the binding pocket of 3CLpro throughout 200 ns MD simulations. The physicochemical and pharmacokinetic features of chetomin were also predicted, and the oral bioavailability of chetomin was demonstrated. Furthermore, the potentiality of chetomin analogues −namely, chetomin A-D− as 3CLpro inhibitors was investigated. These results warrant further in vivo and in vitro assays of chetomin (UMHMNP1403367) as a promising anti-COVID-19 drug candidate.

Design and synthesis of new thiazolidinone/uracil derivatives as antiproliferative agents targeting EGFR and/or BRAFV600E

Thiourea derivatives of uracil were efficiently synthesized via the reaction of 5-aminouracil with isothiocyanates. Then, we prepared uracil-containing thiazoles via condensation of thioureas with diethyl/dimethyl acetylenedicarboxylates. The structures of the products were confirmed by a combination of spectral techniques including infra-red (IR), nuclear magnetic resonance (NMR), mass spectrometry (MS) and elemental analyses. A rationale for the formation of the products is presented. The newly synthesized compounds were evaluated for their in vitro antiproliferative activity against four cancer cell lines. The compounds tested showed promising antiproliferative activity, with GI₅₀ values ranging from 1.10 µM to 10.00 µM. Compounds 3c, 5b, 5c, 5h, 5i, and 5j were the most potent derivatives, with GI₅₀ values ranging from 1.10 µM to 1.80 µM. Compound 5b showed potent inhibitory activity against EGFR and BRAF^V600E with IC₅₀ of 91 ± 07 and 93 ± 08 nM, respectively, indicating that this compound could serve as a dual inhibitor of EGFR and BRAF^V600E with promising antiproliferative properties. Docking computations revealed the great potency of compounds 5b and 5j towards EGFR and BRAF^V600E with docking scores of -8.3 and -9.7 kcal/mol and -8.2 and -9.3 kcal/mol, respectively.

Synthesis, X-ray diffraction analysis, quantum chemical studies and α-amylase inhibition of probenecid derived S-alkylphthalimide-oxadiazole-benzenesulfonamide hybrids

Sulphonamide and 1,3,4-oxadiazole moieties are present as integral structural parts of many drugs and pharmaceuticals. Taking into account the significance of these moieties, we herein present the synthesis, single-crystal X-ray analysis, DFT studies, and α-amylase inhibition of probenecid derived two S-alkylphthalimide-oxadiazole-benzenesulfonamide hybrids. The synthesis has been accomplished in high yields. The final structures of both hybrids have been established completely with the help of different spectro-analytical techniques, including NMR, FTIR, HR-MS, and single-crystal X-ray diffraction analyses. In an effort to confirm the experimental findings, versatile quantum mechanical calculations and Hirshfeld Surface analysis have been performed. α-Amylase inhibition assay has been executed to investigate the enzyme inhibitory potential of both hybrids. The low IC50 value (76.92 ± 0.19 μg/mL) of hybrid 2 shows the good α-amylase inhibition potential of the respective compound. Ultimately, the binding affinities and features of the two hybrids are elucidated utilising a molecular docking technique against the α-amylase enzyme.

Exploring Natural Product Activity and Species Source Candidates for Hunting ABCB1 Transporter Inhibitors: An In Silico Drug Discovery Study

The P-glycoprotein (P-gp/ABCB1) is responsible for a xenobiotic efflux pump that shackles intracellular drug accumulation. Additionally, it is included in the dud of considerable antiviral and anticancer chemotherapies because of the multidrug resistance (MDR) phenomenon. In the search for prospective anticancer drugs that inhibit the ABCB1 transporter, the Natural Product Activity and Species Source (NPASS) database, containing >35,000 molecules, was explored for identifying ABCB1 inhibitors. The performance of AutoDock4.2.6 software to anticipate ABCB1 docking score and pose was first assessed according to available experimental data. The docking scores of the NPASS molecules were predicted against the ABCB1 transporter. Molecular dynamics (MD) simulations were conducted for molecules with docking scores lower than taxol, a reference inhibitor, pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations. On the basis of MM-GBSA calculations, five compounds revealed promising binding affinities as ABCB1 inhibitors with ΔGbinding < −105.0 kcal/mol. The binding affinity and stability of the identified inhibitors were compared to the chemotherapeutic agent. Structural and energetical analyses unveiled great steadiness of the investigated inhibitors within the ABCB1 active site throughout 100 ns MD simulations. Conclusively, these findings point out that NPC104372, NPC475164, NPC2313, NPC197736, and NPC477344 hold guarantees as potential ABCB1 drug candidates and warrant further in vitro/in vivo tests.

Naturally occurring plant-based anticancerous candidates as prospective ABCG2 inhibitors: an in silico drug discovery study

Abstract

ATP-binding cassette transporter G2 (ABCG2) is an efflux transporter related to the clinical multidrug resistance (MDR) phenomenon. Identifying ABCG2 inhibitors could help discover extraordinary curative strategies for carcinoma remediation. Hitherto, there is no medication drug inhibiting ABCG2 transporter, notwithstanding that a considerable number of drugs have been submitted to clinical-trial and investigational phases. In the search for unprecedented chemical compounds that could inhibit the ABCG2 transporter, an in silico screening was conducted on the Naturally Occurring Plant-based Anticancer Compound-Activity-Target (NPACT) database containing 1574 compounds. Inhibitor-ABCG2 binding affinities were estimated based on molecular docking and molecular minimization (MM) calculations and compared to a co-crystallized inhibitor (BWQ) acting as a reference inhibitor. Molecular dynamics (MD) simulations pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations were further executed for compounds with MM-GBSA//MM binding energies lower than BWQ (calc. − 60.5 kcal/mol). NPACT00968 and NPACT01545 demonstrated auspicious inhibitory activities according to binding affinities (ΔGbinding) over the 100 ns MD simulations that were nearly one and a half folds compared to BWQ (− 100.4, − 94.7, and − 62.9 kcal/mol, respectively). Throughout the 100 ns MD simulations, structural and energetical analyses unveiled outstanding stability of the ABCG2 transporter when bound with NPACT00968 and NPACT01545. In silico calculations hold a promise for those two inhibitors as drug candidates of ABCG2 transporter and emphasize that further in vitro and in vivo experiments are guaranteed.

Graphical abstract

NS3 helicase inhibitory potential of the marine sponge Spongia irregularis

In the current study, an investigation of the activity of the total extract of the marine sponge Spongia irregularis and its different fractions against the hepatitis C virus (HCV) was pursued. The results revealed that the ethyl acetate fraction exhibited the highest anti-HCV activity, with an IC50 value of 12.6 ± 0.05 μg ml-1. Chromatographic resolution of the ethyl acetate fraction resulted in the isolation of four known compounds, 1,3,7-trimethylguanine (1), 3,5-dihydroxyfuran-2(5H)-one (2), thymidine (3), and 1H-indazole (4). By using LC-HR-ESI-MS metabolic profiling, compounds 5-14 were also identified in the same fraction. Molecular docking calculations revealed the high binding affinity of compound 14 against the allosteric pocket of HCV NS3-NS4A and the active site of HCV NS3 helicase (-10.1 and -7.4 kcal mol-1, respectively). Molecular dynamics simulations, followed by molecular mechanics-generalized Born surface area energy calculations, demonstrated the structural and energetic stability of compound 14 in complex with HCV targets.