Search

The COVID-19 pandemic continues to pose a substantial threat to human lives and is likely to do so for years to come. Despite the availability of vaccines, searching for efficient small-molecule drugs that are widely available, including in low- and middle-income countries, is an ongoing challenge. In this work, we report the results of a community effort, the “Billion molecules against Covid-19 challenge”, to identify small-molecule inhibitors against SARS-CoV-2 or relevant human receptors. Participating teams used a wide variety of computational methods to screen a minimum of 1 billion virtual molecules against 6 protein targets. Overall, 31 teams participated, and they suggested a total of 639,024 potentially active molecules, which were subsequently ranked to find ‘consensus compounds’. The organizing team coordinated with various contract research organizations (CROs) and collaborating institutions to synthesize and test 878 compounds for activity against proteases (Nsp5, Nsp3, TMPRSS2), nucleocapsid N, RdRP (Nsp12 domain), and (alpha) spike protein S. Overall, 27 potential inhibitors were experimentally confirmed by binding-, cleavage-, and/or viral suppression assays and are presented here. All results are freely available and can be taken further downstream without IP restrictions. Overall, we show the effectiveness of computational techniques, community efforts, and communication across research fields (i.e., protein expression and crystallography, in silico modeling, synthesis and biological assays) to accelerate the early phases of drug discovery.

Nuclear Receptors (NRs) are highly relevant drug targets, for which small molecule modulation goes beyond a simple ligand/receptor interaction. NR–ligands modulate Protein–Protein Interactions (PPIs) with coregulator proteins. Here we bring forward a cooperativity mechanism for small molecule modulation of NR PPIs, using the Peroxisome Proliferator Activated Receptor γ (PPARγ), which describes NR–ligands as allosteric molecular glues. The cooperativity framework uses a thermodynamic model based on three-body binding events, to dissect and quantify reciprocal effects of NR–coregulator binding (KID) and NR–ligand binding (KIID), jointly recapitulated in the cooperativity factor (α) for each specific ternary ligand·NR·coregulator complex formation. These fundamental thermodynamic parameters allow for a conceptually new way of thinking about structure–activity-relationships for NR–ligands and can steer NR modulator discovery and optimization via a completely novel approach.

The EPIC consortium brings together experts from a wide range of fields that include clinical, molecular and basic microbiology, infectious diseases, computational biology and chemistry, drug discovery and design, bioinformatics, biochemistry, biophysics, pharmacology, toxicology, veterinary sciences, environmental sciences, and epidemiology. The main question to be answered by the EPIC alliance is the following: “What is the best approach for data mining on carbapenemase inhibitors and how to translate this data into experiments?” From this forum, we propose that the scientific community think up new strategies to be followed for the discovery of new carbapenemase inhibitors, so that this process is efficient and capable of providing results in the shortest possible time and within acceptable time and economic costs.

Cholesterol biosynthesis is a multistep process in mammals that includes the aerobic removal of three methyl groups from the intermediate lanosterol, one from position 14 and two from position 4. During the demethylations at position 4, a 3-ketosteroid reductase catalyses the conversion of both 4-methylzymosterone and zymosterone to 4-methylzymosterol and zymosterol, respectively, restoring the alcoholic function of lanosterol, which is also maintained in cholesterol. Unlike other eukaryotes, mammals also use the same enzyme as an estrone reductase that can transform estrone (E1) into estradiol (E2). This enzyme, named 17β-hydroxysteroid dehydrogenase type 7 (HSD17B7), is therefore a multifunctional protein in mammals, and one that belongs to both the HSD17B family, which is involved in steroid-hormone metabolism, and to the family of post-squalene cholesterol biosynthesis enzymes. In the present study, a series of known inhibitors of human HSD17B7's E1-reductase activity have been assayed for potential inhibition against 3-ketosteroid reductase activity. Surprisingly, the assayed compounds lost their inhibition activity when tested in HepG2 cells that were incubated with radiolabelled acetate and against the recombinant overexpressed human enzyme incubated with 4-methylzymosterone (both radiolabelled and not). Preliminary kinetic analyses suggest a mixed or non-competitive inhibition on the E1-reductase activity, which is in agreement with Molecular Dynamics simulations. These results raise questions about the mechanism(s) of action of these possible inhibitors, the enzyme dynamic regulation and the interplay between the two activities.

