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 (
P-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 (Pgp) is a member of the ATP-binding cassette family of transporters that confers multidrug resistance to cancer cells and is actively involved in the pharmacokinetics and ...toxicokinetics of a big variety of drugs. Extensive studies have provided insights into the binding of many compounds, but the precise mechanism of translocation across the membrane remains unknown; in this context, the major challenge has been to understand the basis for its polyspecificity. In this study, molecular dynamics (MD) simulations of human P-gp (hP-gp) in an explicit membrane-and-water environment were performed to investigate the dynamic behavior of the transporter in the presence of different compounds (active and inactive) in the binding pocket and ATP molecules within the nucleotide binding domains (NBDs). The complexes studied involve four compounds: cyclosporin A (CSA), amiodarone (AMI), pamidronate (APD), and valproic acid (VPA). While CSA and AMI are known to interact with P-gp, APD and VPA do not. The results highlighted how the presence of ATP notably contributed to increased flexibility of key residues in NBD1 of active systems, indicating potential conformational changes activating the translocation mechanism. MD simulations reveal how these domains adapt and respond to the presence of different substrates, as well as the influence of ATP binding on their flexibility. Furthermore, distinctive behavior was observed in the presence of active and inactive compounds, particularly in the arrangement of ATP between NBDs, supporting the proposed nucleotide sandwich dimer mechanism for ATP binding. This study provides comprehensive insights into P-gp behavior with various ligands and ATP, offering implications for drug development, toxicity assessment and demonstrating the validity of the results derived from the MD simulations.
The emergence of SARS‐CoV‐2 has resulted in nearly 1,280,000 infections and 73,000 deaths globally so far. This novel virus acquired the ability to infect human cells using the SARS‐CoV cell receptor ...hACE2. Because of this, it is essential to improve our understanding of the evolutionary dynamics surrounding the SARS‐CoV‐2 hACE2 interaction. One way theory predicts selection pressures should shape viral evolution is to enhance binding with host cells. We first assessed evolutionary dynamics in select betacoronavirus spike protein genes to predict whether these genomic regions are under directional or purifying selection between divergent viral lineages, at various scales of relatedness. With this analysis, we determine a region inside the receptor‐binding domain with putative sites under positive selection interspersed among highly conserved sites, which are implicated in structural stability of the viral spike protein and its union with human receptor ACE2. Next, to gain further insights into factors associated with recognition of the human host receptor, we performed modeling studies of five different betacoronaviruses and their potential binding to hACE2. Modeling results indicate that interfering with the salt bridges at hot spot 353 could be an effective strategy for inhibiting binding, and hence for the prevention of SARS‐CoV‐2 infections. We also propose that a glycine residue at the receptor‐binding domain of the spike glycoprotein can have a critical role in permitting bat SARS‐related coronaviruses to infect human cells.
New series of 6-substituted-3-arylcoumarins displaying several alkyl, hydroxyl, halogen, and alkoxy groups in the two benzene rings have been designed, synthesized, and evaluated in vitro as human ...monoamine oxidase A and B (hMAO-A and hMAO-B) inhibitors. Most of the studied compounds showed a high affinity and selectivity to the hMAO-B isoenzyme, with IC50 values on nanomolar and picomolar range. Ten of the 22 described compounds displayed higher MAO-B inhibitory activity and selectivity than selegiline. Coumarin 7 is the most active compound of this series, being 64 times more active than selegiline and also showing the highest hMAO-B specificity. In addition, docking experiments were carried out on hMAO-A and h-MAO-B structures. This study provided new information about the enzyme–inhibitor interaction and the potential therapeutic application of this 3-arylcoumarin scaffold.
•We provide a critical review of the current state of the activity cliff research.•We comment on and integrate the opinions of experts working on activity cliffs.•The negative effects of activity ...cliffs on prediction models are discussed.•We provide a machine learning rationale for activity cliffs in chemoinformatics.•Potential solutions to address the negative effects of activity cliffs are proposed.
The impact activity cliffs have on drug discovery is double-edged. For instance, whereas medicinal chemists can take advantage of regions in chemical space rich in activity cliffs, QSAR practitioners need to escape from such regions. The influence of activity cliffs in medicinal chemistry applications is extensively documented. However, the ‘dark side’ of activity cliffs (i.e. their detrimental effect on the development of predictive machine learning algorithms) has been understudied. Similarly, limited amounts of work have been devoted to propose potential solutions to the drawbacks of activity cliffs in similarity-based approaches. In this review, the duality of activity cliffs in medicinal chemistry and computational approaches is addressed, with emphasis on the rationale and potential solutions for handling the ‘ugly face’ of activity cliffs.
Point of view of the current state of the activity cliff phenomenon focusing on the rationale, effects and potential solutions to handle the influence of activity cliffs in drug discovery.
Wuhan, China was the epicenter of the first zoonotic transmission of the severe acute respiratory syndrome coronavirus clade 2 (SARS-CoV-2) in December 2019 and it is the causative agent of the novel ...human coronavirus disease 2019 (COVID-19). Almost from the beginning of the COVID-19 outbreak several attempts were made to predict possible drugs capable of inhibiting the virus replication. In the present work a drug repurposing study is performed to identify potential SARS-CoV-2 protease inhibitors. We created a Quantitative Structure–Activity Relationship (QSAR) model based on a machine learning strategy using hundreds of inhibitor molecules of the main protease (Mpro) of the SARS-CoV coronavirus. The QSAR model was used for virtual screening of a large list of drugs from the DrugBank database. The best 20 candidates were then evaluated in-silico against the Mpro of SARS-CoV-2 by using docking and molecular dynamics analyses. Docking was done by using the Gold software, and the free energies of binding were predicted with the MM-PBSA method as implemented in AMBER. Our results indicate that levothyroxine, amobarbital and ABP-700 are the best potential inhibitors of the SARS-CoV-2 virus through their binding to the Mpro enzyme. Five other compounds showed also a negative but small free energy of binding: nikethamide, nifurtimox, rebimastat, apomine and rebastinib.
•Up to 51 different compounds were identified in the horchata infusion.•The infusion is predicted to have anti-inflammatory and antioxidant effects.•Horchata contains potential inhibitors of the ...SARS-CoV-2 Mpro and PLpro enzymes.
Bioactive plant-derived molecules have emerged as therapeutic alternatives in the fight against the COVID-19 pandemic. In this investigation, principal bioactive compounds of the herbal infusion “horchata” from Ecuador were studied as potential novel inhibitors of the SARS-CoV-2 virus. The chemical composition of horchata was determined through a HPLC-DAD/ESI-MSn and GC–MS analysis while the inhibitory potential of the compounds on SARS-CoV-2 was determined by a computational prediction using various strategies, such as molecular docking and molecular dynamics simulations. Up to 51 different compounds were identified. The computational analysis of predicted targets reveals the compounds’ possible anti-inflammatory (no steroidal) and antioxidant effects. Three compounds were identified as candidates for Mpro inhibition: benzoic acid, 2-(ethylthio)-ethyl ester, l-Leucine-N-isobutoxycarbonyl-N-methyl-heptyl and isorhamnetin and for PLpro: isorhamnetin-3-O-(6-Orhamnosyl-galactoside), dihydroxy-methoxyflavanone and dihydroxyphenyl)-5-hydroxy-4-oxochromen-7-yloxy-3,4,5-trihydroxyoxane-2-carboxylic acid. Our results suggest the potential of Ecuadorian horchata infusion as a starting scaffold for the development of new inhibitors of the SARS-CoV-2 Mpro and PLpro enzymes.
The mosquito
transmits the virus that causes dengue, yellow fever, Zika and Chikungunya viruses, and in several regions of the planet represents a vector of great clinical importance. In terms of ...mortality and morbidity, infections caused by
are among the most serious arthropod transmitted viral diseases. The present study investigated the larvicidal potential of seventeen cinnamic acid derivatives against fourth stage
larvae. The larvicide assays were performed using larval mortality rates to determine lethal concentration (LC
). Compounds containing the medium alkyl chains butyl cinnamate (
) and pentyl cinnamate (
) presented excellent larvicidal activity with LC
values of around 0.21-0.17 mM, respectively. While among the derivatives with aryl substituents, the best LC
result was 0.55 mM for benzyl cinnamate (
). The tested derivatives were natural compounds and in pharmacology and antiparasitic studies, many have been evaluated using biological models for environmental and toxicological safety. Molecular modeling analyses suggest that the larvicidal activity of these compounds might be due to a multi-target mechanism of action involving inhibition of a carbonic anhydrase (CA), a histone deacetylase (HDAC2), and two sodium-dependent cation-chloride co-transporters (CCC2 e CCC3).
Preeclampsia is a hypertensive disorder that occurs during pregnancy. It is a complex disease with unknown pathogenesis and the leading cause of fetal and maternal mortality during pregnancy. Using ...all drugs currently under clinical trial for preeclampsia, we extracted all their possible targets from the DrugBank and ChEMBL databases and labeled them as "targets". The proteins labeled as "off-targets" were extracted in the same way but while taking all antihypertensive drugs which are inhibitors of ACE and/or angiotensin receptor antagonist as query molecules. Classification models were obtained for each of the 55 total proteins (45 targets and 10 off-targets) using the TPOT pipeline optimization tool. The average accuracy of the models in predicting the external dataset for targets and off-targets was 0.830 and 0.850, respectively. The combinations of models maximizing their virtual screening performance were explored by combining the desirability function and genetic algorithms. The virtual screening performance metrics for the best model were: the Boltzmann-Enhanced Discrimination of ROC (BEDROC)
= 0.258, the Enrichment Factor (EF)
= 31.55 and the Area Under the Accumulation Curve (AUAC) = 0.831. The most relevant targets for preeclampsia were: AR, VDR, SLC6A2, NOS3 and CHRM4, while ABCG2, ERBB2, CES1 and REN led to the most relevant off-targets. A virtual screening of the DrugBank database identified estradiol, estriol, vitamins E and D, lynestrenol, mifrepristone, simvastatin, ambroxol, and some antibiotics and antiparasitics as drugs with potential application in the treatment of preeclampsia.
Breast cancer (BC) is the leading cause of cancer-related death among women and the most commonly diagnosed cancer worldwide. Although in recent years large-scale efforts have focused on identifying ...new therapeutic targets, a better understanding of BC molecular processes is required. Here we focused on elucidating the molecular hallmarks of BC heterogeneity and the oncogenic mutations involved in precision medicine that remains poorly defined. To fill this gap, we established an OncoOmics strategy that consists of analyzing genomic alterations, signaling pathways, protein-protein interactome network, protein expression, dependency maps in cell lines and patient-derived xenografts in 230 previously prioritized genes to reveal essential genes in breast cancer. As results, the OncoOmics BC essential genes were rationally filtered to 140. mRNA up-regulation was the most prevalent genomic alteration. The most altered signaling pathways were associated with basal-like and Her2-enriched molecular subtypes. RAC1, AKT1, CCND1, PIK3CA, ERBB2, CDH1, MAPK14, TP53, MAPK1, SRC, RAC3, BCL2, CTNNB1, EGFR, CDK2, GRB2, MED1 and GATA3 were essential genes in at least three OncoOmics approaches. Drugs with the highest amount of clinical trials in phases 3 and 4 were paclitaxel, docetaxel, trastuzumab, tamoxifen and doxorubicin. Lastly, we collected ~3,500 somatic and germline oncogenic variants associated with 50 essential genes, which in turn had therapeutic connectivity with 73 drugs. In conclusion, the OncoOmics strategy reveals essential genes capable of accelerating the development of targeted therapies for precision oncology.