Cheminformatic studies of molecular properties and their influence on bioactivity in the context of drug discovery

Molecular shape has long been known to be an important property for the process of molecular recognition. Previous studies postulated the existence of a drug-like shape space that could be used to artificially bias the composition of screening libraries, with the aim to increase the chance of success in Hit Identification. In this work, it was analysed to which extend this assumption holds true. Normalized Principal Moments of Inertia Ratios (NPRs) have been used to describe the molecular shape of small molecules. It was investigated, whether active molecules of diverse targets are located in preferred subspaces of the NPR shape space. Results illustrated a significantly stronger clustering than could be expected by chance, with parts of the space unlikely to be occupied by active compounds. Furthermore, a strong enrichment of elongated, rather flat shapes could be observed, while globular compounds were highly underrepresented. This was confirmed for a wide range of small molecule datasets from different origins. Active compounds exhibited a high overlap in their shape distributions across different targets, making a purely shape­ based discrimination very difficult. An additional perspective was provided by comparing the shapes of protein binding pockets with those of their respective ligands. Although more globular than their ligands, it was observed that binding sites shapes exhibited a similarly skewed distribution in shape space: spherical shapes were highly underrepresented. This was different for unoccupied binding pockets of smaller size. These were on the contrary identified to possess a more globular shape. The relation between shape complementarity and exhibited bioactivity was analysed; a moderate correlation between bioactivity and parameters including pocket coverage, distance in shape space, and others could be identified, which reflects the importance of shape complementarity. However, this also suggests that other aspects are of relevance for molecular recognition. A subsequent analysis assessed if and how shape and volume information retrieved from pocket or respective reference ligands could be used as a pre-filter in a virtual screening approach. ln Lead Optimization compounds need to get optimized with respect to a variety of pararneters. Here, the availability of past success stories is very valuable, as they can guide medicinal chemists during their analogue synthesis plans. However, although of tremendous interest for the public domain, so far only large corporations had the ability to mine historical knowledge in their proprietary databases. With the aim to provide such information, the SwissBioisostere database was developed and released during this thesis. This database contains information on 21,293,355 performed substructural exchanges, corresponding to 5,586,462 unique replacements that have been measured in 35,039 assays against 1,948 molecular targets representing 30 target classes, and on their impact on bioactivity . A user-friendly interface was developed that provides facile access to these data and is accessible at http//www.swissbioisostere.ch. The ChEMBL database was used as primary data source of bioactivity information. Matched molecular pairs have been identified in the extracted and cleaned data. Success-based scores were developed and integrated into the database to allow re-ranking of proposed replacements by their past outcomes. It was analysed to which degree these scores correlate with chemical similarity of the underlying fragments. An unexpectedly weak relationship was detected and further investigated. Use cases of this database were envisioned, and functionalities implemented accordingly: replacement outcomes are aggregatable at the assay level, and it was shawn that an aggregation at the target or target class level could also be performed, but should be accompanied by a careful case-by-case assessment. It was furthermore observed that replacement success depends on the activity of the starting compound A within a matched molecular pair A-B. With increasing potency the probability to lose bioactivity through any substructural exchange was significantly higher than in low affine binders. A potential existence of a publication bias could be refuted. Furthermore, often performed medicinal chemistry strategies for structure-activity-relationship exploration were analysed using the acquired data. Finally, data originating from pharmaceutical companies were compared with those reported in the literature. It could be seen that industrial medicinal chemistry can access replacement information not available in the public domain. In contrast, a large amount of often-performed replacements within companies could also be identified in literature data. Preferences for particular replacements differed between these two sources. The value of combining different endpoints in an evaluation of molecular replacements was investigated. The performed studies highlighted furthermore that there seem to exist no universal substructural replacement that always retains bioactivity irrespective of the biological environment. A generalization of bioisosteric replacements seems therefore not possible. - La forme tridimensionnelle des molécules a depuis longtemps été reconnue comme une propriété importante pour le processus de reconnaissance moléculaire. Des études antérieures ont postulé que les médicaments occupent préférentiellement un sous-ensemble de l'espace des formes des molécules. Ce sous-ensemble pourrait être utilisé pour biaiser la composition de chimiothèques à cribler, dans le but d'augmenter les chances d'identifier des Hits. L'analyse et la validation de cette assertion fait l'objet de cette première partie. Les Ratios de Moments Principaux d'Inertie Normalisés (RPN) ont été utilisés pour décrire la forme tridimensionnelle de petites molécules de type médicament. Il a été étudié si les molécules actives sur des cibles différentes se co-localisaient dans des sous-espaces privilégiés de l'espace des formes. Les résultats montrent des regroupements de molécules incompatibles avec une répartition aléatoire, avec certaines parties de l'espace peu susceptibles d'être occupées par des composés actifs. Par ailleurs, un fort enrichissement en formes allongées et plutôt plates a pu être observé, tandis que les composés globulaires étaient fortement sous-représentés. Cela a été confirmé pour un large ensemble de compilations de molécules d'origines différentes. Les distributions de forme des molécules actives sur des cibles différentes se recoupent largement, rendant une discrimination fondée uniquement sur la forme très difficile. Une perspective supplémentaire a été ajoutée par la comparaison des formes des ligands avec celles de leurs sites de liaison (poches) dans leurs protéines respectives. Bien que plus globulaires que leurs ligands, il a été observé que les formes des poches présentent une distribution dans l'espace des formes avec le même type d'asymétrie que celle observée pour les ligands: les formes sphériques sont fortement sous­ représentées. Un résultat différent a été obtenu pour les poches de plus petite taille et cristallisées sans ligand: elles possédaient une forme plus globulaire. La relation entre complémentarité de forme et bioactivité a été également analysée; une corrélation modérée entre bioactivité et des paramètres tels que remplissage de poche, distance dans l'espace des formes, ainsi que d'autres, a pu être identifiée. Ceci reflète l'importance de la complémentarité des formes, mais aussi l'implication d'autres facteurs. Une analyse ultérieure a évalué si et comment la forme et le volume d'une poche ou de ses ligands de référence pouvaient être utilisés comme un pré-filtre dans une approche de criblage virtuel. Durant l'optimisation d'un Lead, de nombreux paramètres doivent être optimisés simultanément. Dans ce contexte, la disponibilité d'exemples d'optimisations réussies est précieuse, car ils peuvent orienter les chimistes médicinaux dans leurs plans de synthèse par analogie. Cependant, bien que d'un extrême intérêt pour les chercheurs dans le domaine public, seules les grandes sociétés pharmaceutiques avaient jusqu'à présent la capacité d'exploiter de telles connaissances au sein de leurs bases de données internes. Dans le but de remédier à cette limitation, la base de données SwissBioisostere a été élaborée et publiée dans le domaine public au cours de cette thèse. Cette base de données contient des informations sur 21 293 355 échanges sous-structuraux observés, correspondant à 5 586 462 remplacements uniques mesurés dans 35 039 tests contre 1948 cibles représentant 30 familles, ainsi que sur leur impact sur la bioactivité. Une interface a été développée pour permettre un accès facile à ces données, accessible à http:/ /www.swissbioisostere.ch. La base de données ChEMBL a été utilisée comme source de données de bioactivité. Une version modifiée de l'algorithme de Hussain et Rea a été implémentée pour identifier les Matched Molecular Pairs (MMP) dans les données préparées au préalable. Des scores de succès ont été développés et intégrés dans la base de données pour permettre un reclassement des remplacements proposés selon leurs résultats précédemment observés. La corrélation entre ces scores et la similarité chimique des fragments correspondants a été étudiée. Une corrélation plus faible qu'attendue a été détectée et analysée. Différents cas d'utilisation de cette base de données ont été envisagés, et les fonctionnalités correspondantes implémentées: l'agrégation des résultats de remplacement est effectuée au niveau de chaque test, et il a été montré qu'elle pourrait également être effectuée au niveau de la cible ou de la classe de cible, sous réserve d'une analyse au cas par cas. Il a en outre été constaté que le succès d'un remplacement dépend de l'activité du composé A au sein d'une paire A-B. Il a été montré que la probabilité de perdre la bioactivité à la suite d'un remplacement moléculaire quelconque est plus importante au sein des molécules les plus actives que chez les molécules de plus faible activité. L'existence potentielle d'un biais lié au processus de publication par articles a pu être réfutée. En outre, les stratégies fréquentes de chimie médicinale pour l'exploration des relations structure-activité ont été analysées à l'aide des données acquises. Enfin, les données provenant des compagnies pharmaceutiques ont été comparées à celles reportées dans la littérature. Il a pu être constaté que les chimistes médicinaux dans l'industrie peuvent accéder à des remplacements qui ne sont pas disponibles dans le domaine public. Par contre, un grand nombre de remplacements fréquemment observés dans les données de l'industrie ont également pu être identifiés dans les données de la littérature. Les préférences pour certains remplacements particuliers diffèrent entre ces deux sources. L'intérêt d'évaluer les remplacements moléculaires simultanément selon plusieurs paramètres (bioactivité et stabilité métabolique par ex.) a aussi été étudié. Les études réalisées ont souligné qu'il semble n'exister aucun remplacement sous-structural universel qui conserve toujours la bioactivité quel que soit le contexte biologique. Une généralisation des remplacements bioisostériques ne semble donc pas possible.

