Evidence of distinct gene functional patterns in GC-poor and GC-rich isochores in Bos taurus.

Vertebrate genomes are mosaics of megabase-size DNA segments with a fairly homogeneous base composition, called isochores. They are divided into five families characterized by different guanine-cytosine (GC) levels and linked to several functional and structural properties. The increased availability of fully sequenced genomes allows the investigation of isochores in several species, assessing their level of conservation across vertebrate genomes. In this work, we characterized the isochores in Bos taurus using the ARS-UCD1.2 genome version. The comparison of our results with the well-studied human isochores and those of other mammals revealed a large conservation in isochore families, in number, average GC levels and gene density. Exceptions to the established increase in gene density with the increase in isochores (GC%) were observed for the following gene biotypes: tRNA, small nuclear RNA, small nucleolar RNA and pseudogenes that have their maximum number in H2 and H1 isochores. Subsequently, we assessed the ontology of all gene biotypes looking for functional classes that are statistically over- or under-represented in each isochore. Receptor activity and sensory perception pathways were significantly over-represented in L1 and L2 (GC-poor) isochores. This was also validated for the horse genome. Our analysis of housekeeping genes confirmed a preferential localization in GC-rich isochores, as reported in other species. Finally, we assessed the SNP distribution of a bovine high-density SNP chip across the isochores, finding a higher density in the GC-rich families, reflecting a potential bias in the chip, widely used for genetic selection and biodiversity studies.

Network-guided modeling allows tumor-type independent prediction of sensitivity to all-trans-retinoic acid.

Background: All-trans-retinoic acid (ATRA) is a differentiating agent used in the treatment of acute-promyelocytic-leukemia (APL) and it is under-exploited in other malignancies despite its low systemic toxicity. A rational/personalized use of ATRA requires the development of predictive tools allowing identification of sensitive cancer types and responsive individuals. Materials and methods: RNA-sequencing data for 10 080 patients and 33 different tumor types were derived from the TCGA and Leucegene datasets and completely re-processed. The study was carried out using machine learning methods and network analysis. Results: We profiled a large panel of breast-cancer cell-lines for in vitro sensitivity to ATRA and exploited the associated basal gene-expression data to initially generate a model predicting ATRA-sensitivity in this disease. Starting from these results and using a network-guided approach, we developed a generalized model (ATRA-21) whose validity extends to tumor types other than breast cancer. ATRA-21 predictions correlate with experimentally determined sensitivity in a large panel of cell-lines representative of numerous tumor types. In patients, ATRA-21 correctly identifies APL as the most sensitive acute-myelogenous-leukemia subtype and indicates that uveal-melanoma and low-grade glioma are top-ranking diseases as for average predicted responsiveness to ATRA. There is a consistent number of tumor types for which higher ATRA-21 predictions are associated with better outcomes. Conclusions: In summary, we generated a tumor-type independent ATRA-sensitivity predictor which consists of a restricted number of genes and has the potential to be applied in the clinics. Identification of the tumor types that are likely to be generally sensitive to the action of ATRA paves the way to the design of clinical studies in the context of these diseases. In addition, ATRA-21 may represent an important diagnostic tool for the selection of individual patients who may benefit from ATRA-based therapeutic strategies also in tumors characterized by lower average sensitivity.

ASPicDB: a database resource for alternative splicing analysis.

MOTIVATION: Alternative splicing has recently emerged as a key mechanism responsible for the expansion of transcriptome and proteome complexity in human and other organisms. Although several online resources devoted to alternative splicing analysis are available they may suffer from limitations related both to the computational methodologies adopted and to the extent of the annotations they provide that prevent the full exploitation of the available data. Furthermore, current resources provide limited query and download facilities. RESULTS: ASPicDB is a database designed to provide access to reliable annotations of the alternative splicing pattern of human genes and to the functional annotation of predicted splicing isoforms. Splice-site detection and full-length transcript modeling have been carried out by a genome-wide application of the ASPic algorithm, based on the multiple alignments of gene-related transcripts (typically a Unigene cluster) to the genomic sequence, a strategy that greatly improves prediction accuracy compared to methods based on independent and progressive alignments. Enhanced query and download facilities for annotations and sequences allow users to select and extract specific sets of data related to genes, transcripts and introns fulfilling a combination of user-defined criteria. Several tabular and graphical views of the results are presented, providing a comprehensive assessment of the functional implication of alternative splicing in the gene set under investigation. ASPicDB, which is regularly updated on a monthly basis, also includes information on tissue-specific splicing patterns of normal and cancer cells, based on available EST sequences and their library source annotation. AVAILABILITY: www.caspur.it/ASPicDB

The mepsMAP server. Mapping epitopes on protein surface: mining annotated proteins.

For a growing number of biologists DNA or protein data are typically retrieved and managed on the Web, and not in the laboratory. A large number of bioinformatics datasets from primary and (thousands of) secondary databases are scattered on the Web in various formats. A biologist end-user might need to access and use tens of databases and tools every day. For this reason, the bioinformatics community is developing more and more service-oriented architectures (SOAs): software architecture of loosely coupled software services that can be accessed without knowledge of, or control over, their internal architecture. Data-processing and analysis tasks can be automated by having free access to bioinformatics Web services (WSs) that are the building blocks of the SOAs. In this paper we introduce a new bioinformatics Web server, mepsMAP (mapping epitopes on protein surface: Mining Annotated Proteins), developed to identify the recognition sites between antibodies and their cognate antigens. In some cases, the recognition site is represented by a continuous segment of the antigen sequence, but much more often the epitope is "conformational," i.e., the antibody recognizes the location and type of exposed antigen side chains that are not necessarily contiguous in the antigen's sequence, but brought together by its three-dimensional structure. A facility on the server allows the user to search putative conformational epitopes on protein surface, querying the system for proteins with a given annotation. The mepsMAP server has been implemented as a SOA composed by a database and a set of four WSs. We present here the software architecture of the system with a detailed description of the WS dataflow that has been optimized to provide the best computing performance while maintaining the easiest end-user access to the system via a Web interface.

