Context-aware visual exploration of molecular databases

Facilitating the visual exploration of scientific data has received increasing attention in the past decade or so. Especially in life science related application areas the amount of available data has grown at a breath taking pace. In this paper we describe an approach that allows for visual inspection of large collections of molecular compounds. In contrast to classical visualizations of such spaces we incorporate a specific focus of analysis, for example the outcome of a biological experiment such as high throughout screening results. The presented method uses this experimental data to select molecular fragments of the underlying molecules that have interesting properties and uses the resulting space to generate a two dimensional map based on a singular value decomposition algorithm and a self organizing map. Experiments on real datasets show that the resulting visual landscape groups molecules of similar chemical properties in densely connected regions.

Allergen databases

Allergies represent a significant medical and industrial problem. Molecular and clinical data on allergens are growing exponentially and in this article we have reviewed nine specialized allergen databases and identified data sources related to protein allergens contained in general purpose molecular databases. An analysis of allergens contained in public databases indicates a high level of redundancy of entries and a relatively low coverage of allergens by individual databases. From this analysis we identify current database needs for allergy research and, in particular, highlight the need for a centralized reference allergen database.

Homology modeling and metabolism prediction of human carboxylesterase-2 using docking analyses by GriDock: a parallelized tool based on AutoDock 4.0.

Metabolic problems lead to numerous failures during clinical trials, and much effort is now devoted to developing in silico models predicting metabolic stability and metabolites. Such models are well known for cytochromes P450 and some transferases, whereas less has been done to predict the activity of human hydrolases. The present study was undertaken to develop a computational approach able to predict the hydrolysis of novel esters by human carboxylesterase hCES2. The study involved first a homology modeling of the hCES2 protein based on the model of hCES1 since the two proteins share a high degree of homology (congruent with 73%). A set of 40 known substrates of hCES2 was taken from the literature; the ligands were docked in both their neutral and ionized forms using GriDock, a parallel tool based on the AutoDock4.0 engine which can perform efficient and easy virtual screening analyses of large molecular databases exploiting multi-core architectures. Useful statistical models (e.g., r (2) = 0.91 for substrates in their unprotonated state) were calculated by correlating experimental pK(m) values with distance between the carbon atom of the substrate's ester group and the hydroxy function of Ser228. Additional parameters in the equations accounted for hydrophobic and electrostatic interactions between substrates and contributing residues. The negatively charged residues in the hCES2 cavity explained the preference of the enzyme for neutral substrates and, more generally, suggested that ligands which interact too strongly by ionic bonds (e.g., ACE inhibitors) cannot be good CES2 substrates because they are trapped in the cavity in unproductive modes and behave as inhibitors. The effects of protonation on substrate recognition and the contrasting behavior of substrates and products were finally investigated by MD simulations of some CES2 complexes.

Characterization of a new clinical yeast species, Candida tunisiensis sp. nov., isolated from a strain collection from Tunisian hospitals.

From a collection of yeast isolates isolated from patients in Tunisian hospitals between September 2006 and July 2010, the yeast strain JEY63 (CBS 12513), isolated from a 50-year-old male that suffered from oral thrush, could not be identified to the species level using conventional methods used in clinical laboratories. These methods include matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS), germ tube formation, and the use of CHROMagar Candida and metabolic galleries. Sequence analysis of the nuclear rRNA (18S rRNA, 5.8S rRNA, and 26S rRNA) and internal transcribed spacer regions (ITS1 and ITS2) indicated that the ribosomal DNA sequences of this species were not yet reported. Multiple gene phylogenic analyses suggested that this isolate clustered at the base of the Dipodascaceae (Saccharomycetales, Saccharomycetes, and Ascomycota). JEY63 was named Candida tunisiensis sp. nov. according to several phenotypic criteria and its geographical origin. C. tunisiensis was able to grow at 42°C and does not form chlamydospores and hyphae but could grow as yeast and pseudohyphal forms. C. tunisiensis exhibited most probably a haploid genome with an estimated size of 10 Mb on at least three chromosomes. Using European Committee for Antimicrobial Susceptibility Testing (EUCAST) and Clinical and Laboratory Standards Institute (CLSI) Candida albicans susceptibility breakpoints as a reference, C. tunisiensis was resistant to fluconazole (MIC = 8 μg/ml), voriconazole (MIC = 0.5 μg/ml), itraconazole (MIC = 16 μg/ml), and amphotericin B (MIC = 4 μg/ml) but still susceptible to posaconazole (MIC = 0.008 μg/ml) and caspofungin (MIC = 0.5 μg/ml). In conclusion, MALDI-TOF MS permitted the early selection of an unusual isolate, which was still unreported in molecular databases but could not be unambiguously classified based on phylogenetic approaches.

