Automatic identification of relevant chemical compounds from patents

In commercial research and development projects, public disclosure of new chemical compounds often takes place in patents. Only a small proportion of these compounds are published in journals, usually a few years after the patent. Patent authorities make available the patents but do not provide systematic continuous chemical annotations. Content databases such as Elsevier’s Reaxys provide such services mostly based on manual excerptions, which are time-consuming and costly. Automatic text-mining approaches help overcome some of the limitations of the manual process. Different text-mining approaches exist to extract chemical entities from patents. The majority of them have been developed using sub-sections of patent documents and focus on mentions of compounds. Less attention has been given to relevancy of a compound in a patent. Relevancy of a compound to a patent is based on the patent’s context. A relevant compound plays a major role within a patent. Identification of relevant compounds reduces the size of the extracted data and improves the usefulness of patent resources (e.g. supports identifying the main compounds). Annotators of databases like Reaxys only annotate relevant compounds. In this study, we design an automated system that extracts chemical entities from patents and classifies their relevance. The goldstandard set contained 18 789 chemical entity annotations. Of these, 10% were relevant compounds, 88% were irrelevant and 2% were equivocal. Our compound recognition system was based on proprietary tools. The performance (F-score) of the system on compound recognition was 84% on the development set and 86% on the test set. The relevancy classification system had an F-score of 86% on the development set and 82% on the test set. Our system can extract chemical compounds from patents and classify their relevance with high performance. This enables the extension of the Reaxys database by means of automation

Computational Analysis of Biological networks

Caratterizzare, descrivere ed estrarre informazioni da un network, \ue8 sicuramente uno dei principali obbiettivi della scienza, dato che lo studio dei network interessa differenti campi della ricerca, come la biologia, l'economia, le scienze sociali, l'informatica e cos\uec via. Ci\uf2 che si vuole \ue8 riuscire ad estrarre le propriet\ue0 fondamentali dei network e comprenderne la funzionalit\ue0. Questa tesi riguarda sia l'analisi topologica che l' analisi dinamica dei network biologici, anche se i risultati possono essere applicati a diversi campi. Per quanto riguarda l'analisi topologica viene utilizzato un approccio orientato ai nodi, utilizzando le centralit\ue0 per individuare i nodi pi\uf9 rilevanti e integrando tali risultati con dati da laboratorio. Viene inoltre descritto CentiScaPe, un software implementato per effettuare tale tipo di analisi. Vengono inoltre introdotti i concetti di "interference" e "robustness" che permettono di comprendere come un network si riarrangia in seguito alla rimozione o all'aggiunta di nodi. Per quanto riguarda l'analisi dinamica, si mostra come l'abstract interpretation pu\uf2 essere utilizzata nella simulazione di pathways per ottenere i risultati di migliaia di simulazioni in breve tempo e come possibile soluzione del problema della stima dei parametri mancanti.This thesis, treating both topological and dynamic points of view, concerns several aspects of biological networks analysis. Regarding the topological analysis of biological networks, the main contribution is the node-oriented point of view of the analysis. It means that instead of concentrating on global properties of the networks, we analyze them in order to extract properties of single nodes. An excellent method to face this problem is to use node centralities. Node centralities allow to identify nodes in a network having a relevant role in the network structure. This can not be enough if we are dealing with a biological network, since the role of a protein depends also on its biological activity that can be detected with lab experiments. Our approach is to integrate centralities analysis and data from biological experiments. A protocol of analysis have been produced, and the CentiScaPe tool for computing network centralities and integrating topological analysis with biological data have been designed and implemented. CentiScaPe have been applied to a human kino-phosphatome network and according to our protocol, kinases and phosphatases with highest centralities values have been extracted creating a new subnetwork of most central kinases and phosphatases. A lab experiment established which of this proteins presented high activation level and through CentiScaPe the proteins with both high centrality values and high activation level have been easily identified. The notion of node centralities interference have also been introduced to deal with central role of nodes in a biological network. It allow to identify which are the nodes that are more affected by the remotion of a particular node measuring the variation on their centralities values when such a node is removed from the network. The application of node centralities interference to the human kino-phosphatome revealed that different proteins affect centralities values of different nodes. Similarly to node centralities interference, the notion of centrality robustness of a node is introduced. This notion reveals if the central role of a node depends on other particular nodes in the network or if the node is ``robust'' in the sense that even if we remove or add other nodes the central role of the node remains almost unchanged. The dynamic aspects of biological networks analysis have been treated from an abstract interpretation point of view. Abstract interpretation is a powerful framework for the analysis of software and is excellent in deriving numerical properties of programs. Dealing with pathways, abstract interpretation have been adapted to the analysis of pathways simulation. Intervals domain and constants domain have been succesfully used to automatically extract information about reactants concentration. The intervals domain allow to determine the range of concentration of the proteins, and the constants domain have been used to know if a protein concentration become constant after a certain time. The other domain of analysis used is the congruences domain that, if applied to pathways simulation can easily identify regular oscillating behaviour in reactants concentration. The use of abstract interpretation allows to execute thousands of simulation and to completely and automatically characterize the behaviour of the pathways. In such a way it can be used also to solve the problem of parameters estimation where missing parameters can be detected with a brute force algorithm combined with the abstract interpretation analysis. The abstract interpretation approach have been succesfully applied to the mitotic oscillator pathway, characterizing the behaviour of the pathway depending on some reactants. To help the analysis of relation between reactants in the network, the notions of variables interference and variables abstract interference have been introduced and adapted to biological pathways simulation. They allow to find relations between properties of different reactants of the pathway. Using the abstract interference techniques we can say, for instance, which range of concentration of a protein can induce an oscillating behaviour of the pathway