Iron surface determinant B (IsdB) is a hemoglobin (Hb) receptor essential for hemic iron acquisition by Staphylococcus aureus. Heme transfer to IsdB is possible from oxidized Hb (metHb), but inefficient from Hb either bound to oxygen (oxyHb) or bound to carbon monoxide (HbCO), and encompasses a sequence of structural events that are currently poorly understood. By single-particle cryo-electron microscopy, we determined the structure of two IsdB:Hb complexes, representing key species along the heme extraction pathway. The IsdB:HbCO structure, at 2.9-Å resolution, provides a snapshot of the preextraction complex. In this early stage of IsdB:Hb interaction, the hemophore binds to the β-subunits of the Hb tetramer, exploiting a folding-upon-binding mechanism that is likely triggered by a cis/trans isomerization of Pro173. Binding of IsdB to α-subunits occurs upon dissociation of the Hb tetramer into α/β dimers. The structure of the IsdB:metHb complex reveals the final step of the extraction process, where heme transfer to IsdB is completed. The stability of the complex, both before and after heme transfer from Hb to IsdB, is influenced by isomerization of Pro173. These results greatly enhance current understanding of structural and dynamic aspects of the heme extraction mechanism by IsdB and provide insight into the interactions that stabilize the complex before the heme transfer event. This information will support future efforts to identify inhibitors of heme acquisition by S. aureus by interfering with IsdB:Hb complex formation.

Nuclear receptors (NRs) are transcription factors that play an important role in multiple diseases, such as cancer, inflammation, and metabolic disorders. They share a common structural organization composed of five domains, of which the ligand-binding domain (LBD) can adopt different conformations in response to substrate, agonist, and antagonist binding, leading to distinct transcription effects. A key feature of NRs is, indeed, their intrinsic dynamics that make them a challenging target in drug discovery. This work aims to provide a meaningful investigation of NR structural variability to outline a dynamic profile for each of them. To do that, we propose a methodology based on the computation and comparison of protein cavities among the crystallographic structures of NR LBDs. First, pockets were detected with the FLAPsite algorithm and then an “all against all” approach was applied by comparing each pair of pockets within the same sub-family on the basis of their similarity score. The analysis concerned all the detectable cavities in NRs, with particular attention paid to the active site pockets. This approach can guide the investigation of NR intrinsic dynamics, the selection of reference structures to be used in drug design and the easy identification of alternative binding sites.

Three open-source anti-kinetoplastid chemical boxes derived from a whole-cell phenotypic screening by GlaxoSmithKline (Tres Cantos Anti-Kinetoplastid Screening, TCAKS) were exploited for the discovery of a novel core structure inspiring new treatments of parasitic diseases targeting the trypansosmatidic pteridine reductase 1 (PTR1) and dihydrofolate reductase (DHFR) enzymes. In total, 592 compounds were tested through medium-throughput screening assays. A subset of 14 compounds successfully inhibited the enzyme activity in the low micromolar range of at least one of the enzymes from both Trypanosoma brucei and Lesihmania major parasites (pan-inhibitors), or from both PTR1 and DHFR-TS of the same parasite (dual inhibitors). Molecular docking studies of the protein–ligand interaction focused on new scaffolds not reproducing the well-known antifolate core clearly explaining the experimental data. TCMDC-143249, classified as a benzenesulfonamide derivative by the QikProp descriptor tool, showed selective inhibition of PTR1 and growth inhibition of the kinetoplastid parasites in the 5 μM range. In our work, we enlarged the biological profile of the GSK Kinetobox and identified new core structures inhibiting selectively PTR1, effective against the kinetoplastid infectious protozoans. In perspective, we foresee the development of selective PTR1 and DHFR inhibitors for studies of drug combinations.

