In Silico Prediction of P-glycoprotein Binding: Insights from Molecular Docking Studies

Py-CoMFA, docking, and molecular dynamics simulations of Leishmania (L.) amazonensis arginase inhibitors

Leishmaniasis is a disease caused by a protozoan of the genus Leishmania, affecting millions of people, mainly in tropical countries, due to poor social conditions and low economic development. First-line chemotherapeutic agents involve highly toxic pentavalent antimonials, while treatment failure is mainly due to the emergence of drug-resistant strains. Leishmania arginase (ARG) enzyme is vital in pathogenicity and contributes to a higher infection rate, thus representing a potential drug target. This study helps in designing ARG inhibitors for the treatment of leishmaniasis. Py-CoMFA (3D-QSAR) models were constructed using 34 inhibitors from different chemical classes against ARG from L. (L.) amazonensis (LaARG). The 3D-QSAR predictions showed an excellent correlation between experimental and calculated pIC50 values. The molecular docking study identified the favorable hydrophobicity contribution of phenyl and cyclohexyl groups as substituents in the enzyme allosteric site. Molecular dynamics simulations of selected protein-ligand complexes were conducted to understand derivatives' interaction modes and affinity in both active and allosteric sites. Two cinnamide compounds, 7g and 7k, were identified, with similar structures to the reference 4h allosteric site inhibitor. These compounds can guide the development of more effective arginase inhibitors as potential antileishmanial drugs.

ADMET-PrInt: Evaluation of ADMET Properties: Prediction and Interpretation

Great progress in the development of computational strategies for drug design applications has revolutionized the process of searching for new drugs. Although the focus of in silico strategies is still put on the provision of the desired activity of a compound to the considered target, characterization of a compound in terms of its physicochemical and ADMET properties becomes an indispensable element of computer-aided drug design protocols. In the study, an online application ADMET-PrInt for in silico assessment of selected compound features: cardiotoxicity, solubility, genotoxicity, membrane permeability, and plasma protein binding was prepared. In addition to the prediction of particular property, ADMET-PrInt enables also the identification of compound features influencing this property thanks to the application of two explainability approaches: local interpretabile model-agnostic explanations and counterfactual analysis. It is an important factor for medicinal chemists, as it greatly facilitates the process of optimization of the compound structure in terms of the evaluated properties. The intuitive webpage, available at admet.if-pan.krakow.pl, allows making use of all predictive and interpretability models also by nonexperts and nonprogrammers.

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.

Advances in Developing Small Molecule Drugs for Alzheimer's Disease

Alzheimer's disease (AD) is the most common type of dementia among middle-aged and elderly individuals. Accelerating the prevention and treatment of AD has become an urgent problem. New technology including Computer-aided drug design (CADD) can effectively reduce the medication cost for patients with AD, reduce the cost of living, and improve the quality of life of patients, providing new ideas for treating AD. This paper reviews the pathogenesis of AD, the latest developments in CADD and other small-molecule docking technologies for drug discovery and development; the current research status of small-molecule compounds for AD at home and abroad from the perspective of drug action targets; the future of AD drug development.

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.

In silico approaches and in vitro assays identify a coumarin derivative as antiviral potential against SARS-CoV-2

COVID-19, a disease caused by SARS-CoV-2, was declared a pandemic in 2020 and created a global crisis in health systems, with more than 545 million confirmed cases and 6.33 million deaths. In this sense, this work aims to identify possible inhibitors of the SARS-CoV-2 RdRp enzyme using in silico approaches. RdRp is a crucial enzyme in the replication and assembly cycle of new viral particles and a critical pharmacological target in the treatment of COVID-19. We performed a virtual screening based on molecular docking from our in-house chemical library, which contains a diversity of 313 structures from different chemical classes. Nine compounds were selected since they showed important interactions with the active site from RdRp. Next, the ADME-Tox in silico predictions served as a filter and selected the three most promising compounds: a coumarin LMed-052, a hydantoin LMed-087, and a guanidine LMed-250. Molecular dynamics simulations revealed details such as changes in the positions of ligands and catalytic residues during the simulations compared to the complex from molecular docking studies. Binding free energy analysis was performed using the MMGBSA method, demonstrating that LMed-052 and LMed-087 have better affinities for the RdRp by energetic contributions to the stability of the complexes when compared to LMed-250. Furthermore, LMed-052 showed significant in vitro inhibition against MHV-3, decreasing 99% of viral titers. Finally, these findings are useful to guide structural modifications aiming to improve the potential of these compounds to act as inhibitors of SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Simulated Microgravity Alters P-Glycoprotein Efflux Function and Expression via the Wnt/β-Catenin Signaling Pathway in Rat Intestine and Brain