Prominent use of distal 5' transcription start sites and discovery of a large number of additional exons in ENCODE regions.

This report presents systematic empirical annotation of transcript products from 399 annotated protein-coding loci across the 1% of the human genome targeted by the Encyclopedia of DNA elements (ENCODE) pilot project using a combination of 5' rapid amplification of cDNA ends (RACE) and high-density resolution tiling arrays. We identified previously unannotated and often tissue- or cell-line-specific transcribed fragments (RACEfrags), both 5' distal to the annotated 5' terminus and internal to the annotated gene bounds for the vast majority (81.5%) of the tested genes. Half of the distal RACEfrags span large segments of genomic sequences away from the main portion of the coding transcript and often overlap with the upstream-annotated gene(s). Notably, at least 20% of the resultant novel transcripts have changes in their open reading frames (ORFs), most of them fusing ORFs of adjacent transcripts. A significant fraction of distal RACEfrags show expression levels comparable to those of known exons of the same locus, suggesting that they are not part of very minority splice forms. These results have significant implications concerning (1) our current understanding of the architecture of protein-coding genes; (2) our views on locations of regulatory regions in the genome; and (3) the interpretation of sequence polymorphisms mapping to regions hitherto considered to be "noncoding," ultimately relating to the identification of disease-related sequence alterations.