In the search for new chemical scaffolds able to afford NLRP3 inflammasome inhibitors, we used a pharmacophore-hybridization strategy by combining the structure of the acrylic acid derivative INF39 with the 1-(piperidin-4-yl)1,3-dihydro-2H-benzo[d]imidazole-2-one substructure present in HS203873, a recently identified NLRP3 binder. A series of differently modulated benzo[d]imidazole-2-one derivatives were designed and synthesised. The obtained compounds were screened in vitro to test their ability to inhibit NLRP3-dependent pyroptosis and IL-1β release in PMA-differentiated THP-1 cells stimulated with LPS/ATP. The selected compounds were evaluated for their ability to reduce the ATPase activity of human recombinant NLRP3 using a newly developed assay. From this screening, compounds 9, 13 and 18, able to concentration-dependently inhibit IL-1β release in LPS/ATP-stimulated human macrophages, emerged as the most promising NLRP3 inhibitors of the series. Computational simulations were applied for building the first complete model of the NLRP3 inactive state and for identifying possible binding sites available to the tested compounds. The analyses led us to suggest a mechanism of protein–ligand binding that might explain the activity of the compounds.

The generation of reactive oxygen species (ROS) and reactive nitrogen species (RNS) as “unconventional” therapeutics with precise spatiotemporal control by using light stimuli may open entirely new horizons for innovative therapeutic modalities. Among ROS and RNS, peroxynitrite (ONOO−) plays a dominant role in chemistry and biology in view of its potent oxidizing power and cytotoxic action. We have designed and synthesized a molecular hybrid based on benzophenothiazine as a red light-harvesting antenna joined to an N-nitroso appendage through a flexible spacer. Single photon red light excitation of this molecular construct triggers the release of nitric oxide (˙NO) and simultaneously produces superoxide anions (O2˙−). The diffusion-controlled reaction between these two radical species generates ONOO−, as confirmed by the use of fluorescein-boronate as a highly selective chemical probe. Besides, the red fluorescence of the hybrid allows its tracking in different types of cancer cells where it is well-tolerated in the dark but induces remarkable cell mortality under irradiation with red light in a very low concentration range, with very low light doses (ca. 1 J cm−2). This ONOO− generator activatable by highly biocompatible and tissue penetrating single photon red light can open up intriguing prospects in biomedical research, where precise and spatiotemporally controlled concentrations of ONOO− are required., Excitation of a molecular hybrid with highly biocompatible red light generates cytotoxic peroxynitrite, produces red fluorescence useful for cell tracking and induces remarkable cancer cell death at very low concentrations and very low light doses.

The aldo-keto reductase 1C3 (AKR1C3) enzyme is considered an attractive target in Castration Resistant Prostate Cancer (CRPC) because of its role in the biosynthesis of androgens. Flufenamic acid, a non-selective AKR1C3 inhibitor, has previously been subjected to bioisosteric modulation to give rise to a series of compounds with the hydroxytriazole core. In this work, the hit compound of the previous series has been modulated further, and new, more potent, and selective derivatives have been obtained. The poor solubility of the most active compound (cpd 5) has been improved by substituting the triazole core with an isoxazole heteronucleous, with similar enzymatic activity being retained. Potent AKR1C3 inhibition is translated into antiproliferative effects against the 22RV1 CRPC cellular model, and the in-silico design, synthesis and biological activity of new compounds are described herein. Compounds have also been assayed in combination with two approved antitumor drugs, abiraterone and enzalutamide.