The drug efflux transporter permeability glycoprotein (P-gp) plays an important role in oral drug absorption and distribution. Under microgravity (MG), the changes in P-gp efflux function may alter the efficacy of oral drugs or lead to unexpected effects. Oral drugs are currently used to protect and treat multisystem physiological damage caused by MG; whether P-gp efflux function changes under MG remains unclear. This study aimed to investigate the alteration of P-gp efflux function, expression, and potential signaling pathway in rats and cells under different simulated MG (SMG) duration. The altered P-gp efflux function was verified by the in vivo intestinal perfusion and the brain distribution of P-gp substrate drugs. Results showed that the efflux function of P-gp was inhibited in the 7 and 21 day SMG-treated rat intestine and brain and 72 h SMG-treated human colon adenocarcinoma cells and human cerebral microvascular endothelial cells. P-gp protein and gene expression levels were continually down-regulated in rat intestine and up-regulated in rat brain by SMG. P-gp expression was regulated by the Wnt/β-catenin signaling pathway under SMG, verified by a pathway-specific agonist and inhibitor. The elevated intestinal absorption and brain distribution of acetaminophen levels also confirmed the inhibited P-gp efflux function in rat intestine and brain under SMG. This study revealed that SMG alters the efflux function of P-gp and regulates the Wnt/β-catenin signaling pathway in the intestine and the brain. These findings may be helpful in guiding the use of P-gp substrate drugs during spaceflight.

Antileishmanial Activity of 4,8-Dimethoxynaphthalenyl Chalcones on Leishmania amazonensis

Leishmaniasis is a neglected tropical disease caused by Leishmania species. Available therapeutic options have several limitations. The drive to develop new, more potent, and selective antileishmanial agents is thus a major goal. Herein we report the synthesis and the biological activity evaluation against promastigote and amastigote forms of Leishmania amazonensis of nine 4,8-dimethoxynaphthalenyl chalcones. Compound ((E)-1-(4,8-dimethoxynaphthalen-1-yl)-3-(4-nitrophenyl)prop-2-en-1-one), 4f, was the most promising with an IC₅₀ = 3.3 ± 0.34 μM (promastigotes), a low cytotoxicity profile (CC₅₀ = 372.9 ± 0.04 μM), and a high selectivity index (SI = 112.6). Furthermore, 4f induced several morphological and ultrastructural changes in the free promastigote forms, loss of plasma membrane integrity, and increased reactive oxygen species (ROS). An in silico analysis of drug-likeness and ADME parameters suggested high oral bioavailability and intestinal absorption. Compound 4f reduced the number of infected macrophages and the number of amastigotes per macrophage, with an IC₅₀ value of 18.5 ± 1.19 μM. Molecular docking studies with targets, ARG and TR, showed that compound 4f had more hydrogen bond interactions with the ARG enzyme, indicating a more stable protein-ligand binding. These results suggest that 4,8-dimethoxynaphthalenyl chalcones are worthy of further study as potential antileishmanial drugs.

Prediction of potential drug interactions between repurposed COVID-19 and antitubercular drugs: an integrational approach of drug information software and computational techniques data

Introduction:

Tuberculosis is a major respiratory disease globally with a higher prevalence in Asian and African countries than rest of the world. With a larger population of tuberculosis patients anticipated to be co-infected with COVID-19 infection, an ongoing pandemic, identifying, preventing and managing drug–drug interactions is inevitable for maximizing patient benefits for the current repurposed COVID-19 and antitubercular drugs.

Methods:

We assessed the potential drug–drug interactions between repurposed COVID-19 drugs and antitubercular drugs using the drug interaction checker of IBM Micromedex®. Extensive computational studies were performed at a molecular level to validate and understand the drug–drug interactions found from the Micromedex drug interaction checker database at a molecular level. The integrated knowledge derived from Micromedex and computational data was collated and curated for predicting potential drug–drug interactions between repurposed COVID-19 and antitubercular drugs.

Results:

