The TBC1D31/praja2 complex controls primary ciliogenesis through PKA-directed OFD1 ubiquitylation.

The primary cilium is a microtubule-based sensory organelle that dynamically links signalling pathways to cell differentiation, growth, and development. Genetic defects of primary cilia are responsible for genetic disorders known as ciliopathies. Orofacial digital type I syndrome (OFDI) is an X-linked congenital ciliopathy caused by mutations in the OFD1 gene and characterized by malformations of the face, oral cavity, digits and, in the majority of cases, polycystic kidney disease. OFD1 plays a key role in cilium biogenesis. However, the impact of signalling pathways and the role of the ubiquitin-proteasome system (UPS) in the control of OFD1 stability remain unknown. Here, we identify a novel complex assembled at centrosomes by TBC1D31, including the E3 ubiquitin ligase praja2, protein kinase A (PKA), and OFD1. We show that TBC1D31 is essential for ciliogenesis. Mechanistically, upon G-protein-coupled receptor (GPCR)-cAMP stimulation, PKA phosphorylates OFD1 at ser735, thus promoting OFD1 proteolysis through the praja2-UPS circuitry. This pathway is essential for ciliogenesis. In addition, a non-phosphorylatable OFD1 mutant dramatically affects cilium morphology and dynamics. Consistent with a role of the TBC1D31/praja2/OFD1 axis in ciliogenesis, alteration of this molecular network impairs ciliogenesis in vivo in Medaka fish, resulting in developmental defects. Our findings reveal a multifunctional transduction unit at the centrosome that links GPCR signalling to ubiquitylation and proteolysis of the ciliopathy protein OFD1, with important implications on cilium biology and development. Derangement of this control mechanism may underpin human genetic disorders.

Ubiquitylation of BBSome is required for ciliary assembly and signaling.

Bardet-Biedl syndrome (BBS) is a ciliopathy characterized by retinal degeneration, obesity, renal abnormalities, postaxial polydactyly, and developmental defects. Genes mutated in BBS encode for components and regulators of the BBSome, an octameric complex that controls the trafficking of cargos and receptors within the primary cilium. Although both structure and function of the BBSome have been extensively studied, the impact of ubiquitin signaling on BBSome is largely unknown. We identify the E3 ubiquitin ligase PJA2 as a novel resident of the ciliary compartment and regulator of the BBSome. Upon GPCR-cAMP stimulation, PJA2 ubiquitylates BBSome subunits. We demonstrate that ubiquitylation of BBS1 at lysine 143 increases the stability of the BBSome and promotes its binding to BBS3, an Arf-like GTPase protein controlling the targeting of the BBSome to the ciliary membrane. Downregulation of PJA2 or expression of a ubiquitylation-defective BBS1 mutant (BBS1K143R ) affects the trafficking of G-protein-coupled receptors (GPCRs) and Shh-dependent gene transcription. Expression of BBS1K143R in vivo impairs cilium formation, embryonic development, and photoreceptors' morphogenesis, thus recapitulating the BBS phenotype in the medaka fish model.

Discovery of Bile Acid Derivatives as Potent ACE2 Activators by Virtual Screening and Essential Dynamics.

The angiotensin-converting enzyme II (ACE2) is a key molecular player in the regulation of vessel contraction, inflammation, and reduction of oxidative stress. In addition, ACE2 has assumed a prominent role in the fight against the COVID-19 pandemic-causing virus SARS-CoV-2, as it is the very first receptor in the host of the viral spike protein. The binding of the spike protein to ACE2 triggers a cascade of events that eventually leads the virus to enter the host cell and initiate its life cycle. At the same time, SARS-CoV-2 infection downregulates ACE2 expression especially in the lung, altering the biochemical signals regulated by the enzyme and contributing to the poor clinical prognosis characterizing the late stage of the COVID-19 disease. Despite its important biological role, a very limited number of ACE2 activators are known. Here, using a combined in silico and experimental approach, we show that ursodeoxycholic acid (UDCA) derivatives work as ACE2 activators. In detail, we have identified two potent ACE2 ligands, BAR107 and BAR708, through a docking virtual screening campaign and elucidated their mechanism of action from essential dynamics of the enzyme observed during microsecond molecular dynamics calculations. The in silico results were confirmed by in vitro pharmacological assays with the newly identified compounds showing ACE2 activity comparable to that of DIZE, the most potent ACE2 activator known so far. Our work provides structural insight into ACE2/ligand-binding interaction useful for the design of compounds with therapeutic potential against SARS-CoV-2 infection, inflammation, and other ACE2-related diseases.