C6O4 is a water soluble perfluoroether carboxylic acid ammonium salt used as surfactant in the synthesis of fluoropolymers. Available experimental data in rats exposed by the oral route indicate it is eliminated in urine. Previous studies with various linear perfluorocarboxylic acids have suggested that these compounds are substrates of renal membrane transporters in rats and humans, and that the interaction with basal and apical membrane transporters can influence the elimination kinetic by these organisms and explain, in part, the observed differences in the respective half-lives. In particular, apical transporters may contribute to the reuptake of these exogenous compounds from the tubule lumen. The present study was designed to investigate the uptake of C6O4 in two renal cell lines transiently transfected with the human apical membrane transporters, organic anion transporter 4 (OAT4), and urate transporter 1 (URAT1). The uptake of the linear perfluorohexanoic acid (PFC6) was evaluated in parallel. While the uptake of the conjugated steroid estrone-3-sulfate (E3S), a known substrate for renal transporters, and of PFC6 was clearly observed in both cell types transfected with either OAT4 or URAT1, no significant uptake of C6O4 was measured under the same test conditions. The results of the transporter's functionality measured in vitro were consistent with molecular docking simulations. Both outward and inward models of the transporters showed a reduced interaction between C6O4 and URAT1 or OAT4. In contrast, more stable interactions were predicted for PFC6 and PFOA, as well as for the E3S substrate, as shown by the respective docking scores reflecting the binding strength and by the poses assumed in the transporter channels. Altogether, the in vitro and in silico modeling results showed a low reuptake potential and limited interactions of C6O4 molecule with two human apical membrane transporters, contrasting with the more efficient reuptake of PFC6 from the tubule lumen. These results suggest reabsorption from the proximal tubule by apical renal transporters is not likely to interfere with the elimination pathway of C6O4 in humans.

Drug-target interaction, cellular internalization, and target engagement should be addressed to design a lead with high chances of success in further optimization stages. Accordingly, we have designed conjugates of folic acid with anticancer peptides able to bind human thymidylate synthase (hTS) and enter cancer cells through folate receptor α (FRα) highly expressed by several cancer cells. Mechanistic analyses and molecular modeling simulations have shown that these conjugates bind the hTS monomer-monomer interface with affinities over 20 times larger than the enzyme active site. When tested on several cancer cell models, these conjugates exhibited FRα selectivity at nanomolar concentrations. A similar selectivity was observed when the conjugates were delivered in synergistic or additive combinations with anticancer agents. At variance with 5-fluorouracil and other anticancer drugs that target the hTS catalytic pocket, these conjugates do not induce overexpression of this protein and can thus help combating drug resistance associated with high hTS levels.

Molecules online 26(6), 1613 (2021). doi:10.3390/molecules26061613 special issue: "Special Issue "Computational Methods in Drug Design and Food Chemistry" / Special Issue Editor: Prof. Dr. Giosuè Costa, Guest Editor", Published by Springer, Berlin ; Heidelberg

Strigolactones (SLs) are a class of sesquiterpenoid plant hormones that play a role in the response of plants to various biotic and abiotic stresses. When released into the rhizosphere, they are perceived by both beneficial symbiotic mycorrhizal fungi and parasitic plants. Due to their multiple roles, SLs are potentially interesting agricultural targets. Indeed, the use of SLs as agrochemicals can favor sustainable agriculture via multiple mechanisms, including shaping root architecture, promoting ideal branching, stimulating nutrient assimilation, controlling parasitic weeds, mitigating drought and enhancing mycorrhization. Moreover, over the last few years, a number of studies have shed light onto the effects exerted by SLs on human cells and on their possible applications in medicine. For example, SLs have been demonstrated to play a key role in the control of pathways related to apoptosis and inflammation. The elucidation of the molecular mechanisms behind their action has inspired further investigations into their effects on human cells and their possible uses as anti-cancer and antimicrobial agents.

The SARS-CoV-2 coronavirus outbreak continues to spread at a rapid rate worldwide. The main protease (Mpro) is an attractive target for anti-COVID-19 agents. Unfortunately, unexpected difficulties have been encountered in the design of specific inhibitors. Here, by analyzing an ensemble of ~30,000 SARS-CoV-2 Mpro conformations from crystallographic studies and molecular simulations, we show that small structural variations in the binding site dramatically impact ligand binding properties. Hence, traditional druggability indices fail to adequately discriminate between highly and poorly druggable conformations of the binding site. By performing ~200 virtual screenings of compound libraries on selected protein structures, we redefine the protein’s druggability as the consensus chemical space arising from the multiple conformations of the binding site formed upon ligand binding. This procedure revealed a unique SARS-CoV-2 Mpro blueprint that led to a definition of a specific structure-based pharmacophore. The latter explains the poor transferability of potent SARS-CoV Mpro inhibitors to SARS-CoV-2 Mpro, despite the identical sequences of the active sites. Importantly, application of the pharmacophore predicted novel high affinity inhibitors of SARS-CoV-2 Mpro, that were validated by in vitro assays performed here and by a newly solved X-ray crystal structure. These results provide a strong basis for effective rational drug design campaigns against SARS-CoV-2 Mpro and a new computational approach to screen protein targets with malleable binding sites.