A steroidomic approach for biomarkers discovery in doping control.

Anti-doping authorities have high expectations of the athlete steroidal passport (ASP) for anabolic-androgenic steroids misuse detection. However, it is still limited to the monitoring of known well-established compounds and might greatly benefit from the discovery of new relevant biomarkers candidates. In this context, steroidomics opens the way to the untargeted simultaneous evaluation of a high number of compounds. Analytical platforms associating the performance of ultra-high pressure liquid chromatography (UHPLC) and the high mass-resolving power of quadrupole time-of-flight (QTOF) mass spectrometers are particularly adapted for such purpose. An untargeted steroidomic approach was proposed to analyse urine samples from a clinical trial for the discovery of relevant biomarkers of testosterone undecanoate oral intake. Automatic peak detection was performed and a filter of reference steroid metabolites mass-to-charge ratio (m/z) values was applied to the raw data to ensure the selection of a subset of steroid-related features. Chemometric tools were applied for the filtering and the analysis of UHPLC-QTOF-MS(E) data. Time kinetics could be assessed with N-way projections to latent structures discriminant analysis (N-PLS-DA) and a detection window was confirmed. Orthogonal projections to latent structures discriminant analysis (O-PLS-DA) classification models were evaluated in a second step to assess the predictive power of both known metabolites and unknown compounds. A shared and unique structure plot (SUS-plot) analysis was performed to select the most promising unknown candidates and receiver operating characteristic (ROC) curves were computed to assess specificity criteria applied in routine doping control. This approach underlined the pertinence to monitor both glucuronide and sulphate steroid conjugates and include them in the athletes passport, while promising biomarkers were also highlighted.

Biomarqueurs associés à la dyspnée et douleurs thoraciques: quelle utilité clinique [Biomarkers associated with dyspnea and chest pain: which cinical utility?].

Dyspnea and chest pain are typical reasons for consultations. biomarkers (CRP, procalcitonin, NT-proBNP, troponins, D-dimers) can have an interest for the diagnosis, the prognosis and the follow-up of several pathologies. There are however numerous pitfalls and limitations between the discovery of a biomarker and the utility in clinical practice. It is essential to always estimate a pre-test probability based on an attentive history and a careful physical examination, to know the intrinsic and extrinsic qualities of a test, and to determine a threshold of care. A biomarker should be used only if it modifies the patient's care and if it brings him a benefit compared to the patient who has no biomarker.