Development and validation of in silico tools for efficient library design and data analysis in high throughput screening campaigns

My PhD project findings have their major application in the early phase of the drug discovery process, in particular we have developed and validated two computational tools (Molecular Assembles and LiGen) to support the hit finding and the hit to lead phases. I have reported here novel methods to first design chemical libraries optimized for HTS and then profile them for a specific target receptor or enzyme. I also analyzed the generated bio-chemical data in order to obtain robust SARs and to select the most promising hits for the follow up. The described methods support the iterative process of validated hit series optimization up to the identification of a lead. In chapter 3, Ligand generator (LiGen), a de novo tool for structure based virtual screening, is presented. The development of LiGen is a project based on a collaboration among Dompé Farmaceutici SpA, CINECA and the University of Parma. In this multidisciplinary group, the integration of different skills has allowed the development, from scratch, of a virtual screening tool, able to compete in terms of performance with long standing, well-established molecular docking tools such as Glide, Autodock and PLANTS. LiGen, using a novel docking algorithm, is able to perform ligand flexible docking without performing a conformational sampling. LiGen also has other distinctive features with respect to other molecular docking programs: • LiGen uses the inverse pharmacophore derived from the binding site to identify the putative bioactive conformation of the molecules, thus avoiding the evaluation of molecular conformations which do not match the key features of the binding site. • LiGen implemenst a de novo molecule builder based on the accurate definition of chemical rules taking account of building block (reagents) reactivity. • LiGen is natively a multi-platform C++ portable code designed for HPC applications and optimized for the most recent hardware architectures like the Xeon Phi Accelerators. Chapter 3 also reports the further development and optimization of the software starting from the results obtained in the first optimization step performed to validate the software and to derive the default parameters. In chapter 4, the application of LiGen in the discovery and optimization of novel inhibitors of the complement factor 5 receptor (C5aR) is reported. Briefly, the C5a anaphylatoxin acting on its cognate G protein-coupled receptor C5aR is a potent pronociceptive mediator in several models of inflammatory and neuropathic pain. Although there has long been interest in the identification of C5aR inhibitors, their development has been complicated, as is the case with many peptidomimetic drugs, mostly due to the poor drug-like properties of these molecules. Herein, we report the de novo design of a potent and selective C5aR noncompetitive allosteric inhibitor, DF2593A. DF2593A design was guided by the hypothesis that an allosteric site, the “minor pocket”, previously characterized in CXCR1 and CXCR2, could be functionally conserved in the GPCR class.DF2593A potently inhibited C5a-induced migration of human and rodent neutrophils in vitro. Moreover, oral administration of DF2593A effectively reduced mechanical hyperalgesia in several models of acute and chronic inflammatory and neuropathic pain in vivo, without any apparent side effects. Chapter 5 describes another tool: Molecular Assemblies (MA), a novel metrics based on a hierarchical representation of the molecule based on different representations of the scaffold of the molecule and pruning rules. The algorithm used by MA, defining a priori a metrics (a set of rules), creates a representation of the chemical structure through hierarchical decomposition of the scaffold in fragments, in a pathway invariant way (this feature is novel with respect to the other algorithms reported in literature). Such structure decomposition is applied to nine hierarchical representation of the scaffold of the reference molecule, differing for the content of structural information: atom typing and bond order (this feature is novel with respect to the other algorithms reported in literature) The algorithm (metrics) generates a multi-dimensional hierarchical representation of the molecule. This descriptor applied to a library of compounds is able to extract structural (molecule having the same scaffold, wireframe or framework) and sub structural (molecule having the same fragments in common) relations among all the molecules. At least, this method generates relations among molecules based on identities (scaffolds or fragments). Such an approach produces a unique representation of the reference chemical space not biased by the threshold used to define the similarity cut-off between two molecules. This is in contrast to other methods which generate representations based in similarities. MA procedure, retrieving all scaffold representation, fragments and fragmentation’s patterns (according to the predefined rules) from a molecule, creates a molecular descriptor useful for several cheminformatics applications: • Visualization of the chemical space. The scaffold relations (Figure 7) and the fragmentation patterns can be plotted using a network representation. The obtained graphs are useful depictions of the chemical space highlighting the relations that occur among the molecule in a two dimensional space. • Clustering of the chemical space. The relations among the molecules are based on identities. This means that the scaffold representations and their fragments can be used as a hierarchical clustering method. This descriptor produces clusters that are independent from the number and similarity among closest neighbors because belonging to a cluster is a property of the single molecule (Figure 8). This intrinsic feature makes the scaffold based clustering much faster than other methods in producing “stable” clusters in fact, adding and removing molecules increases and decreases the number of clusters while adding or removing relations among the clusters. However these changes do not affect the cluster number and the relation of the other molecules in dataset. • Generate scaffold-based fingerprints. The descriptor can be used as a fingerprint of the molecule and to generate a similarity index able to compare single molecules or also to compare the diversity of two libraries as a whole. Chapter 6 reports an application of MA in the design of a diverse drug-like scaffold based library optimized for HTS campaigns. A well designed, sizeable and properly organized chemical library is a fundamental prerequisite for any HTS project. To build a collection of chemical compounds with high chemical diversity was the aim of the Italian Drug Discovery Network (IDDN) initiative. A structurally diverse collection of about 200,000 chemical molecules was designed and built taking into account practical aspects related to experimental HTS procedures. Algorithms and procedures were developed and implemented to address compound filtering, selection, clusterization and plating. Chapter 7 collects concluding remarks and plans for the further development of the tools

Text Mining for Chemical Compounds

Exploring the chemical and biological space covered by patent and journal publications is crucial in early- stage medicinal chemistry activities. The analysis provides understanding of compound prior art, novelty checking, validation of biological assays, and identification of new starting points for chemical exploration. Extracting chemical and biological entities from patents and journals through manual extraction by expert curators can take substantial amount of time and resources. Text mining methods can help to ease this process. In this book, we addressed the lack of quality measurements for assessing the correctness of structural representation within and across chemical databases; lack of resources to build text-mining systems; lack of high performance systems to extract chemical compounds from journals and patents; and lack of automated systems to identify relevant compounds in patents. The consistency and ambiguity of chemical identifiers was analyzed within and between small- molecule databases in Chapter 2 and Chapter 3. In Chapter 4 and Chapter 7 we developed resources to enable the construction of chemical text-mining systems. In Chapter 5 and Chapter 6, we used community challenges (BioCreative V and BioCreative VI) and their corresponding resources to identify mentions of chemical compounds in journal abstracts and patents. In Chapter 7 we used our findings in previous chapters to extract chemical named entities from patent full text and to classify the relevancy of chemical compounds