A novel molecular hybrid has been designed and synthesized in which acridine orange (AO) is covalently linked to an N-nitrosoaniline derivative through an alkyl spacer. Photoexcitation of the AO antenna with the highly biocompatible green light results in intense fluorescence emission and triggers NO detachment from the N-nitroso appendage via an intramolecular electron transfer. The presence of the AO moiety encourages the binding with DNA through both external and partially intercalative fashions, depending on the DNA:molecular hybrid molar ratio. Importantly, this dual-mode binding interaction with the biopolymer does not preclude the NO photoreleasing performances of the molecular hybrid, permitting NO to be photogenerated nearby DNA with an efficiency similar to that of the free molecule. These properties make the presented compound an intriguing candidate for fundamental and potential applicative research studies where NO delivery in the DNA proximity precisely regulated by harmless green light is required.

The emergence of bacteria that co-express serine- and metallo- carbapenemases is a threat to the efficacy of the available &beta, lactam antibiotic armamentarium. The 4-amino-1,2,4-triazole-3-thione scaffold has been selected as the starting chemical moiety in the design of a small library of &beta, Lactamase inhibitors (BLIs) with extended activity profiles. The synthesised compounds have been validated in vitro against class A serine &beta, &minus, Lactamase (SBLs) KPC-2 and class B1 metallo &beta, Lactamases (MBLs) VIM-1 and IMP-1. Of the synthesised derivatives, four compounds showed cross-class micromolar inhibition potency and therefore underwent in silico analyses to elucidate their binding mode within the catalytic pockets of serine- and metallo-BLs. Moreover, several members of the synthesised library have been evaluated, in combination with meropenem (MEM), against clinical strains that overexpress BLs for their ability to synergise carbapenems.

We report for the first time a NO photodonor (NOPD) operating with the widely used chemotherapeutic agent doxorubicin (DOX) as the light-harvesting antenna. This permits NO uncaging from an N-nitroso appendage upon selective excitation of DOX with highly biocompatible green light, without precluding its typical red emission. This NOPD effectively binds DNA and photodelivers NO nearby, representing an intriguing candidate for potential multimodal therapeutic applications based on the combination of DOX and NO.

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

International audience; Despite prophylactic vaccination campaigns, human papillomavirus (HPV)-induced cancers still represent a major medical issue for global population, thus specific anti-HPV drugs are needed. Since the ability of HPV E6 oncoprotein to promote p53 degradation is linked to tumor progression, E6 has been proposed as an ideal target for cancer treatment. Using the crystal structure of the E6/E6AP/p53 complex, we performed an in silico screening of small-molecule libraries against a highly conserved alpha-helix in the N-terminal domain of E6 involved in the E6-p53 interaction. We discovered a compound able to inhibit the E6-mediated degradation of p53 through disruption of E6-p53 binding both in vitro and in cells. This compound could restore p53 intracellular levels and transcriptional activity, reduce the viability and proliferation of HPV-positive cancer cells, and block 3D cervospheres formation. Mechanistic studies revealed that the compound anti-tumor activity mainly relies on induction of cell cycle arrest and senescence. Our data demonstrate that the disruption of the direct E6-p53 interaction can be obtained with a small-molecule compound leading to specific antitumoral activity in HPV-positive cancer cells and thus represents a new approach for anti-HPV drug development.