Discovery and molecular characterization of an ambisense densovirus from South American populations of Solenopsis invicta

In an effort to discover viruses as classical biological control agents, a metatranscriptomics/pyrosequencing approach was used to survey native Solenopsis invicta collected exclusively in Argentina. A new virus was discovered with characteristics consistent with the family Parvoviridae, subfamily Densovirinae. The virus, tentatively named Solenopsis invicta densovirus (SiDNV), represents the first DNA virus discovered in ants (Formicidae) and the first densovirus in a hymenopteran insect. The ambisense genome was 5280 nucleotides in length and the termini possessed asymmetrically positioned inverted terminal repeats, formed hairpin loops, and had transcriptional regulatory elements including CAAT and TATA sites. Phylogenetic analysis revealed that SiDNV belongs to a group that includes two other densoviruses found in insects (Acheta domestica densovirus and Planococcus citri densovirus). SiDNV was prevalent in fire ants from Argentina but completely absent in fire ants found in the USA indicating that this virus has potential for biological control of introduced S. invicta.

Anti-60-kDa heat shock protein antibodies in fetal serum: a biomarker for unexplained small for gestational age fetuses.

Introduction: Small for gestational age (SGA) is an important problem affecting 10% of pregnancies and is associated with significant perinatal morbidity. In about 80% of cases, a probable etiology or a major risk factor can be identified. But almost 20% of SGA cases are considered unexplained. The 60-kDa heat shock protein (HSP60) is a highly immunogenic protein whose synthesis is greatly upregulated under nonphysiological conditions. Bacterial and human HSP60 share a high degree of sequence homology, and immunity to conserved epitopes may result in development of autoimmunity following a bacterial infection. We hypothesized that unexplained SGA could be the consequence of immune sensitization to human HSP60. Methods: Unexplained SGA fetuses were identified by ultrasound biometry with normal Doppler velocimetry and with no detectable maternal or fetal abnormalities. Fetal sera were obtained by cordocentesis performed for a karyotype analysis in cases of unexplained SGA (study group) or for screening of Rhesus incompatibility (control group). Fetal sera were tested for HSP60 antigen and for IgG and IgM anti-HSP60 by ELISA as well as for other immune and hematological parameters. Results: Maternal parameters were similar between the 12 study cases and the 23 control cases. The mean gestational age at cordocentesis was 29 weeks. IgM anti-HSP60 was detected in 12 cases (100%) and in no controls (p < 0.00017), while IgG anti-HSP60 was detected in 7 cases (58%) and only 1 control (p < 0.001). Three of the 4 cases with the highest IgM antibody levels died. There were no differences in fetal serum levels of HSP60 antigen or other immune and hematological markers between the two groups. Conclusion: Fetuses with unexplained SGA are positive for IgM and IgG antibody to human HSP60 and the specific IgM antibody level is predictive of fetal mortality. Detection of these antibodies indicates that a placental perturbation and a fetal autoimmune reaction to HSP60 are associated with this developmental delay.

Copeptin as a diagnostic biomarker for sepsis-related deaths.

Copeptin has been shown to increase in patients with sepsis, severe sepsis and septic shock. In the study herein described, copeptin was measured in a series of sepsis-related fatalities and control cases that underwent medico-legal investigations. No age-dependent differences in copeptin levels in either sepsis or control cases were observed. No correlation between copeptin concentrations and postmortem interval was identified in either group. Copeptin levels were significantly higher in sepsis cases. Moreover, copeptin concentrations in septic cases correlated with procalcitonin, C-reactive protein and interleukin 6 values. These preliminary findings seem to indicate that copeptin can be reliably measured in biological samples collected during postmortem investigations. These results also suggest that hemodynamic instability associated with sepsis and septic shock can be characterized by copeptin measurement also in the forensic casework.

A process for mining science & technology documents databases, illustrated for the case of "knowledge discovery and data mining"

This paper presents a process of mining research & development abstract databases to profile current status and to project potential developments for target technologies, The process is called "technology opportunities analysis." This article steps through the process using a sample data set of abstracts from the INSPEC database on the topic o "knowledge discovery and data mining." The paper offers a set of specific indicators suitable for mining such databases to understand innovation prospects. In illustrating the uses of such indicators, it offers some insights into the status of knowledge discovery research*.