A total of 91 potential drug–drug interactions along with their severity and level of documentation were identified from Micromedex between repurposed COVID-19 drugs and antitubercular drugs. We identified 47 pharmacodynamic, 42 pharmacokinetic and 2 unknown DDIs. The majority of our molecular modelling results were in line with drug–drug interaction data obtained from the drug information software. QT prolongation was identified as the most common type of pharmacodynamic drug–drug interaction, whereas drug–drug interactions associated with cytochrome P450 3A4 (CYP3A4) and P-glycoprotein (P-gp) inhibition and induction were identified as the frequent pharmacokinetic drug–drug interactions. The results suggest antitubercular drugs, particularly rifampin and second-line agents, warrant high alert and monitoring while prescribing with the repurposed COVID-19 drugs.

Conclusion:

Predicting these potential drug–drug interactions, particularly related to CYP3A4, P-gp and the human Ether-à-go-go-Related Gene proteins, could be used in clinical settings for screening and management of drug–drug interactions for delivering safer chemotherapeutic tuberculosis and COVID-19 care. The current study provides an initial propulsion for further well-designed pharmacokinetic-pharmacodynamic-based drug–drug interaction studies.

Plain Language Summary

Comparative in silico prediction of P-glycoprotein-mediated transport for 2010-2020 US FDA-approved drugs using six Web-tools

P-glycoprotein (P-gp) is an efflux pump implicated in pharmacokinetics and drug-drug interactions. The identification of its substrates is consequently an important issue, notably for drugs under development. For such a purpose, various in silico methods have been developed, but their relevance remains to be fully established. The present study was designed to get insight about this point, through determining the performance values of six freely accessible Web-tools (ADMETlab, AdmetSAR2.0, PgpRules, pkCSM, SwissADME and vNN-ADMET), computationally predicting P-gp-mediated transport. Using an external test set of 231 marketed drugs, approved over the 2010-2020 period by the US Food and Drug Administration and fully in vitro characterized for their P-gp substrate status, various performance parameters (including sensitivity, specificity, accuracy, Matthews correlation coefficient and area under the receiver operating characteristics curve) were determined. They were found to rather poorly meet criteria commonly required for acceptable prediction, whatever the Web-tools were used alone or in combination. Predictions of being P-gp substrate or non-substrate by these online in silico methods may therefore be considered with caution.

Unique indolizidine alkaloid securinine is a promising scaffold for the development of neuroprotective and antitumor drugs

Alkaloids, secondary plant metabolites, are used in traditional medicine in many countries to treat various pathological conditions. Securinine, a unique indolizidine alkaloid combining four cycles, "6-azobicyclo[3.2.1]octane" as a key structure fused with α,β-unsaturated-γ-lactone and piperidine ring, has a broad spectrum of actions including anti-inflammatory, antibacterial, neuroprotective and antitumor, and has been previously used in medical practice. It has several reactive centers, which are double bonds at positions 12-13 and 14-15, and this is a challenging scaffold for the synthesis of biologically active compounds. In this review, works on the production of modified securinine derivatives and their biological activity are addressed. Both monovalent and bivalent derivatives that are most promising in our opinion, and have potential for further research, are considered.

Advances in Applying Computer-Aided Drug Design for Neurodegenerative Diseases

Neurodegenerative diseases (NDs) including Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, and Huntington's disease are incurable and affect millions of people worldwide. The development of treatments for this unmet clinical need is a major global research challenge. Computer-aided drug design (CADD) methods minimize the huge number of ligands that could be screened in biological assays, reducing the cost, time, and effort required to develop new drugs. In this review, we provide an introduction to CADD and examine the progress in applying CADD and other molecular docking studies to NDs. We provide an updated overview of potential therapeutic targets for various NDs and discuss some of the advantages and disadvantages of these tools.

Surface plasmon resonance biosensor combined with lentiviral particle stabilization strategy for rapid and specific screening of P-Glycoprotein ligands

A novel surface plasmon resonance-based P-gp ligand screening system (SPR-PLSS) combined with lentiviral particle (LVP) stabilization strategy was constructed to screen out potential P-gp inhibitors from natural products. Firstly, we constructed LVPs with high and low expression levels of P-gp. The LVPs can ensure the natural conformation of P-gp based on the principle that LVPs germinated from packaging cells will contain cell membrane fragments and P-gp they carry. Then the LVPs with high P-gp expression for active channel and LVPs with low P-gp expression for reference channel were immobilized on CM5 chip respectively. The affinity detection was thus carried out with the signal reduction on the two channels. The P-gp inhibitors, Valspodar (Val) and cyclosporin (CsA), as positive compounds, were detected to characterize the chip's activity, and the K_D of Val and CsA were 14.09 μM and 16.41 μM, respectively. Forty compounds from natural product library were screened using the SPR CM5 chip, and magnolol (Mag), honokiol (Hon), and resveratrol (Res) were screened out as potential P-gp ligands, showing a significant response signal. This work presented a novel P-gp ligand screening system based on LVP-immobilized biosensor to rapidly screen out P-gp ligands from natural product library. Compared with traditional cell experiments which the screening time may take up to several days, our method only takes several hours. Furthermore, this study has also provided solid evidences to support that some complicated membrane proteins would apply to the lentivirus-based SPR screening system.