Molecular Basis for Non-Covalent, Non-Competitive FAAH Inhibition.

Fatty acid amide hydrolase (FAAH) plays a key role in the control of cannabinoid signaling and it represents a promising therapeutic strategy for the treatment of a wide range of diseases, including neuropathic pain and chronic inflammation. Starting from kinetics experiments carried out in our previous work for the most potent inhibitor 2-amino-3-chloropyridine amide (TPA14), we have investigated its non-competitive mechanism of action using molecular dynamics, thermodynamic integration and QM-MM/GBSA calculations. The computational studies highlighted the impact of mutations on the receptor binding pockets and elucidated the molecular basis of the non-competitive inhibition mechanism of TPA14, which prevents the endocannabinoid anandamide (AEA) from reaching its pro-active conformation. Our study provides a rationale for the design of non-competitive potent FAAH inhibitors for the treatment of neuropathic pain and chronic inflammation.

Design, synthesis and in vitro and in vivo biological evaluation of flurbiprofen amides as new fatty acid amide hydrolase/cyclooxygenase-2 dual inhibitory potential analgesic agents.

Compounds combining dual inhibitory action against FAAH and cyclooxygenase (COX) may be potentially useful analgesics. Here, we describe a novel flurbiprofen analogue, N-(3-bromopyridin-2-yl)-2-(2-fluoro-(1,1'-biphenyl)-4-yl)propanamide (Flu-AM4). The compound is a competitive, reversible inhibitor of FAAH with a Ki value of 13 nM and which inhibits COX activity in a substrate-selective manner. Molecular modelling suggested that Flu-AM4 optimally fits a hydrophobic pocket in the ACB region of FAAH, and binds to COX-2 similarly to flurbiprofen. In vivo studies indicated that at a dose of 10 mg/kg, Flu-AM4 was active in models of prolonged (formalin) and neuropathic (chronic constriction injury) pain and reduced the spinal expression of iNOS, COX-2, and NFκB in the neuropathic model. Thus, the present study identifies Flu-AM4 as a dual-action FAAH/substrate-selective COX inhibitor with anti-inflammatory and analgesic activity in animal pain models. These findings underscore the potential usefulness of such dual-action compounds.

Exploring the fatty acid amide hydrolase and cyclooxygenase inhibitory properties of novel amide derivatives of ibuprofen.

Inhibition of fatty acid amide hydrolase (FAAH) reduces the gastrointestinal damage produced by non-steroidal anti-inflammatory agents such as sulindac and indomethacin in experimental animals, suggesting that a dual-action FAAH-cyclooxygenase (COX) inhibitor could have useful therapeutic properties. Here, we have investigated 12 novel amide analogues of ibuprofen as potential dual-action FAAH/COX inhibitors. N-(3-Bromopyridin-2-yl)-2-(4-isobutylphenyl)propanamide (Ibu-AM68) was found to inhibit the hydrolysis of [3H]anandamide by rat brain homogenates by a reversible, mixed-type mechanism of inhibition with a Ki value of 0.26 µM and an α value of 4.9. At a concentration of 10 µM, the compound did not inhibit the cyclooxygenation of arachidonic acid by either ovine COX-1 or human recombinant COX-2. However, this concentration of Ibu-AM68 greatly reduced the ability of the COX-2 to catalyse the cyclooxygenation of the endocannabinoid 2-arachidonoylglycerol. It is concluded that Ibu-AM68 is a dual-acting FAAH/substrate-selective COX inhibitor.

Ligand binding to telomeric G-quadruplex DNA investigated by funnel-metadynamics simulations.

G-quadruplexes (G4s) are higher-order DNA structures typically present at promoter regions of genes and telomeres. Here, the G4 formation decreases the replicative DNA at each cell cycle, finally leading to apoptosis. The ability to control this mitotic clock, particularly in cancer cells, is fascinating and passes through a rational understanding of the ligand/G4 interaction. We demonstrate that an accurate description of the ligand/G4 binding mechanism is possible using an innovative free-energy method called funnel-metadynamics (FM), which we have recently developed to investigate ligand/protein interaction. Using FM simulations, we have elucidated the binding mechanism of the anticancer alkaloid berberine to the human telomeric G4 (d[AG3(T2AG3)3]), computing also the binding free-energy landscape. Two ligand binding modes have been identified as the lowest energy states. Furthermore, we have found prebinding sites, which are preparatory to reach the final binding mode. In our simulations, the ions and the water molecules have been explicitly represented and the energetic contribution of the solvent during ligand binding evaluated. Our theoretical results provide an accurate estimate of the absolute ligand/DNA binding free energy ([Formula: see text] = -10.3 ± 0.5 kcal/mol) that we validated through steady-state fluorescence binding assays. The good agreement between the theoretical and experimental value demonstrates that FM is a most powerful method to investigate ligand/DNA interaction and can be a useful tool for the rational design also of G4 ligands.

Multi-Targeting Bioactive Compounds Extracted from Essential Oils as Kinase Inhibitors.

Essential oils (EOs) are popular in aromatherapy, a branch of alternative medicine that claims their curative effects. Moreover, several studies reported EOs as potential anti-cancer agents by inducing apoptosis in different cancer cell models. In this study, we have considered EOs as a potential resource of new kinase inhibitors with a polypharmacological profile. On the other hand, computational methods offer the possibility to predict the theoretical activity profile of ligands, discovering dangerous off-targets and/or synergistic effects due to the potential multi-target action. With this aim, we performed a Structure-Based Virtual Screening (SBVS) against X-ray models of several protein kinases selected from the Protein Data Bank (PDB) by using a chemoinformatics database of EOs. By evaluating theoretical binding affinity, 13 molecules were detected among EOs as new potential kinase inhibitors with a multi-target profile. The two compounds with higher percentages in the EOs were studied more in depth by means Induced Fit Docking (IFD) protocol, in order to better predict their binding modes taking into account also structural changes in the receptor. Finally, given its good binding affinity towards five different kinases, cinnamyl cinnamate was biologically tested on different cell lines with the aim to verify the antiproliferative activity. Thus, this work represents a starting point for the optimization of the most promising EOs structure as kinase inhibitors with multi-target features.

Benzylamides and piperazinoarylamides of ibuprofen as fatty acid amide hydrolase inhibitors.

Fatty Acid Amide Hydrolase (FAAH) is a serine hydrolase that plays a key role in controlling endogenous levels of endocannabinoids. FAAH inhibition is considered a powerful approach to enhance the endocannabinoid signalling, and therefore it has been largely studied as a potential target for the treatment of neurological disorders such as anxiety or depression, or of inflammatory processes. We present two novel series of amide derivatives of ibuprofen designed as analogues of our reference FAAH inhibitor Ibu-AM5 to further explore its structure-activity relationships. In the new amides, the 2-methylpyridine moiety of Ibu-AM5 was substituted by benzylamino and piperazinoaryl moieties. The obtained benzylamides and piperazinoarylamides showed FAAH inhibition ranging from the low to high micromolar potency. The binding of the new amides in the active site of FAAH, estimated using the induced fit protocol, indicated arylpiperazinoamides binding the ACB channel and the cytosolic port, and benzylamides binding the ACB channel.

Identification of G-quadruplex DNA/RNA binders: Structure-based virtual screening and biophysical characterization.

BACKGROUND: Recent findings demonstrated that, in mammalian cells, telomere DNA (Tel) is transcribed into telomeric repeat-containing RNA (TERRA), which is involved in fundamental biological processes, thus representing a promising anticancer target. For this reason, the discovery of dual (as well as selective) Tel/TERRA G-quadruplex (G4) binders could represent an innovative strategy to enhance telomerase inhibition. METHODS: Initially, docking simulations of known Tel and TERRA active ligands were performed on the 3D coordinates of bimolecular G4 Tel DNA (Tel2) and TERRA (TERRA2). Structure-based pharmacophore models were generated on the best complexes and employed for the virtual screening of ~257,000 natural compounds. The 20 best candidates were submitted to biophysical assays, which included circular dichroism and mass spectrometry at different K+ concentrations. RESULTS: Three hits were here identified and characterized by biophysical assays. Compound 7 acts as dual Tel2/TERRA2 G4-ligand at physiological KCl concentration, while hits 15 and 17 show preferential thermal stabilization for Tel2 DNA. The different molecular recognition against the two targets was also discussed. CONCLUSIONS: Our successful results pave the way to further lead optimization to achieve both increased selectivity and stabilizing effect against TERRA and Tel DNA G4s. GENERAL SIGNIFICANCE: The current study combines for the first time molecular modelling and biophysical assays applied to bimolecular DNA and RNA G4s, leading to the identification of innovative ligand chemical scaffolds with a promising anticancer profile. This article is part of a Special Issue entitled "G-quadruplex" Guest Editor: Dr. Concetta Giancola and Dr. Daniela Montesarchio.

Chemoinformatic Database Building and in Silico Hit-Identification of Potential Multi-Targeting Bioactive Compounds Extracted from Mushroom Species.

Mushrooms are widely-consumed fungi which contain natural compounds that can be used both for their nutritive and medicinal properties, i.e., taking advantage of their antimicrobial, antiviral, antitumor, anti-allergic, immunomodulation, anti-inflammatory, anti-atherogenic, hypoglycemic, hepatoprotective and antioxidant effects. Currently, scientific interest in natural compounds extracted from the fungal species is increasing because these compounds are also known to have pharmacological/biological activity. Unfortunately, however, their mechanisms of action are often unknown, not well understood or have not been investigated in their entirety. Given the poly-pharmacological properties of bioactive fungal compounds, it was decided to carry out a multi-targeted approach to predict possible interactions occurring among bioactive natural fungal extracts and several macromolecular targets that are therapeutically interesting, i.e., proteins, enzymes and nucleic acids. A chemical database of compounds extracted from both edible and no-edible mushrooms was created. This database was virtually screened against 43 macromolecular targets downloaded from the Protein Data Bank website. The aim of this work is to provide a molecular description of the main interactions involving ligand/multi-target recognition in order to understand the polypharmacological profile of the most interesting fungal extracts and to suggest a design strategy of new multi-target agents.

Macrocyclic naphthalene diimides as G-quadruplex binders.

The synthesis, biological and molecular modeling evaluation of a series of macrocyclic naphthalene diimides is reported. The present investigation expands on the study of structure-activity relationships of prototype compound 2 by constraining the molecule into a macrocyclic structure with the aim of improving its G-quadruplex binding activity and selectivity. The new derivatives, compounds 4-7 carry spermidine- and spermine-like linkers while in compound 8 the inner basic nitrogen atoms of spermine have been replaced with oxygen atoms. The design strategy has led to potent compounds stabilizing both human telomeric (F21T) and c-KIT2 quadruplex sequences, and high selectivity for quadruplex in comparison to duplex DNA. Antiproliferative effects of the new derivatives 4-8 have been evaluated in a panel of cancer cell lines and all the tested compounds showed activity in the low micromolar or sub-micromolar range of concentrations. In order to rationalize the molecular basis of the DNA G-quadruplex versus duplex recognition preference, docking and molecular dynamics studies have been performed. The computational results support the observation that the main driving force in the recognition is due to electrostatic factors.

Structure-based virtual screening of novel natural alkaloid derivatives as potential binders of h-telo and c-myc DNA G-quadruplex conformations.

Several ligands can bind to the non-canonical G-quadruplex DNA structures thereby stabilizing them. These molecules can act as effective anticancer agents by stabilizing the telomeric regions of DNA or by regulating oncogene expression. In order to better interact with the quartets of G-quadruplex structures, G-binders are generally characterized by a large aromatic core involved in π-π stacking. Some natural flexible cyclic molecules from Traditional Chinese Medicine have shown high binding affinity with G-quadruplex, such as berbamine and many other alkaloids. Using the structural information available on G-quadruplex structures, we performed a high throughput in silico screening of commercially available alkaloid derivative databases by means of a structure-based approach based on docking and molecular dynamics simulations against the human telomeric sequence d[AG3(T2AG3)3] and the c-myc promoter structure. We identified 69 best hits reporting an improved theoretical binding affinity with respect to the active set. Among them, a berberine derivative, already known to remarkably inhibit telomerase activity, was related to a better theoretical affinity versus c-myc.

Exploring the DNA2-PNA heterotriplex formation in targeting the Bcl-2 gene promoter: A structural insight by physico-chemical and microsecond-scale MD investigation.

Peptide Nucleic Acids (PNAs) represent a promising tool for gene modulation in anticancer treatment. The uncharged peptidyl backbone and the resistance to chemical and enzymatic degradation make PNAs highly advantageous to form stable hybrid complexes with complementary DNA and RNA strands, providing higher stability than the corresponding natural analogues. Our and other groups' research has successfully shown that tailored PNA sequences can effectively downregulate the expression of human oncogenes using antigene, antisense, or anti-miRNA approaches. Specifically, we identified a seven bases-long PNA sequence, complementary to the longer loop of the main G-quadruplex structure formed by the bcl2midG4 promoter sequence, capable of downregulating the expression of the antiapoptotic Bcl-2 protein and enhancing the anticancer activity of an oncolytic adenovirus. Here, we extended the length of the PNA probe with the aim of including the double-stranded Bcl-2 promoter among the targets of the PNA probe. Our investigation primarily focused on the structural aspects of the resulting DNA2-PNA heterotriplex that were determined by employing conventional and accelerated microsecond-scale molecular dynamics simulations and chemical-physical analysis. Additionally, we conducted preliminary biological experiments using cytotoxicity assays on human A549 and MDA-MB-436 adenocarcinoma cell lines, employing the oncolytic adenovirus delivery strategy.

Bile acids serve as endogenous antagonists of the Leukemia inhibitory factor (LIF) receptor in oncogenesis.

The leukemia inhibitory factor (LIF) is member of interleukin (IL)-6 family of cytokines involved immune regulation, morphogenesis and oncogenesis. In cancer tissues, LIF binds a heterodimeric receptor (LIFR), formed by a LIFRß subunit and glycoprotein(gp)130, promoting epithelial mesenchymal transition and cell growth. Bile acids are cholesterol metabolites generated at the interface of host metabolism and the intestinal microbiota. Here we demonstrated that bile acids serve as endogenous antagonist to LIFR in oncogenesis. The tissue characterization of bile acids content in non-cancer and cancer biopsy pairs from gastric adenocarcinomas (GC) demonstrated that bile acids accumulate within cancer tissues, with glyco-deoxycholic acid (GDCA) functioning as negative regulator of LIFR expression. In patient-derived organoids (hPDOs) from GC patients, GDCA reverses LIF-induced stemness and proliferation. In summary, we have identified the secondary bile acids as the first endogenous antagonist to LIFR supporting a development of bile acid-based therapies in LIF-mediated oncogenesis.

Molecular and Structural Aspects of Clinically Relevant Mutations of SARS-CoV-2 RNA-Dependent RNA Polymerase in Remdesivir-Treated Patients.

(1) Background: SARS-CoV-2 RNA-dependent RNA polymerase (RdRp) is a promising therapeutic target to fight COVID-19, and many RdRp inhibitors nucleotide/nucleoside analogs, such as remdesivir, have been identified or are in clinical studies. However, the appearance of resistant mutations could reduce their efficacy. In the present work, we structurally evaluated the impact of RdRp mutations found at baseline in 39 patients treated with remdesivir and associated with a different degree of antiviral response in vivo. (2) Methods: A refined bioinformatics approach was applied to assign SARS-CoV-2 clade and lineage, and to define RdRp mutational profiles. In line with such a method, the same mutations were built and analyzed by combining docking and thermodynamics evaluations with both molecular dynamics and representative pharmacophore models. (3) Results: Clinical studies revealed that patients bearing the most prevalent triple mutant P323L+671S+M899I, which was present in 41% of patients, or the more complex mutational profile P323L+G671S+L838I+D738Y+K91E, which was found with a prevalence of 2.6%, showed a delayed reduced response to remdesivir, as confirmed by the increase in SARS-CoV-2 viral load and by a reduced theoretical binding affinity versus RdRp (ΔGbindWT = -122.70 kcal/mol; ΔGbindP323L+671S+M899I = -84.78 kcal/mol; ΔGbindP323L+G671S+L838I+D738Y+K91E = -96.74 kcal/mol). Combined computational approaches helped to rationalize such clinical observations, offering a mechanistic understanding of the allosteric effects of mutants on the global motions of the viral RNA synthesis machine and in the changes of the interactions patterns of remdesivir during its binding.

Discovery of BAR502, as potent steroidal antagonist of leukemia inhibitory factor receptor for the treatment of pancreatic adenocarcinoma.

Introduction: The leukemia inhibitory factor (LIF), is a cytokine belonging to IL-6 family, whose overexpression correlate with poor prognosis in cancer patients, including pancreatic ductal adenocarcinoma (PDAC). LIF signaling is mediate by its binding to the heterodimeric LIF receptor (LIFR) complex formed by the LIFR receptor and Gp130, leading to JAK1/STAT3 activation. Bile acids are steroid that modulates the expression/activity of membrane and nuclear receptors, including the Farnesoid-X-Receptor (FXR) and G Protein Bile Acid Activated Receptor (GPBAR1). Methods: Herein we have investigated whether ligands to FXR and GPBAR1 modulate LIF/LIFR pathway in PDAC cells and whether these receptors are expressed in human neoplastic tissues. Results: The transcriptome analysis of a cohort of PDCA patients revealed that expression of LIF and LIFR is increased in the neoplastic tissue in comparison to paired non-neoplastic tissues. By in vitro assay we found that both primary and secondary bile acids exert a weak antagonistic effect on LIF/LIFR signaling. In contrast, BAR502 a non-bile acid steroidal dual FXR and GPBAR1 ligand, potently inhibits binding of LIF to LIFR with an IC50 of 3.8 µM. Discussion: BAR502 reverses the pattern LIF-induced in a FXR and GPBAR1 independent manner, suggesting a potential role for BAR502 in the treatment of LIFR overexpressing-PDAC.

The leukemia inhibitory factor regulates fibroblast growth factor receptor 4 transcription in gastric cancer.

PURPOSE: The gastric adenocarcinoma (GC) represents the third cause of cancer-related mortality worldwide, and available therapeutic options remain sub-optimal. The Fibroblast growth factor receptors (FGFRs) are oncogenic transmembrane tyrosine kinase receptors. FGFR inhibitors have been approved for the treatment of various cancers and a STAT3-dependent regulation of FGFR4 has been documented in the H.pylori infected intestinal GC. Therefore, the modulation of FGFR4 might be useful for the treatment of GC. METHODS: To investigate wich factors could modulate FGFR4 signalling in GC, we employed RNA-seq analysis on GC patients biopsies, human patients derived organoids (PDOs) and cancer cell lines. RESULTS: We report that FGFR4 expression/function is regulated by the leukemia inhibitory factor (LIF) an IL-6 related oncogenic cytokine, in JAK1/STAT3 dependent manner. The transcriptomic analysis revealed a direct correlation between the expression of LIFR and FGFR4 in the tissue of an exploratory cohort of 31 GC and confirmed these findings by two external validation cohorts of GC. A LIFR inhibitor (LIR-201) abrogates STAT3 phosphorylation induced by LIF as well as recruitment of pSTAT3 to the promoter of FGFR4. Furthermore, inhibition of FGFR4 by roblitinib or siRNA abrogates STAT3 phosphorylation and oncogentic effects of LIF in GC cells, indicating that FGFR4 is a downstream target of LIF/LIFR complex. Treating cells with LIR-201 abrogates oncogenic potential of FGF19, the physiological ligand of FGFR4. CONCLUSIONS: Together these data unreveal a previously unregnized regulatory mechanism of FGFR4 by LIF/LIFR and demonstrate that LIF and FGF19 converge on the regulation of oncogenic STAT3 in GC cells.

GRID-based three-dimensional pharmacophores II: PharmBench, a benchmark data set for evaluating pharmacophore elucidation methods.

To date, published pharmacophore elucidation approaches typically use a handful of data sets for validation: here, we have assembled a data set for 81 targets, containing 960 ligands aligned using their cocrystallized protein targets, to provide the experimental "gold standard". The two-dimensional structures are also assembled to remove conformational bias; an ideal method would be able to take these structures as input, find the common features, and reproduce the bioactive conformations and their alignments to correspond with the X-ray-determined gold standard alignments. Here we present this data set and describe three objective measures to evaluate performance: the ability to identify the bioactive conformation, the ability to identify and correctly align this conformation for 50% of the molecules in each data set, and the pharmacophoric field similarity. We have applied this validation methodology to our pharmacophore elucidation method FLAPpharm, that is published in the first paper of this series and discuss the limitations of the data set and objective success criteria. Starting from two-dimensional structures and producing unbiased models, FLAPpharm was able to identify the bioactive conformations for 67% of the ligands and also to produce successful models according to the second metric for 67% of the Pharmbench data sets. Inspection of the unsuccessful models highlighted the limitation of this root mean square (rms)-derived metric, since many were found to be pharmacophorically reasonable, increasing the overall success rate to 83%. The PharmBench data set is available at http://www.moldiscovery.com/PharmBench , along with a web service to enable users to score model alignments coming from external methods in the same way that we have presented here and, therefore, establishes a pharmacophore elucidation benchmark data set available to be used by the community.

Chromene Derivatives as Selective TERRA G-Quadruplex RNA Binders with Antiproliferative Properties.

In mammalian cells, telomerase transcribes telomeres in large G-rich non-coding RNA, known as telomeric repeat-containing RNA (TERRA), which folds into noncanonical nucleic acid secondary structures called G-quadruplexes (G4s). Since TERRA G4 has been shown to be involved in telomere length and translation regulation, it could provide valuable insight into fundamental biological processes, such as cancer growth, and TERRA G4 binders could represent an innovative strategy for cancer treatment. In this work, the three best candidates identified in our previous virtual screening campaign on bimolecular DNA/RNA G4s were investigated on the monomolecular Tel DNA and TERRA G4s by means of molecular modelling simulations and in vitro and in cell analysis. The results obtained in this work highlighted the stabilizing power of all the three candidates on TERRA G4. In particular, the two compounds characterized by a chromene scaffold were selective TERRA G4 binders, while the compound with a naphthyridine core acted as a dual Tel/TERRA G4-binder. A biophysical investigation by circular dichroism confirmed the relative stabilization efficiency of the compounds towards TERRA and Tel G4s. The TERRA G4 stabilizing hits showed good antiproliferative activity against colorectal and lung adenocarcinoma cell lines. Lead optimization to increase TERRA G4 stabilization may provide new powerful tools against cancer.