MOTIVATION: Pathway modeling requires the integration of multiple data including prior knowledge. In this study, we quantitatively assess the application of Gene Ontology (GO)-derived similarity measures for the characterization of direct and indirect interactions within human regulatory pathways. The characterization would help the integration of prior pathway knowledge for the modeling. RESULTS: Our analysis indicates information content-based measures outperform graph structure-based measures for stratifying protein interactions. Measures in terms of GO biological process and molecular function annotations can be used alone or together for the validation of protein interactions involved in the pathways. However, GO cellular component-derived measures may not have the ability to separate true positives from noise. Furthermore, we demonstrate that the functional similarity of proteins within known regulatory pathways decays rapidly as the path length between two proteins increases. Several logistic regression models are built to estimate the confidence of both direct and indirect interactions within a pathway, which may be used to score putative pathways inferred from a scaffold of molecular interactions.

MOTIVATION: In eukaryotic genomes, histone acetylation and thereafter departure from the chromatin are essential for gene transcription initiation. Because gene transcription is tightly regulated by transcription factors, there are some speculations on the cooperation of histone acetylation and transcription factor binding. However, systematic statistical analyses of this relationship on a genomic scale have not been reported. RESULTS: We apply several statistical methods to explore this relationship on two recent genomic datasets: acetylation levels on 11 histone lysines and binding activities of 203 transcription factors, both in promoter regions across the yeast genome. By canonical correlation analysis, we find that a histone acetylation pattern is correlated with a certain profile of transcription factor binding in the genome. Furthermore, after clustering the genes by their acetylation levels on the 11 histone lysines, the genes within clusters show distinct transcription factor binding profiles, as discovered by principle component analysis. Even after applying fairly stringent statistical measurement, most of these clusters have transcription factors with binding activities significantly deviated from the overall genome. We conclude that in the yeast genome, there is a strong correlation between histone acetylation and transcription factor binding in the promoter regions.

MOTIVATION: Current methods for identifying sequence specific binding sites in DNA sequence using position specific weight matrices are limited in both sensitivity and specificity. Double strand DNA helix exhibits sequence dependent variations in conformation. Interactions between macromolecules result from complementarity of the two tertiary structures. We hypothesize that this conformational variation plays a role in transcription factor binding site recognition, and that the use of this structure information will improve the predictive power of transcription factor binding site models. RESULTS: Conformation models for the sequence dependence of DNA helix distortion have been developed. Using our conformational models, we defined a tertiary structure template for the met operon repressor MetJ binding site. Both naturally occurring sites and precursor binding sites identified through in vitro selection were used as the basis for template definition. The conformational model appears to recognize features of protein binding sites that are distinct from the features recognized by primary sequence based profiles. Combining the conformational model and primary sequence profile yields a hybrid model with improved discriminatory power compared with either the conformational model or sequence profile alone. Using our hybrid model, we searched the E.coli genome. We are able to identify the documented MetJ sites in the promoter regions of metA, metB, metC, metR and metF. In addition, we find several novel loci with characteristics suggesting that they are functional MetJ repressor binding sites. Novel MetJ binding sites are found upstream of the metK gene, as well as upstream of a gene, abc, a gene that encodes for a component of a multifunction transporter which may transport amino acids across the membrane. The false positive rate is significantly lower than the sequence profile method. AVAILABILITY: The programs of implementation of this algorithm are available upon request. The list of crystal structures used for compiling the mean base step parameters of DNA is available by anonymous ftp at http://stateslab.wustl.edu/pub/helix/StructureList.

Double strand DNA exhibits sequence dependent variations in three-dimensional conformation.Crystallographic studies have shown that regulatory proteins bind to DNA and may use these structural effects to indirectly recognize binding sites in a sequence specific manner.In this paper we describe a sequence similarity scoring measure based on the energetics of helix distortion.A set of preferred conformations is compiled from data in the Nucleic Acid Database ( NDB ) for 16 dinucleotide steps.These conformational parameters are combined with elastic force constants derived from solution studies of DNA to compute a log likelihood score for one sequence adopting the preferred conformation of another.