Synthesis, Molecular Docking and Molecular Dynamics Simulation of 2- Thioxothiazolidin-4-One Derivatives against Gp41

INTRODUCTION

Gp41 and its conserved hydrophobic groove on the N-terminal heptad repeat region are attractive targets in the design of HIV-1 entry inhibitors. Linearly extended molecules have shown potent anti-HIV-1 activity for their effective interactions with the gp41 binding pocket. Rhodanine ring attached to substituted pyrrole or furan rings has been proved a preferred moiety to be inserted inside the molecular structure of the gp41 inhibitors.

OBJECTIVES

Based on the previous findings we are going to describe some rhodanine derivatives in which a substituted imidazole ring is introduced in place of the pyrrole or furan rings. The compounds' flexibility is increased by inserting methylene groups inside the main scaffold.

METHODS

Molecular docking and molecular dynamics simulations approaches were exploited to investigate the chemical interactions and the stability of the designed ligands-gp41 complex. All compounds were synthesized and their chemical structures were elucidated by 1HNMR, 13CNMR, FTIR and Mass spectroscopy. Biological activities of the compounds against HIV-1 and HIV-2 and their cellular toxicities against the T-lymphocyte (MT-4) cell line were determined.

RESULTS

All the designed compounds showed proper and stable chemical interactions with gp41 according to the in silico studies. The results of the biological tests proved none of the compounds active against HIV-1 replication in cell cultures.

CONCLUSION

Since all the studied compounds were potently toxic for the host cell; it was therefore not possible to assess their anti-HIV activities.

In silico structural homology modelling of EST073 motif coding protein of tea Camellia sinensis (L)

Background

Tea (Camellia sinensis (L). O. Kuntze) is known as the oldest, mild stimulating caffeine containing non-alcoholic beverage. One of the major threats in south Asian tea industry is the blister blight leaf disease (BB), caused by the fungus Exobasidium vexans Masse. SSR DNA marker EST SSR 073 is used as a molecular marker to tag blister blight disease resistance trait of tea. The amino acid sequences were derived from cDNA sequences related to EST SSR 073 of BB susceptible (TRI 2023) and BB resistant (TRI 2043) cultivars. An attempt has been made to understand the structural characteristics and variations of EST SSR 073 locus that may reveal the factors influencing the BB resistance of tea with multiple bioinformatics tools such as ORF finder, ExPasy ProtParam tools, modeler V 9.17, Rampage server, UCSF-Chimera, and HADDOCK docking server.

Results

The primary, secondary, and tertiary structures of EST SSR 073 coding protein were analyzed using the amino acid sequences of both BB resistant TRI 2043 and BB susceptible TRI 2023 tea cultivars. The coding amino acid sequences of both the cultivars were homologous to photosystem I subunit protein (PsaD I) of Pisum sativum. The predicted 3D structures of proteins were validated and considered as an acceptable overall stereochemical quality. The BB resistant protein showed CT repeat extension and did not involve in topology of the PsaD I subunit. The C terminal truncation of BB resistance caused the formation of hydrogen bonds interacting with PsaD I and other subunits of photosystem I in the modeled three-dimensional protein structure.

Conclusions

Camellia sinensis EST 073 SSR motif coding protein was identified as the PsaD I subunit of photosystem I. The exact mechanism of PsaD I conferring the resistance for blister blight in tea needs to be further investigated.

Search for ABCB1 Modulators Among 2-Amine-5-Arylideneimidazolones as a New Perspective to Overcome Cancer Multidrug Resistance