MGMT methylation analysis of glioblastoma on the Infinium methylation BeadChip identifies two distinct CpG regions associated with gene silencing and outcome, yielding a prediction model for comparisons across datasets, tumor grades, and CIMP-status.

The methylation status of the O(6)-methylguanine-DNA methyltransferase (MGMT) gene is an important predictive biomarker for benefit from alkylating agent therapy in glioblastoma. Recent studies in anaplastic glioma suggest a prognostic value for MGMT methylation. Investigation of pathogenetic and epigenetic features of this intriguingly distinct behavior requires accurate MGMT classification to assess high throughput molecular databases. Promoter methylation-mediated gene silencing is strongly dependent on the location of the methylated CpGs, complicating classification. Using the HumanMethylation450 (HM-450K) BeadChip interrogating 176 CpGs annotated for the MGMT gene, with 14 located in the promoter, two distinct regions in the CpG island of the promoter were identified with high importance for gene silencing and outcome prediction. A logistic regression model (MGMT-STP27) comprising probes cg1243587 and cg12981137 provided good classification properties and prognostic value (kappa = 0.85; log-rank p < 0.001) using a training-set of 63 glioblastomas from homogenously treated patients, for whom MGMT methylation was previously shown to be predictive for outcome based on classification by methylation-specific PCR. MGMT-STP27 was successfully validated in an independent cohort of chemo-radiotherapy-treated glioblastoma patients (n = 50; kappa = 0.88; outcome, log-rank p < 0.001). Lower prevalence of MGMT methylation among CpG island methylator phenotype (CIMP) positive tumors was found in glioblastomas from The Cancer Genome Atlas than in low grade and anaplastic glioma cohorts, while in CIMP-negative gliomas MGMT was classified as methylated in approximately 50 % regardless of tumor grade. The proposed MGMT-STP27 prediction model allows mining of datasets derived on the HM-450K or HM-27K BeadChip to explore effects of distinct epigenetic context of MGMT methylation suspected to modulate treatment resistance in different tumor types.

Hands-On: Interlinking Institutional Repositories and Europe PubMed Central

Presentation at Open Repositories 2014, Helsinki, Finland, June 9-13, 2014

Deduction of functional peptide motifs in scorpion toxins

Scorpion toxins are important physiological probes for characterizing ion channels. Molecular databases have limited functional annotation of scorpion toxins. Their function can be inferred by searching for conserved motifs in sequence signature databases that are derived statistically but are not necessarily biologically relevant. Mutation studies provide biological information on residues and positions important for structure-function relationship but are not normally used for extraction of binding motifs. 3D structure analyses also aid in the extraction of peptide motifs in which non-contiguous residues are clustered spatially. Here we present new, functionally relevant peptide motifs for ion channels, derived from the analyses of scorpion toxin native and mutant peptides. Copyright (c) 2006 European Peptide Society and John Wiley & Sons, Ltd.

Molecular immunology databases and data repositories