Virtual screening and computational optimization for the discovery of covalent prolyl oligopeptidase inhibitors with activity in human cells.

Our docking program, Fitted, implemented in our computational platform, Forecaster, has been modified to carry out automated virtual screening of covalent inhibitors. With this modified version of the program, virtual screening and further docking-based optimization of a selected hit led to the identification of potential covalent reversible inhibitors of prolyl oligopeptidase activity. After visual inspection, a virtual hit molecule together with four analogues were selected for synthesis and made in one-five chemical steps. Biological evaluations on recombinant POP and FAPα enzymes, cell extracts, and living cells demonstrated high potency and selectivity for POP over FAPα and DPPIV. Three compounds even exhibited high nanomolar inhibitory activities in intact living human cells and acceptable metabolic stability. This small set of molecules also demonstrated that covalent binding and/or geometrical constraints to the ligand/protein complex may lead to an increase in bioactivity.

Methylation of the hTERT Promoter: A Novel Cancer Biomarker for Leptomeningeal Metastasis Detection in Cerebrospinal Fluids.

PURPOSE: The diagnosis of leptomeningeal metastases is usually confirmed by the finding of malignant cells by cytologic examination in the cerebrospinal fluid (CSF). More sensitive and specific cancer biomarkers may improve the detection of tumor cells in the CSF. Promoter methylation of the human telomerase reverse transcriptase (hTERT) gene characterizes most cancer cells. The aim of this study was to develop a sensitive method to detect hTERT methylation and to explore its use as a cancer biomarker in CSF. EXPERIMENTAL DESIGN: In 77 CSF specimens from 67 patients, hTERT promoter methylation was evaluated using real-time methylation-sensitive high-resolution melting (MS-HRM) and real-time TaqMan PCR and MS-HRM in a single-tube assay. RESULTS: Real-time MS-HRM assay was able to detect down to 1% hTERT-methylated DNA in a background of unmethylated DNA. PCR products were obtained from 90% (69/77) of CSF samples. No false positive hTERT was detected in the 21 non-neoplastic control cases, given to the method a specificity of 100%. The sensitivity of the real-time MS-HRM compared with the cytologic gold standard analysis was of 92% (11/12). Twenty-six CSFs from 22 patients with an hTERT-methylated primary tumor showed cytologic results suspicious for malignancy; in 17 (65%) of them, a diagnosis of leptomeningeal metastases could be confirmed by the hTERT methylation test. CONCLUSION: The hTERT real-time MS-HRM approach is fast, specific, sensitive, and could therefore become a valuable tool for diagnosis of leptomeningeal metastases as an adjunct to the traditional examination of CSF. Clin Cancer Res; 19(8); 2216-23. ©2013 AACR.

Gas-chromatography mass-spectrometry determination of phthalic acid in human urine as a biomarker of folpet exposure.

Agricultural workers are exposed to folpet, but biomonitoring data are limited. Phthalimide (PI), phthalamic acid (PAA), and phthalic acid (PA) are the ring metabolites of this fungicide according to animal studies, but they have not yet been measured in human urine as metabolites of folpet, only PA as a metabolite of phthalates. The objective of this study was thus to develop a reliable gas chromatography-tandem mass spectrometry (GC-MS) method to quantify the sum of PI, PAA, and PA ring-metabolites of folpet in human urine. Briefly, the method consisted of adding p-methylhippuric acid as an internal standard, performing an acid hydrolysis at 100 °C to convert ring-metabolites into PA, purifying samples by ethyl acetate extraction, and derivatizing with N,O-bis(trimethylsilyl)trifluoro acetamide prior to GC-MS analysis. The method had a detection limit of 60.2 nmol/L (10 ng/mL); it was found to be accurate (mean recovery, 97%), precise (inter- and intra-day percentage relative standard deviations <13%), and with a good linearity (R (2) > 0.98). Validation was conducted using unexposed peoples urine spiked at concentrations ranging from 4.0 to 16.1 μmol/L, along with urine samples of volunteers dosed with folpet, and of exposed workers. The method proved to be (1) suitable and accurate to determine the kinetic profile of PA equivalents in the urine of volunteers orally and dermally administered folpet and (2) relevant for the biomonitoring of exposure in workers.