Multidrug resistance (MDR) is a severe problem in the treatment of cancer with overexpression of glycoprotein P (Pgp, ABCB1) as a reason for chemotherapy failure. A series of 14 novel 5-arylideneimidazolone derivatives containing the morpholine moiety, with respect to two different topologies (groups A and B), were designed and obtained in a three- or four-step synthesis, involving the Dimroth rearrangement. The new compounds were tested for their inhibition of the ABCB1 efflux pump in both sensitive (parental (PAR)) and ABCB1-overexpressing (MDR) T-lymphoma cancer cells in a rhodamine 123 accumulation assay. Their cytotoxic and antiproliferative effects were investigated by a thiazolyl blue tetrazolium bromide (MTT) assay. For active compounds, an insight into the mechanisms of action using either the luminescent Pgp-Glo™ Assay in vitro or docking studies to human Pgp was performed. The safety profile in vitro was examined. Structure–activity relationship (SAR) analysis was discussed. The most active compounds, representing both 2-substituted- (11) and Dimroth-rearranged 3-substituted (18) imidazolone topologies, displayed 1.38–1.46 fold stronger efflux pump inhibiting effects than reference verapamil and were significantly safer than doxorubicin in cell-based toxicity assays in the HEK-293 cell line. Results of mechanistic studies indicate that active imidazolones are substrates with increasing Pgp ATPase activity, and their dye-efflux inhibition via competitive action on the Pgp verapamil binding site was predicted in silico.

Green and cytocompatible carboxyl modified gold-lysozyme nanoantibacterial for combating multidrug-resistant superbugs

The dissemination of multi-drug resistant (MDR) superbugs in hospital environments, communities and food animals and the very dynamic bacterial mutation frequency require the development of prolonged therapeutic strategies to gain mastery over antibiotic resistance. A AuNP-lysozyme nanoantibacterial was fabricated by the conjugation of AuNPs-C6H4-4-COOH with lysozyme via green reduction of aryldiazonium gold(iii) salt [HOOC-4-C6H4N[triple bond, length as m-dash]N]AuCl4. Results from molecular docking calculations aimed at revealing the binding mode of benzoic acid with the lysozyme structure clearly showed the lowest energy conformation with benzoic acid bound in the deep buried hydrophobic cavity of the protein active site through strong hydrogen bonding and hydrophobic interactions, thus validating the experimental outcomes of the current study which also exhibited the binding of -COOH functional groups in the interior of the protein structure. The superiority of the lysozyme bioconjugate against superbugs was demonstrated by the enhanced and broadened lysozyme antibacterial activities of 98-99% against extended spectrum beta lactamase (ESBL) producing Escherichia coli and imipenem-resistant Pseudomonas aeruginosa clinical isolates and a selection of Gram-negative and Gram-positive standard ATCC strains. Selective toxicity against bacteria was confirmed by the high viability of HeLa and fibroblast cell lines and the outstanding hemocompatibility at the minimum bacterial inhibitory concentrations (MICs). Turbidimetric enzyme kinetic assay showed the enhancement of the lysozyme hydrolytic activity by gold nanoparticles on the Micrococcus lysodeikticus bacterial substrate. Using gel electrophoresis, the induced cell wall breakdown was confirmed by detecting the leaked-out bacterial genomic DNA. The integrity and morphology changes of the E. coli bacteria were investigated using a scanning electron microscope after one hour of contact with the lysozyme-gold bioconjugate. The antibacterial functionalities showed little or no damage to healthy human cells and can be applied to wound dressings and medical devices.

Replacing the eleven native tryptophans by directed evolution produces an active P-glycoprotein with site-specific, non-conservative substitutions

P-glycoprotein (Pgp) pumps an array of hydrophobic compounds out of cells, and has major roles in drug pharmacokinetics and cancer multidrug resistance. Yet, polyspecific drug binding and ATP hydrolysis-driven drug export in Pgp are poorly understood. Fluorescence spectroscopy using tryptophans (Trp) inserted at strategic positions is an important tool to study ligand binding. In Pgp, this method will require removal of 11 endogenous Trps, including highly conserved Trps that may be important for function, protein-lipid interactions, and/or protein stability. Here, we developed a directed evolutionary approach to first replace all eight transmembrane Trps and select for transport-active mutants in Saccharomyces cerevisiae. Surprisingly, many Trp positions contained non-conservative substitutions that supported in vivo activity, and were preferred over aromatic amino acids. The most active construct, W(3Cyto), served for directed evolution of the three cytoplasmic Trps, where two positions revealed strong functional bias towards tyrosine. W(3Cyto) and Trp-less Pgp retained wild-type-like protein expression, localization and transport function, and purified proteins retained drug stimulation of ATP hydrolysis and drug binding affinities. The data indicate preferred Trp substitutions specific to the local context, often dictated by protein structural requirements and/or membrane lipid interactions, and these new insights will offer guidance for membrane protein engineering.