Over recent years databases have become an extremely important resource for biomedical research. Immunology research is increasingly dependent on access to extensive biological databases to extract existing information, plan experiments, and analyse experimental results. This review describes 15 immunological databases that have appeared over the last 30 years. In addition, important issues regarding database design and the potential for misuse of information contained within these databases are discussed. Access pointers are provided for the major immunological databases and also for a number of other immunological resources accessible over the World Wide Web (WWW). (C) 2000 Elsevier Science B.V. All rights reserved.

Molecular modeling databases: A new way in the search of protein targets for drug development

DBMODELING is a relational database of annotated comparative protein structure models and their metabolic, pathway characterization. It is focused on enzymes identified in the genomes of Mycobacterium tuberculosis and Xylella fastidiosa. The main goal of the present database is to provide structural models to be used in docking simulations and drug design. However, since the accuracy of structural models is highly dependent on sequence identity between template and target, it is necessary to make clear to the user that only models which show high structural quality should be used in such efforts. Molecular modeling of these genomes generated a database, in which all structural models were built using alignments presenting more than 30% of sequence identity, generating models with medium and high accuracy. All models in the database are publicly accessible at http://www.biocristalografia.df.ibilce.unesp.br/tools. DBMODELING user interface provides users friendly menus, so that all information can be printed in one stop from any web browser. Furthermore, DBMODELING also provides a docking interface, which allows the user to carry out geometric docking simulation, against the molecular models available in the database. There are three other important homology model databases: MODBASE, SWISSMODEL, and GTOP. The main applications of these databases are described in the present article. © 2007 Bentham Science Publishers Ltd.

Chimeric 16S rDNA sequences of diverse origin are accumulating in the public databases

A significant number of chimeric 16S rDNA sequences of diverse origin were identified in the public databases by partial treeing analysis. This suggests that chimeric sequences, representing phylogenetically novel non-existent organisms, are routinely being overlooked in molecular phylogenetic surveys despite a general awareness of PCR-generated artefacts amongst researchers.

Protein interaction data curation: the International Molecular Exchange (IMEx) consortium.