[Image: see text] The value of including protein flexibility in structure-based drug design (SBDD) is widely documented, and currently, molecular dynamics (MD) simulations represent a powerful tool to investigate protein dynamics. Yet, the inclusion of MD-derived information in pre-existing SBDD workflows is still far from trivial. We recently published an integrated MD-FLAP (Fingerprints for Ligands and Proteins) approach combining MD, clustering and Linear Discriminant Analysis (LDA) for enhancing accuracy, efficacy, and for protein conformational selection in virtual screening (VS) campaigns. Here we prospectively applied the MD-FLAP workflow to discover novel chemotypes inhibiting the Casein Kinase 1 delta (CSNK1D) enzyme. We first obtained a VS model able to separate active from inactive compounds, with a global AUC of 0.9 and a partial ROC enrichment at 0.5% of 0.18, and use it to mine the internal Pfizer screening database. Seven active molecules sharing a phenyl-indazole scaffold, not yet reported among CSNK1D inhibitors, were found. The most potent inhibitor showed an IC(50) of 134 nM.

Paracetamol has been one of the most commonly used and prescribed analgesic drugs for more than a hundred years. Despite being generally well tolerated, it can result in high liver toxicity when administered in specific conditions, such as overdose, or in vulnerable individuals. We have synthesized and characterized a paracetamol galactosylated prodrug (PARgal) with the aim of improving both the pharmacodynamic and pharmacological profile of paracetamol. PARgal shows a range of physicochemical properties, solubility, lipophilicity, and chemical stability at differing physiological pH values and in human serum. PARgal could still be preclinically detected 2 h after administration, meaning that it displays reduced hepatic metabolism compared to paracetamol. In overdose conditions, PARgal has not shown any cytotoxic effect in in vitro analyses performed on human liver cells. Furthermore, when tested in an animal pain model, PARgal demonstrated a sustained analgesic effect up to the 12th hour after oral administration. These findings support the use of galactose as a suitable carrier in the development of prodrugs for analgesic treatment.

Due to the poor aqueous solubility of retinoids, evolution has tuned their binding to cellular proteins to address specialized physiological roles by modulating uptake, storage, and delivery to specific targets. With the aim to disentangle the structure-function relationships in these proteins and disclose clues for engineering selective carriers, the binding mechanism of the two most abundant retinol-binding isoforms was explored by using enhanced sampling molecular dynamics simulations and surface plasmon resonance. The distinctive dynamics of the entry portal site in the holo species was crucial to modulate retinol dissociation. Remarkably, this process is controlled at large extent by the replacement of Ile by Leu in the two isoforms, thus suggesting that a fine control of ligand release can be achieved through a rigorous selection of conservative mutations in accessory sites.

The worldwide emergence of New Delhi metallo-β-lactamase-1 (NDM-1) as a carbapenemase able to hydrolyze nearly all available β-lactam antibiotics has characterized the past decade, endangering efficacious antibacterial treatments. No inhibitors for NDM-1 are available in therapy, nor are promising compounds in the pipeline for future NDM-1 inhibitors. We report the studies dedicated to the design and development of effective NDM-1 inhibitors. The discussion for each agent moves from the employed design strategy to the ability of the identified inhibitor to synergize β-lactam antibiotics. A structural analysis of NDM-1 mechanism of action based on selected X-ray complexes is also reported: the intrinsic flexibility of the binding site and the comparison between penicillin/cephalosporin and carbapenem mechanisms of hydrolysis are evaluated. Despite the valuable progress in terms of structural and mechanistic information, the design of a potent NDM-1 inhibitor to be introduced in therapy remains challenging. Certainly, only the deep knowledge of NDM-1 architecture and of the variable mechanism of action that NDM-1 employs against different classes of substrates could orient a successful drug discovery campaign.

Two novel NO photodonors (NOPDs) based on BODIPY and Rhodamine antennae activatable with the highly biocompatible green light are reported. Both NOPDs exhibit considerable fluorescence emission and release NO with remarkable quantum efficiencies. The combination of the photoreleasing and emissive performance for both compounds is superior to those exhibited by other NOPDs based on similar light-harvesting centres, making them very intriguing for image-guided phototherapeutic applications. Preliminary biological data prove their easy visualization in cell environment due to the intense green and orange-red fluorescence and their photodynamic action on cancer cells due to the NO photo-liberated.