Structure and hydration of the DNA-human topoisomerase I covalent complex.

The structure and hydration of reconstituted human topoisomerase I comprising the core and the carboxyl-terminal domains in covalent complex with 22-basepair DNA duplex has been investigated by molecular dynamics simulation. The structure and the intermolecular interactions were found to be well maintained over the simulation. The complex displays a high degree of flexibility of the contact area, confirmed by the presence of numerous water-mediated protein-DNA hydrogen bonds comparable in quantity and distribution to the direct ones. The interaction between the enzyme and the solvent also provides the key for interpreting the experimental reduction of activity or affinity observed upon single residue mutation. Finally, four long lasting water molecules are observed in the proximity of the active site, one of which in the appropriate position to accept a proton from the active Tyr723.

Structure and hydration of BamHI DNA recognition site: a molecular dynamics investigation.

The results of a 3-ns molecular dynamics simulation of the dodecamer duplex d(TATGGATCCATA)(2) recognized by the BamHI endonuclease are presented here. The DNA has been simulated as a flexible molecule using an AMBER force field and the Ewald summation method, which eliminates the undesired effects of truncation and permits evaluation of the full effects of electrostatic forces. The starting B conformation evolves toward a configuration quite close to that observed through x-ray diffraction in its complex with BamHI. This configuration is fairly stable and the Watson-Crick hydrogen bonds are well maintained over the simulation trajectory. Hydration analysis indicates a preferential hydration for the phosphate rather than for the ester oxygens. Hydration shells in both the major and minor groove were observed. In both grooves the C-G pairs were found to be more hydrated than A-T pairs. The "spine of hydration" in the minor groove was clear. Water residence times are longer in the minor groove than in the major groove, although relatively short in both cases. No special long values are observed for sites where water molecules were observed by x-ray diffraction, indicating that water molecules having a high probability of being located in a specific site are also fast-exchanging.

A novel parameter to estimate the minimum number of bound ligands needed to activate an ion channel.

We consider the Hill coefficient obtained from concentration-response curves and we stress the differences between the use of binding and the use of current for the response. We first show that in order to estimate the cooperativity of binding sites one must use the coefficient obtained from the binding curves and not that obtained from the current curves. For what concerns the estimation of the minimum number of bound ligands needed to activate the channel, while the coefficient obtained from the binding is simply not applicable, we show that the coefficient obtained from the current is not appropriate. We finally introduce a novel parameter that predicts the number of ligand molecules required to open the channel from current measurements at very low ligand concentration. The above considerations are exemplified by a few theoretical models.

Are only repeated triplets guilty?

It is well known that in some places of the human genome one finds a variable number of tandem repeats of trinucleotides; it is now commonly acknowledged that in many cases an excessive expansion of such a number is the cause of nervous system diseases. Moreover there exist cases of genetic disorders linked with loci where a variable number of tandem repeats of sequences longer than three bases has been found. The abnormal number of these repeats in few cases has been associated with the onset of the disease. Considering the above facts, we have performed an extensive study of published sequences of genes connected with various diseases. We have examined, inside or near those genes, all possible tandem repeats. The analysis has led to the detection of a large number of repeats of both triplets and longer sequences, many of which, as far as we know, had not been pointed out before. The results of our analysis lead us to put forward the hypothesis that in more cases than those till now established, a variable number of tandem repeats of generic sequences, not only of triplets, could be associated with disease onset. Finally we suggest to allocate experimental researches for all the possible tandem repeats and their possible correlation with the neurodegenerative disorders and with other kinds of syndromes.

A study of oligonucleotide occurrence distributions in DNA coding segments.

In this paper we present a general strategy designed to study the occurrence frequency distributions of oligonucleotides in DNA coding segments and to deal with the problem of detecting possible patterns of genomic compositional inhomogeneities and disuniformities. Identifying specific tendencies or peculiar deviations in the distributions of the effective occurrence frequencies of oligonucleotides, with respect to what can be a priori expected, is of the greatest importance in biology. Differences between expected and actual distributions may in fact suggest or confirm the existence of specific biological mechanisms related to them. Similarly, a marked deviation in the occurrence frequency of an oligonucleotide may suggest that it belongs to the class of so-called "DNA signal (target) sequences". The approach we have elaborated is innovative in various aspects. Firstly, the analysis of the genomic data is carried out in the light of the observation that the distribution of the four nucleotides along the coding regions of the genoma is biased by the existence of a well-defined "reading frame". Secondly, the "experimental" numbers found by counting the occurrences of the various oligonucleotide sequences are appropriately corrected for the many kinds of mistakes and redundancies present in the available genetic Data Bases. A methodologically significant further improvement of our approach over the existing searching strategies is represented by the fact that, in order to decide whether or not the (corrected) "experimental" value of the occurrence frequency of a given oligonucleotide is within statistical expectations, a measure of the strength of the selective pressure, having acted on it in the course of the evolution, is assigned to the sequence, in a way that takes into account both the value of the "experimental" occurrence frequency of the sequence and the magnitude of the probability that this number might be the result of statistical fluctuations. If the strength of the selective pressure evaluated in this fashion turns out to be sufficiently large, the corresponding sequence will be considered to have an occurrence frequency beyond expectations and, hence, to be statistically and biologically interesting.