The International Molecular Exchange (IMEx) consortium is an international collaboration between major public interaction data providers to share literature-curation efforts and make a nonredundant set of protein interactions available in a single search interface on a common website (http://www.imexconsortium.org/). Common curation rules have been developed, and a central registry is used to manage the selection of articles to enter into the dataset. We discuss the advantages of such a service to the user, our quality-control measures and our data-distribution practices.

Obtenció de nous anàlegs amb activitat brassinoesteroide mitjançant modelització molecular i síntesi

Els brassinoesteroides són productes naturals que actuen com a potents reguladors del creixement vegetal. Presenten aplicacions prometedores en l’agricultura degut a que, aplicats exògenament, augmenten la qualitat i la quantitat de les collites. Ara bé, el seu ús s’ha vist restringit degut a la seva costosa obtenció. Aquest fet ha motivat la recerca de nous compostos actius més assequibles. En aquest projecte es planteja el disseny i obtenció de nous anàlegs seguint diferents estratègies que impliquen tant l’ús de mètodes de modelització molecular com de síntesi orgànica. La primera d’aquestes estratègies consisteix en buscar compostos actius en bases de dades de compostos comercials a través de processos de Virtual Screening desenvolupats amb mètodes computacionals basats en Camps d’Interacció Molecular. Així, es van establir i interpretar models de Relacions Quantitatives Estructura-Activitat (QSAR) emprant descriptors independents de l’alineament (GRIND) i, amb col•laboració amb la Universitat de Perugia, aquest criteri de cerca es va ampliar amb l’aplicació de descriptors FLAP de nova generació. Una altra estratègia es va basar en intentar substituir l’esquelet esteroide dels brassinoesteroides per una estructura equivalent, fixant com a cadena lateral el grup (R)-hexahidromandelil. S’han aplicat dos criteris: mètodes computacionals basats en models QSAR establerts amb descriptors GRIND i també en la metodologia SHOP (scaffold hopping), i, per altra banda, anàlegs proposats racionalment a partir d’un estudi efectuat sobre disruptors endocrins no esteroïdals. Sobre les estructures trobades s’hi va unir la cadena lateral comercial esmentada per via sintètica, en la qual s’ha hagut de fer un èmfasi especial en grups protectors. En total, 49 estructures es proposen per a ser obtingudes sintèticament. També s’ha treballat en l’obtenció un agonista derivat de l’hipotètic antagonista KM-01. Totes les molècules candidates, ja siguin comercials o obtingudes sintèticament, estant sent avaluades en el test d’inclinació de la làmina d’arròs (RLIT).

trEST, trGEN and Hits: access to databases of predicted protein sequences.

High throughput genome (HTG) and expressed sequence tag (EST) sequences are currently the most abundant nucleotide sequence classes in the public database. The large volume, high degree of fragmentation and lack of gene structure annotations prevent efficient and effective searches of HTG and EST data for protein sequence homologies by standard search methods. Here, we briefly describe three newly developed resources that should make discovery of interesting genes in these sequence classes easier in the future, especially to biologists not having access to a powerful local bioinformatics environment. trEST and trGEN are regularly regenerated databases of hypothetical protein sequences predicted from EST and HTG sequences, respectively. Hits is a web-based data retrieval and analysis system providing access to precomputed matches between protein sequences (including sequences from trEST and trGEN) and patterns and profiles from Prosite and Pfam. The three resources can be accessed via the Hits home page (http://hits. isb-sib.ch).

Molecular evolutionary patterns in the Glomeromycota;

Abstract Arbuscular mycorrhizal fungi (AMF) form symbiosis with roots of approximately 80% of known land plants. These fungi play a key role in the ecology and adaptation of plants to various ecosystems.by increasing the plant resources for various nutrients. Despite their important ecological role, we still have poor understanding of their genetic structure and their molecular evolution. The work presented in this thesis aims to isolate and analyse AMF genes with various molecular techniques, in order to obtain new insights about their genetics, phylogeny and molecular evolution. Some AMF genes were shown through phylogenetic analyses to be more related with plants or mycoparasites than with other fungal organisms. These results led to the prediction that lateral gene transfers (LGT) occurred between AMF and plants during their long-term co-évolution. By phylogenetic and molecular analyses, in the chapter 2 I demonstrate that the hypothesis of LGT is most likely a consequence of analyses carried out on contaminant non AMF-DNA. In addition, various features characteristic of AMF genes have been determined, allowing researchers to scan their own sequence databases for potential non-AMF contaminants. Phylogenetic relationships of AMF with other fungi has been mostly analysed using molecular markers of ribosomal origin. In chapter 2 I successfully isolated gene encoding α- and ß-tubulins from several AMF genera. Consequently, phylogenetic analyses showed that AMF possess an unexpected relationship with ancestral aquatic fungi (chytrids). These results are consistent with the prediction stating that AMF may have played an important role in the colonisation of land by green plants through the establishment of a symbiosis and after the divergence of AMF from aquatic ancestors. In Chapter 4 I tried to isolate the entire AMF gene family encoding P-Type II ATPases, in order to determine their molecular evolution with the fungal kingdom. These genes were further analysed to detect the level of sequence polymorphism that is present within an AMF population. The results obtained show that mutational events previously thought as occurring only among divergent evolutionary lineages (gene duplications, indel mutations in coding regions) can occur within a single population of AMF. These results have far reaching consequences for our understanding of the genetics and ecology of AMF. Résumé Les champignons endomycorrhiziens arbusculaires (CEA) forment une symbiose racinaire avec environ 80% des plantes vasculaires connues. Ces champignons possèdent un rôle important dans l'écologie et l'adaptation des plantes au sein de différents écosystèmes en .augmentant leurs ressources en nutriments. Le travail présenté dans cette thèse se propose d'isoler et d'analyser certains gènes de CEA avec différentes techniques moléculaires à fin d'obtenir de pÌus amples informations concernant l'évolution moléculaire, la phylogénie et leur diversité génétique à diverses échelles taxonomiques. Certaines analyses phylogénétiques des CEA ont conduit à l'hypothèse que des transferts horizontaux de gènes (THG) ont pu avoir lieu durant leur longue co-évolution avec les plantes vasculaires. Dans le chapitre 2 de cette thèse nous démontrons par analyses moléculaire et phylogénétique que l'hypothèse de THG est une conséquence de contaminations à partir d'ADN de plante ou d'autres micro-organismes. De plus, de nombreuses caractéristiques moléculaires de CEA ont pu être déterminées, permettant la mise en place d'un plan à suivre lors de l'analyse de gènes de CEA dans les études futures. Les relations évolutives des. CEA avec d'autres champignons ont été analysées majoritairement à l'aide de marqueurs moléculaires d'origine ribosomiale. Dans les chapitres 2 et 3 j'ai isolé des gènes codant pour l'a- et la ß-tubuline chez différents genres, de CEA. Les analyses phylogénétiques ont démontré une parenté entre les CEA et des champignons aquatiques ancestraux (chytrides). Ces résultats sont en accord avec l'hypothèse selon laquelle les CEA ont probablement joué un rôle primordial dans l'établissement des plantes sur terre à travers une symbiose et suite à leur évolution à partir d'ancêtres vivant dans des milieux aquatiques: Dans le chapitre 4 j'ai isolé une entière famille de gènes chez les CEA codant des ATPases de la membrane plasmique, et étudié leur évolution moléculaire dans le règne des champignons. Ces mêmes gènes ont été analysés ultérieurement à fin de déterminer le degré de polymorphisme de séquence qui peut être présent au sein d'une population de CEA. Les résultats obtenus montrent que des évènements mutationnels considérés comme apparaissant exclusivement dans des lignées évolutives très divergentes (duplication de gènes, insertions/délétions dans des régions transcrites du génome) ont lieu sein d'une même population de CEA. Cette découverte a un impact important sur nos connaissances concernant la génétique des populations des CEA et leur écologie.

Use of the frc gene as a molecular marker to characterize oxalate-oxidizing bacterial abundance and diversity structure in soil.

Oxalate catabolism, which can have both medical and environmental implications, is performed by phylogenetically diverse bacteria. The formyl-CoA-transferase gene was chosen as a molecular marker of the oxalotrophic function. Degenerated primers were deduced from an alignment of frc gene sequences available in databases. The specificity of primers was tested on a variety of frc-containing and frc-lacking bacteria. The frc-primers were then used to develop PCR-DGGE and real-time SybrGreen PCR assays in soils containing various amounts of oxalate. Some PCR products from pure cultures and from soil samples were cloned and sequenced. Data were used to generate a phylogenetic tree showing that environmental PCR products belonged to the target physiological group. The extent of diversity visualised on DGGE pattern was higher for soil samples containing carbonate resulting from oxalate catabolism. Moreover, the amount of frc gene copies in the investigated soils was detected in the range of 1.64x10(7) to 1.75x10(8)/g of dry soil under oxalogenic tree (representing 0.5 to 1.2% of total 16S rRNA gene copies), whereas the number of frc gene copies in the reference soil was 6.4x10(6) (or 0.2% of 16S rRNA gene copies). This indicates that oxalotrophic bacteria are numerous and widespread in soils and that a relationship exists between the presence of the oxalogenic trees Milicia excelsa and Afzelia africana and the relative abundance of oxalotrophic guilds in the total bacterial communities. This is obviously related to the accomplishment of the oxalate-carbonate pathway, which explains the alkalinization and calcium carbonate accumulation occurring below these trees in an otherwise acidic soil. The molecular tools developed in this study will allow in-depth understanding of the functional implication of these bacteria on carbonate accumulation as a way of atmospheric CO(2) sequestration.