Does Ligand Symmetry Play a Role in the Stabilization of DNA G-Quadruplex Host-Guest Complexes?
In efforts to find agents with improved biological activity against cancer cells, recent years have seen an increased interest in the study of small molecules able to bind the deoxyribonucleic acid (DNA) when it assumes secondary structures known as G-quadruplexes (G4s) preferring them over the B form. Currently, several compounds reported in literature have already shown to be good candidates as G4s DNA stabilizers. Even though some specific features for the G4s affinity are known, such as a π-delocalized system able to stack at the top/end of a G-tetrad and positively charged substituents able to interact with the grooves, it is not clear yet what kind of structural features affect more t…
Micelles of the chiral biocompatible surfactant (1R,2S)-dodecyl(2-hydroxy-1-methyl-2-phenylethyl)dimethylammonium bromide (DMEB): molecular dynamics and fragmentation patterns in the gas phase.
RATIONALE: The study of self-assembly processes of surfactant molecules in the gas phase is of great interest for several theoretical and technological reasons related to their possible exploitation as drug carriers, protein shields and cleaning agents in the gas phase. METHODS: The stability and fragmentation patterns of singly and multiply charged (either positively or negatively) aggregates of the surfactant (1R,2S)-dodecyl(2-hydroxy-1-methyl-2-phenylethyl)dimethyl ammonium bromide (DMEB) in the gas phase have been studied by ion mobility mass spectrometry and tandem mass spectrometry. Molecular dynamics (MD) simulations of positively and negatively singly and multiply charged DMEB aggre…
Drugs Polypharmacology by in Silico Methods: New Opportunities in Drug Discovery
Background Polypharmacology, defined as the modulation of multiple proteins rather than a single target to achieve a desired therapeutic effect, has been gaining increasing attention since 1990s, when industries had to withdraw several drugs due to their adverse effects, leading to permanent injuries or death, with multi-billiondollar legal damages. Therefore, if up to then the "one drug one target" paradigm had seen many researchers interest focused on the identification of selective drugs, with the strong expectation to avoid adverse drug reactions (ADRs), very recently new research strategies resulted more appealing even as attempts to overcome the decline in productivity of the drug dis…
2-methoxyestradiol impacts on amino acids-mediated metabolic reprogramming in osteosarcoma cells by interaction with NMDA receptor
Deregulation of serine and glycine metabolism, have been identified to function as metabolic regulators in supporting tumor cell growth. The role of serine and glycine in regulation of cancer cell proliferation is complicated, dependent on concentrations of amino acids and tissue-specific. D-serine and glycine are coagonists of N-methyl-D-aspartate receptor subunit GRIN1. Importantly, NMDA receptors are widely expressed in cancer cells and play an important role in regulation of cell death, proliferation and metabolism of numerous malignancies. The aim of the present work was to associate the metabolism of glycine and D-serine with the anticancer activity of 2-methoxyestradiol. 2-methoxyest…
Targeting SARS-CoV-2 RBD Interface: a Supervised Computational Data-Driven Approach to Identify Potential Modulators
Coronavirus disease 2019 (COVID-19) has spread out as a pandemic threat affecting over 2 million people. The infectious process initiates via binding of SARS-CoV-2 Spike (S) glycoprotein to host angiotensin-converting enzyme 2 (ACE2). The interaction is mediated by the receptor-binding domain (RBD) of S glycoprotein, promoting host receptor recognition and binding to ACE2 peptidase domain (PD), thus representing a promising target for therapeutic intervention. Herein, we present a computational study aimed at identifying small molecules potentially able to target RBD. Although targeting PPI remains a challenge in drug discovery, our investigation highlights that interaction between SARS-CoV…
EMBER—Embedding Multiple Molecular Fingerprints for Virtual Screening
In recent years, the debate in the field of applications of Deep Learning to Virtual Screening has focused on the use of neural embeddings with respect to classical descriptors in order to encode both structural and physical properties of ligands and/or targets. The attention on embeddings with the increasing use of Graph Neural Networks aimed at overcoming molecular fingerprints that are short range embeddings for atomic neighborhoods. Here, we present EMBER, a novel molecular embedding made by seven molecular fingerprints arranged as different “spectra” to describe the same molecule, and we prove its effectiveness by using deep convolutional architecture that assesses ligands&…
Dynamic-shared Pharmacophore Approach as Tool to Design New Allosteric PRC2 Inhibitors, Targeting EED Binding Pocket.
Abstract: The Polycomb Repressive complex 2 (PRC2) maintains a repressive chromatin state and silences many genes, acting as methylase on histone tails. This enzyme was found overexpressed in many types of cancer. In this work, we have set up a Computer-Aided Drug Design approach based on the allosteric modulation of PRC2. In order to minimize the possible bias derived from using a single set of coordinates within the protein-ligand complex, a dynamic workflow was developed. In details, molecular dynamic was used as tool to identify the most significant ligand-protein interactions from several crystallized protein structures. The identified features were used for the creation of dynamic pha…
A dynamic multiple receptor conformations (MD-MRC) approach to enhance early enrichment in virtual screening. A case study on PPAR-alpha
A dynamic-common pharmacophore approach to improve virtual screening. A case study on PPAR-alpha
Micelles, Rods, Liposomes, and Other Supramolecular Surfactant Aggregates: Computational Approaches
Surfactants are an interesting class of compounds characterized by the segregation of polar and apolar domains in the same molecule. This peculiarity makes possible a whole series of microscopic and macroscopic effects. Among their features, their ability to segregate particles (fluids or entire domains) and to reduce the surface/interfacial tension is the utmost important. The interest in the chemistry of surfactants never weakened; instead, waves of increasing interest have occurred every time a new field of application of these molecules has been discovered. All these special characteristics depend largely on the ability of surfactants to self-assemble and constitute supramolecular struc…
Conf-VLKA: A structure-based revisitation of the Virtual Lock-and-key Approach
In a previous work, we developed the in house Virtual Lock-and-Key Approach (VLKA) in order to evaluate target assignment starting from molecular descriptors calculated on known inhibitors used as an information source. This protocol was able to predict the correct biological target for the whole dataset with a good degree of reliability (80%), and proved experimentally, which was useful for the target fishing of unknown compounds. In this paper, we tried to remodel the previous in house developed VLKA in a more sophisticated one in order to evaluate the influence of 3D conformation of ligands on the accuracy of the prediction. We applied the same previous algorithm of scoring and ranking b…
In Silico Insights towards the Identification of NLRP3 Druggable Hot Spots
NLRP3 (NOD-like receptor family, pyrin domain-containing protein 3) activation has been linked to several chronic pathologies, including atherosclerosis, type-II diabetes, fibrosis, rheumatoid arthritis, and Alzheimer’s disease. Therefore, NLRP3 represents an appealing target for the development of innovative therapeutic approaches. A few companies are currently working on the discovery of selective modulators of NLRP3 inflammasome. Unfortunately, limited structural data are available for this target. To date, MCC950 represents one of the most promising noncovalent NLRP3 inhibitors. Recently, a possible region for the binding of MCC950 to the NLRP3 protein was described but no details were …
3-(6-Phenylimidazo [2,1-b][1,3,4]thiadiazol-2-yl)-1H-Indole Derivatives as New Anticancer Agents in the Treatment of Pancreatic Ductal Adenocarcinoma
A new series of imidazo[2,1-b][1,3,4]thiadiazole derivatives was efficiently synthesized and screened for their in vitro antiproliferative activity on a panel of pancreatic ductal adenocarcinoma (PDAC) cells, including SUIT-2, Capan-1 and Panc-1. Compounds 9c and 9l, showed relevant in vitro antiproliferative activity on all three pre-clinical models with half maximal inhibitory concentration (IC50) ranging from 5.11 to 10.8 µ
Design, synthesis and preliminary evaluation of dopamine-amino acid conjugates as potential D1 dopaminergic modulators.
Abstract The dopamine-amino acid conjugate DA-Phen was firstly designed to obtain a useful prodrug for the therapy of Parkinson's disease, but experimental evidence shows that it effectively interacts with D1 dopamine receptors (D1DRs), leading to an enhancement in cognitive flexibility and to the development of adaptive strategies in aversive mazes, together with a decrease in despair-like behavior. In this paper, homology modelling, molecular dynamics, and site mapping of D1 receptor were carried out with the aim of further performing docking studies on other dopamine conjugates compared with D1 agonists, in the attempt to identify new compounds with potential dopaminergic activity. Two n…
Evaluating the stability of pharmacophore features using molecular dynamics simulations.
Abstract Molecular dynamics simulations of twelve protein—ligand systems were used to derive a single, structure based pharmacophore model for each system. These merged models combine the information from the initial experimental structure and from all snapshots saved during the simulation. We compared the merged pharmacophore models with the corresponding PDB pharmacophore models, i.e., the static models generated from an experimental structure in the usual manner. The frequency of individual features, of feature types and the occurrence of features not present in the static model derived from the experimental structure were analyzed. We observed both pharmacophore features not visible in …
Convolutional architectures for virtual screening
Abstract Background A Virtual Screening algorithm has to adapt to the different stages of this process. Early screening needs to ensure that all bioactive compounds are ranked in the first positions despite of the number of false positives, while a second screening round is aimed at increasing the prediction accuracy. Results A novel CNN architecture is presented to this aim, which predicts bioactivity of candidate compounds on CDK1 using a combination of molecular fingerprints as their vector representation, and has been trained suitably to achieve good results as regards both enrichment factor and accuracy in different screening modes (98.55% accuracy in active-only selection, and 98.88% …
A Molecular Dynamics-Shared Pharmacophore Approach to Boost Early-Enrichment Virtual Screening: A Case Study on Peroxisome Proliferator-Activated Receptor α
Molecular dynamics (MD) simulations can be used, prior to virtual screening, to add flexibility to proteins and study them in a dynamic way. Furthermore, the use of multiple crystal structures of the same protein containing different co-crystallized ligands can help elucidate the role of the ligand on a protein's active conformation, and then explore the most common interactions between small molecules and the receptor. In this work, we evaluated the contribution of the combined use of MD on crystal structures containing the same protein but different ligands to examine the crucial ligand-protein interactions within the complexes. The study was carried out on peroxisome proliferator-activat…
Pharmacophore Models Derived from Molecular Dynamics Simulations of Protein-Ligand Complexes: A Case Study
A single, merged pharmacophore hypothesis is derived combining 2000 pharmacophore models obtained during a 20 ns molecular dynamics simulation of a protein-ligand complex with one pharmacophore model derived from the initial PDB structure. This merged pharmacophore model contains all features that are present during the simulation and statistical information about the dynamics of the pharmacophore features. Based on the dynamics of the pharmacophore features we derive two distinctive feature patterns resulting in two different pharmacophore models for the analyzed system – the first model consists of features that are obtained from the PDB structure and the second uses two features that ca…
A convolutional neural network for virtual screening of molecular fingerprints
In the last few years, Deep Learning (DL) gained more and more impact on drug design because it allows a huge increase of the prediction accuracy in many stages of such a complex process. In this paper a Virtual Screening (VS) procedure based on Convolutional Neural Networks (CNN) is presented, that is aimed at classifying a set of candidate compounds as regards their biological activity on a particular target protein. The model has been trained on a dataset of active/inactive compounds with respect to the Cyclin-Dependent Kinase 1 (CDK1) a very important protein family, which is heavily involved in regulating the cell cycle. One qualifying point of the proposed approach is the use of molec…
The Repurposing of Old Drugs or Unsuccessful Lead Compounds by in Silico Approaches: New Advances and Perspectives
Have you a compound in your lab, which was not successful against the designed target, or a drug that is no more attractive? The drug repurposing represents the right way to reconsider them. It can be defined as the modern and rationale approach of the traditional methods adopted in drug discovery, based on the knowledge, insight and luck, alias known as serendipity. This repurposing approach can be applied both in silico and in wet. In this review we report the molecular modeling facilities that can be of huge support in the repurposing of drugs and/or unsuccessful lead compounds. In the last decades, different methods were proposed to help the scientists in drug design and in drug repurpo…
Synthesis, antiproliferative activity, and in silico insights of new 3-benzoylamino-benzo[ b ]thiophene derivatives
A new series of 3-benzoylamino-5-imidazol-5-yl-benzo[b]thiophenes and the parent amino derivatives were synthesized and screened as antitumor agents. All tested compounds showed concentration-dependent antiproliferative activity profile against HeLa cell line, exhibiting GI50 values in the low micromolar range. The most active compounds were tested in cell cycle perturbation experiments. A rapid accumulation of cells in the G2/M phase, with a concomitant reduction of cells in both the S and G0/G1 phases, was observed, suggesting that cell exposure to selected derivatives produces mitotic failure. To rationalize the biological results, the 3-benzoylamino-benzo[b]thiophenes were analyzed thro…
SF3B1 modulators affect key genes in metastasis and drug influx: a new approach to fight pancreatic cancer chemoresistance.
Aim: Because mutations of splicing factor 3B subunit-1 (SF3B1) have been identified in 4% of pancreatic ductal adenocarcinoma (PDAC) patients, we investigated the activity of new potential inhibitors of SF3B1 in combination with gemcitabine, one of the standard drugs, in PDAC cell lines. Methods: One imidazo[2,1-b][1,3,4]thiadiazole derivative (IS1) and three indole derivatives (IS2, IS3 and IS4), selected by virtual screening from an in-house library, were evaluated by the sulforhodamine-B and wound healing assay for their cytotoxic and antimigratory activity in the PDAC cells SUIT-2, Hs766t and Panc05.04, the latter harbouring the SF3B1 mutations. The effects on the splicing pattern of pr…
Exploring the SARS-CoV-2 Proteome in the Search of Potential Inhibitors via Structure-based Pharmacophore Modeling/Docking Approach
To date, SARS-CoV-2 infectious disease, named COVID-19 by the World Health Organization (WHO) in February 2020, has caused millions of infections and hundreds of thousands of deaths. Despite the scientific community efforts, there are currently no approved therapies for treating this coronavirus infection. The process of new drug development is expensive and time-consuming, so that drug repurposing may be the ideal solution to fight the pandemic. In this paper, we selected the proteins encoded by SARS-CoV-2 and using homology modeling we identified the high-quality model of proteins. A structure-based pharmacophore modeling study was performed to identify the pharmacophore features for each…
2-Methoxyestradiol Affects Mitochondrial Biogenesis Pathway and Succinate Dehydrogenase Complex Flavoprotein Subunit A in Osteosarcoma Cancer Cells.
Background/aim Dysregulation of mitochondrial pathways is implicated in several diseases, including cancer. Notably, mitochondrial respiration and mitochondrial biogenesis are favored in some invasive cancer cells, such as osteosarcoma. Hence, the aim of the current work was to investigate the effects of 2-methoxyestradiol (2-ME), a potent anticancer agent, on the mitochondrial biogenesis of osteosarcoma cells. Materials and methods Highly metastatic osteosarcoma 143B cells were treated with 2-ME separately or in combination with L-lactate, or with the solvent (non-treated control cells). Protein levels of α-syntrophin and peroxisome proliferator-activated receptor gamma, coactivator 1 alph…
Micelles of the chiral biocompatible surfactant (1R,2S)-dodecyl(2-hydroxy-1-methyl-2-phenylethyl)dimethylammonium bromide (DMEB): molecular dynamics and fragmentation patterns in the gas phase
Rationale The study of self-assembly process of surfactant molecules in gas phase is of actually interest for several theoretical and technological reasons related to their possible exploitation as drug carriers, protein shields and cleaning agents in gas phase. Methods Stability and fragmentation patterns of singly and multiply charged (either positively or negatively) aggregates of the surfactant (1R,2S)-dodecyl(2-hydroxy-1-methyl-2-phenylethyl) dimethyl ammonium bromide (DMEB) in gas phase have been studied by ion mobility mass spectrometry and tandem mass spectrometry. Molecular dynamics (MD) simulations of positively and negatively singly and multiply charged DMEB aggregates have been …
Common Hits Approach: Combining Pharmacophore Modeling and Molecular Dynamics Simulations.
We present a new approach that incorporates flexibility based on extensive MD simulations of protein-ligand complexes into structure-based pharmacophore modeling and virtual screening. The approach uses the multiple coordinate sets saved during the MD simulations and generates for each frame a pharmacophore model. Pharmacophore models with the same pharmacophore features are pooled. In this way the high number of pharmacophore models that results from the MD simulation is reduced to only a few hundred representative pharmacophore models. Virtual screening runs are performed with every representative pharmacophore model; the screening results are combined and rescored to generate a single hi…
An overview of recent molecular dynamics applications as medicinal chemistry tools for the undruggable site challenge
Molecular dynamics (MD) has become increasingly popular due to the development of hardware and software solutions and the improvement in algorithms, which allowed researchers to scale up calculations in order to speed them up. MD simulations are usually used to address protein folding issues or protein-ligand complex stability through energy profile analysis over time. In recent years, the development of new tools able to deeply explore a potential energy surface (PES) has allowed researchers to focus on the dynamic nature of the binding recognition process and binding-induced protein conformational changes. Moreover, modern approaches have been demonstrated to be effective and reliable in …
Development and optimisation of computational tools for drug discovery
The aim of my PhD project was the development, optimisation, and implementation of new in silico virtual screening protocols. Specifically, this thesis manuscript is divided into three main parts, presenting some of the papers published during my doctoral work. The first one, here named CHEMOMETRIC PROTOCOLS IN DRUG DISCOVERY, is about the optimisation and application of an in house developed chemometric protocol. This part has been entirely developed at the University of Palermo - STEBICEF Department - under the guide of my supervisors. During the development of this part I have personally worked on the tuning and optimisation of the algorithm and on the docking campaigns to obtain molecul…