Our group produced the best predictions overall in the DREAM3 signaling response challenge, being... more Our group produced the best predictions overall in the DREAM3 signaling response challenge, being tops by a substantial margin in the cytokine sub-challenge and nearly tied for best in the phosphoprotein sub-challenge. We achieved this success using a simple interpolation strategy. For each combination of a stimulus and inhibitor for which predictions were required, we had noted there were six other datasets using the same stimulus (but different inhibitor treatments) and six other datasets using the same inhibitor (but different stimuli). Therefore, for each treatment combination for which values were to be predicted, we calculated rank correlations for the data that were in common between the treatment combination and each of the 12 related combinations. The data from the 12 related combinations were then used to calculate missing values, weighting the contributions from each experiment based on the rank correlation coefficients. The success of this simple method suggests that the missing data were largely over-determined by similarities in the treatments. We offer some thoughts on the current state and future development of DREAM that are based on our success in this challenge, our success in the earlier DREAM2 transcription factor target challenge, and our experience as the data provider for the gene expression challenge in DREAM3.
DNA repair methyltransferases (MTases) remove methyl or other alkyl groups from the O 6 position ... more DNA repair methyltransferases (MTases) remove methyl or other alkyl groups from the O 6 position of guanine or the O 4 position of thymine by transfering the group to an active site cysteine. In order to trap an MTase-DNA complex via a disulfide bond, 2′-deoxy-6-(cystamine)-2-aminopurine (d 6 Cys 2 AP) was synthesized and incorporated into oligonucleotides. d 6 Cys 2 AP has a disulfide bond within an alkyl chain linked to the 6 position of 2,6-diaminopurine, which disulfide can be reduced to form a free thiol. Addition of human MTase to reduced oligonucleotide resulted in a protein-DNA complex that was insensitive to denaturation by SDS and high salt, but which readily dissociated in the presence of dithiothreitol. Formation of this complex was prevented by methylation of the active site cysteine. Evidence that the active site cysteine is directly involved in disulfide bond formation was obtained by N-terminal sequencing of peptides that remained associated with DNA after proteolysis of the complex.
Transcription factors that belong to the same family typically have similar, but not identical, b... more Transcription factors that belong to the same family typically have similar, but not identical, binding specificities. As such, they can be expected to compete differentially for binding to different variants of their binding sites. Pho4 is a yeast factor whose nuclear concentration is up-regulated in low phosphate, while the related factor, Cbf1, is constitutively expressed. We constructed 16 GFPreporter genes containing all palindromic variants of the motif NNCACGTGNN, and determined their activities at a range of phosphate concentrations. Pho4 affinity did not explain expression data well except under fully induced conditions. However, reporter activity was quantitatively well explained under all conditions by a model in which Cbf1 itself has modest activating activity, and Pho4 and Cbf1 compete with one another. Chromatin immunoprecipitation and computational analyses of natural Pho4 target genes, along with the activities of the reporter constructs, indicates that genes differ in their sensitivity to intermediate induction signals in part because of differences in their affinity for Cbf1. The induction sensitivity of both natural Pho4 target genes and reporter genes was well explained only by a model that assumes a role for Cbf1 in remodeling chromatin. Our analyses highlight the importance of taking into account the activities of related transcription factors in explaining system-wide gene expression data.
Several recent papers show that differences in histone modification and the use of histone varian... more Several recent papers show that differences in histone modification and the use of histone variants at the 5′ and 3′ ends of genes influence the location and kinetics of transcriptional initiation. The ultimate target of most epigenetic mechanisms may be the regulation of nucleosome occupancy, which in turn controls access to DNA at specific genomic locations.
Homeodomains are 60 amino acid DNA binding domains found in numerous eukaryotic transcription fac... more Homeodomains are 60 amino acid DNA binding domains found in numerous eukaryotic transcription factors. The homeodomain family is a useful system for studying sequence-structure relationships because several hundred sequences are known and the structures of several homeodomains have been determined. Covariation of amino acid residues in the homeodomain family has been investigated to see whether strongly covariant residue pairs can be understood in terms of the structure and function of these domains. Among 16 strongly covariant pairs examined, 2 are explained by the ability to form salt bridges, and 9 appear related to the DNA binding function of the proteins. For the remaining 5 pairs, the rationale for covariance remains unclear and the likelihood of artifactual correlations is discussed in the context of experimental and evolutionary biases in the selection of sequences. No significant correlation was found between covariance and structural proximity in the hydrophobic core.
Introduction: Opsins are a large and sequence-diverse family of light-responsive G-protein couple... more Introduction: Opsins are a large and sequence-diverse family of light-responsive G-protein coupled receptors involved in vision, circadian rhythm, and other processes. Numerous subfamilies have been defined based on sequence similarity, cell-type localization, signal transduction mechanism, or biological function, but there is no consensus classification system. Methods: We used multiple hidden Markov models (HMMs) to identify opsins in the UniProt Reference Proteomes database. Opsin-specific HMMs were also used in an annotation procedure that represents sequences as a vector of HMM scores and assess the similarity of these vectors to those of annotated sequences. UniProt Reference Proteomes are built from genome sequences, allowing us to make meaningful comparisons of the number of opsins in each of the 260 species available at the time of the survey in absolute terms and relative to a larger superfamily of which opsins are a member. Results: More than 2,000 opsins were retrieved from 262 species (all metazoans). Discussion: Merging opsin counts into higher order taxa paints a broad view of the taxonomic distribution of opsins, and of opsin subfamilies, annotated according to three different schemes.
Hidden Markov models representing 167 protein sequence families were used to infer the presence o... more Hidden Markov models representing 167 protein sequence families were used to infer the presence or absence of homologs within the transcriptomes of 183 algal species/strains. Statistical analyses of the distribution of HMM hits across major clades of algae, or at branch points on the phylogenetic tree of 98 chlorophytes, confirmed and extended known cases of metabolic loss and gain, most notably the loss of the mevalonate pathway for terpenoid synthesis in green algae but not, as we show here, in the streptophyte algae. Evidence for novel events was found as well, most remarkably in the recurrent and coordinated gain or loss of enzymes for the glyoxylate shunt. We find, as well, a curious pattern of retention (or re-gain) of HMG-CoA synthase in chlorophytes that have otherwise lost the mevalonate pathway, suggesting a novel, co-opted function for this enzyme in select lineages. Finally, we find striking, phylogenetically linked distributions of coding sequences for three pathways that synthesize the major membrane lipid phosphatidylcholine, and a complementary phylogenetic distribution pattern for the non-phospholipid DGTS (diacyl-glyceryl-trimethylhomoserine). Mass spectrometric analysis of lipids from 25 species was used to validate the inference of DGTS synthesis from sequence data. Photosynthetic eukaryotes arose from the endosymbiosis of a photosynthetic bacterium by a eukaryotic cell 1 . The descendants of this endosymbiotic event, with the exception of the land plants (Embryophyta), are all referred to as 'algae' . The algae include the Glaucophyta (a small group of freshwater species that have features thought to resemble those of the common ancestor), the Rhodophyta ('red algae' , some of which are unicellular and others, complex seaweeds), and the 'green algae' . The 'green algae' are essentially the Viridiplantae ("green plants") that are not embryophytes (land plants). They are extraordinarily diverse organisms, and by definition consist of lineages that are as old or older than the land plants. They consist of multicellular seaweeds (Ulvophyceae) as well as numerous classes of mostly unicellular organisms, Chlorophyceae and Trebouxiophyceae being the two largest groups. The green algae also include prasinophytes, a paraphyletic group that includes several classes of mostly planktonic marine species. Together, all of these groups are considered to be members of the Chlorophyta. The green algae also include the streptophyte algae, a paraphyletic group, that, together with embryophytes, form the sister group to the Chlorophyta. Most algae are unicellular, but many are not; indeed, multicellularity evolved several times in this group 2 . Given this diversity of algal morphology, and the presence of species in nearly every imaginable ecological niche, it is reasonable to suppose that there have been interesting adaptations in metabolism. The potential for metabolic novelty among the algae is of particular interest because these organisms are poised to be of increasing importance to biotechnology 3 .
Background: Effects on gene expression due to environmental or genetic changes can be easily meas... more Background: Effects on gene expression due to environmental or genetic changes can be easily measured using microarrays. However, indirect effects on expression can be substantial. The indirect effects of a perturbation need to be distinguished from the direct effects if we are to understand the structure and behavior of regulatory networks. The most direct way to perturb a transcriptional network is to alter transcription factor activity. Here, for the first time, we compare expression changes and genomic binding in a simple regulon under conditions of both low and high transcription factor activity. Specifically, we assessed the effects on expression and binding due to deletion of the yeast LEU3 transcription factor gene and effects due to elevation of Leu3 activity. Leu3 activity was elevated through overexpression and the introduction of a mutation that renders the protein constitutively active. Genes that are bound and/or regulated by Leu3 under one or both conditions were characterized in terms of their functional annotations and their predicted potential to be bound by Leu3. We also assessed the evolutionary conservation of the predicted binding potential using a novel alignment-independent method. Both perturbations yield genes that are likely to be direct targets of Leu3, including most of the classically defined targets. Additional direct targets are identified by each of the methods. However, experimental and computational criteria suggest that most genes whose expression is affected by the Leu3 genotype are unlikely to be regulated by binding of the protein. Most genes that are differentially expressed by Leu3 are not direct targets despite the exceptional simplicity of the regulon, and the unusually direct nature of the perturbations investigated. These conclusions are reached through computational analyses that support and extend chromatin immunoprecipitation data on the identities of direct targets. These results have implications for the interpretation of expression experiments, especially in cases for which chromatin immunoprecipitation data are unavailable, incomplete, or ambiguous.
Carolina Digital Repository (University of North Carolina at Chapel Hill), 2008
The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recog... more The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNA-binding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes, one-third of which have not been previously reported. Several newly discovered binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially ~100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.
The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement an... more The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement and refined to an R-factor of 19.7% at a resolution of 2.1 A. This structure offers a high-resolution view of an important family of DNA-binding proteins and allows comparison to the structure of the same protein bound to DNA. The most significant difference between the current structure and that of the 2.8-A engrailed-DNA complex is the close packing of an extended strand against the rest of the protein in the unbound protein. Structural features of the protein not previously noted include a "herringbone" packing of 4 aromatic residues in the core of the protein and an extensive network of salt bridges that covers much of the helix 1-helix 2 surface. Other features that may play a role in stabilizing the native state include the interaction of buried carbonyl oxygen atoms with the edge of Phe 49 and a bias toward statistically preferred side-chain dihedral angles. There is substantial disorder at both ends of the 61 amino acid protein. A 51-amino acid variant of engrailed (residues 6-56) was synthesized and shown by CD and thermal denaturation studies to be structurally and thermodynamically similar to the full-length domain.
Carolina Digital Repository (University of North Carolina at Chapel Hill), 2006
Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in g... more Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in genomic DNA than do observed in vivo protein-DNA interactions. To determine how chromatin influences the utilization of particular DNA-binding sites, we compared the in vivo genome-wide binding location of the yeast transcription factor Leu3 to the binding location observed on the same genomic DNA in the absence of any protein cofactors. We found that the DNA-sequence motif recognized by Leu3 in vitro and in vivo was functionally indistinguishable, but Leu3 bound different genomic locations under the two conditions. Accounting for nucleosome occupancy in addition to DNA-sequence motifs significantly improved the prediction of protein-DNA interactions in vivo, but not the prediction of sites bound by purified Leu3 in vitro. Use of histone modification data does not further improve binding predictions, presumably because their effect is already manifest in the global histone distribution. Measurements of nucleosome occupancy in strains that differ in Leu3 genotype show that low nucleosome occupancy at loci bound by Leu3 is not a consequence of Leu3 binding. These results permit quantitation of the epigenetic influence that chromatin exerts on DNA binding-site selection, and provide evidence for an instructive, functionally important role for nucleosome occupancy in determining patterns of regulatory factor targeting genome-wide.
Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in g... more Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in genomic DNA than do observed in vivo protein-DNA interactions. To determine how chromatin influences the utilization of particular DNA-binding sites, we compared the in vivo genome-wide binding location of the yeast transcription factor Leu3 to the binding location observed on the same genomic DNA in the absence of any protein cofactors. We found that the DNA-sequence motif recognized by Leu3 in vitro and in vivo was functionally indistinguishable, but Leu3 bound different genomic locations under the two conditions. Accounting for nucleosome occupancy in addition to DNA-sequence motifs significantly improved the prediction of protein-DNA interactions in vivo, but not the prediction of sites bound by purified Leu3 in vitro. Use of histone modification data does not further improve binding predictions, presumably because their effect is already manifest in the global histone distribution. Measurements of nucleosome occupancy in strains that differ in Leu3 genotype show that low nucleosome occupancy at loci bound by Leu3 is not a consequence of Leu3 binding. These results permit quantitation of the epigenetic influence that chromatin exerts on DNA binding-site selection, and provide evidence for an instructive, functionally important role for nucleosome occupancy in determining patterns of regulatory factor targeting genome-wide.
Hidden Markov models representing 167 protein sequence families were used to infer the presence o... more Hidden Markov models representing 167 protein sequence families were used to infer the presence or absence of homologs within the transcriptomes of 183 algal species/strains. Statistical analyses of the distribution of HMM hits across major clades of algae, or at branch points on the phylogenetic tree of 98 chlorophytes, confirmed and extended known cases of metabolic loss and gain, most notably the loss of the mevalonate pathway for terpenoid synthesis in green algae but not, as we show here, in the streptophyte algae. Evidence for novel events was found as well, most remarkably in the recurrent and coordinated gain or loss of enzymes for the glyoxylate shunt. We find, as well, a curious pattern of retention (or re-gain) of HMG-CoA synthase in chlorophytes that have otherwise lost the mevalonate pathway, suggesting a novel, co-opted function for this enzyme in select lineages. Finally, we find striking, phylogenetically linked distributions of coding sequences for three pathways th...
More than 3 percent of the protein sequences inferred from the Caenorhabditis elegans genome cont... more More than 3 percent of the protein sequences inferred from the Caenorhabditis elegans genome contain sequence motifs characteristic of zinc-binding structural domains, and of these more than half are believed to be sequence-specific DNA-binding proteins. The distribution of these zinc-binding domains among the genomes of various organisms offers insights into the role of zinc-binding proteins in evolution. In addition, the complete genome sequence of C. elegans provides an opportunity to analyze, and perhaps predict, pathways of transcriptional regulation.
Using a physically principled method of scoring genomic sequences for the potential to be bound b... more Using a physically principled method of scoring genomic sequences for the potential to be bound by transcription factors, we have developed an algorithm for assessing the conservation of predicted binding occupancy that does not rely on sequence alignment of promoters. The method, which we call orthologweighting, assesses the degree to which the predicted binding occupancy of a transcription factor in a reference gene is also predicted in the promoters of orthologous genes. The analysis was performed separately for over 100 different transcription factors in S. cerevisiae. Statistical significance was evaluated by simulation using permuted versions of the position weight matrices. Orthologweighting produced about twice as many significantly high scoring genes as were obtained from the S. cerevisiae genome alone. Gene Ontology analysis found a similar twofold enrichment of genes. Both analyses suggest that orthologweighting improves the prediction of true regulatory targets. Interestingly, the method has only a marginal effect on the prediction of binding by chromatin immunoprecipitation (ChIP) assays. We suggest several possibilities for reconciling this result with the improved enrichment that we observe for functionally related promoters and for promoters that are under positive selection.
Journal of integrative bioinformatics, Jan 18, 2012
Advances in high throughput sequencing technology have enabled the identification of transcriptio... more Advances in high throughput sequencing technology have enabled the identification of transcription factor (TF) binding sites in genome scale. TF binding studies are important for medical applications and stem cell research. Somatic cells can be reprogrammed to a pluripotent state by the combined introduction of factors such as Oct4, Sox2, c-Myc, Klf4. These reprogrammed cells share many characteristics with embryonic stem cells (ESCs) and are known as induced pluripotent stem cells (iPSCs). The signaling requirements for maintenance of human and murine embryonic stem cells (ESCs) differ considerably. Genome wide ChIP-seq TF binding maps in mouse stem cells include Oct4, Sox2, Nanog, Tbx3, Smad2 as well as group of other factors. ChIP-seq allows study of new candidate transcription factors for reprogramming. It was shown that Nr5a2 could replace Oct4 for reprogramming. Epigenetic modifications play important role in regulation of gene expression adding additional complexity to transc...
The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recog... more The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNAbinding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes. One-third of the binding specificities have not been previously reported. Several binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially $100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.
The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement an... more The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement and refined to an R-factor of 19.7% at a resolution of 2.1 A. This structure offers a high-resolution view of an important family of DNA-binding proteins and allows comparison to the structure of the same protein bound to DNA. The most significant difference between the current structure and that of the 2.8-A engrailed-DNA complex is the close packing of an extended strand against the rest of the protein in the unbound protein. Structural features of the protein not previously noted include a "herringbone" packing of 4 aromatic residues in the core of the protein and an extensive network of salt bridges that covers much of the helix 1-helix 2 surface. Other features that may play a role in stabilizing the native state include the interaction of buried carbonyl oxygen atoms with the edge of Phe 49 and a bias toward statistically preferred side-chain dihedral angles. There is substantial disorder at both ends of the 61 amino acid protein. A 51-amino acid variant of engrailed (residues 6-56) was synthesized and shown by CD and thermal denaturation studies to be structurally and thermodynamically similar to the full-length domain.
Homeodomains are 60 amino acid DNA binding domains found in numerous eukaryotic transcription fac... more Homeodomains are 60 amino acid DNA binding domains found in numerous eukaryotic transcription factors. The homeodomain family is a useful system for studying sequence-structure relationships because several hundred sequences are known and the structures of several homeodomains have been determined. Covariation of amino acid residues in the homeodomain family has been investigated to see whether strongly covariant residue pairs can be understood in terms of the structure and function of these domains. Among 16 strongly covariant pairs examined, 2 are explained by the ability to form salt bridges, and 9 appear related to the DNA binding function of the proteins. For the remaining 5 pairs, the rationale for covariance remains unclear and the likelihood of artifactual correlations is discussed in the context of experimental and evolutionary biases in the selection of sequences. No significant correlation was found between covariance and structural proximity in the hydrophobic core.
Our group produced the best predictions overall in the DREAM3 signaling response challenge, being... more Our group produced the best predictions overall in the DREAM3 signaling response challenge, being tops by a substantial margin in the cytokine sub-challenge and nearly tied for best in the phosphoprotein sub-challenge. We achieved this success using a simple interpolation strategy. For each combination of a stimulus and inhibitor for which predictions were required, we had noted there were six other datasets using the same stimulus (but different inhibitor treatments) and six other datasets using the same inhibitor (but different stimuli). Therefore, for each treatment combination for which values were to be predicted, we calculated rank correlations for the data that were in common between the treatment combination and each of the 12 related combinations. The data from the 12 related combinations were then used to calculate missing values, weighting the contributions from each experiment based on the rank correlation coefficients. The success of this simple method suggests that the missing data were largely over-determined by similarities in the treatments. We offer some thoughts on the current state and future development of DREAM that are based on our success in this challenge, our success in the earlier DREAM2 transcription factor target challenge, and our experience as the data provider for the gene expression challenge in DREAM3.
DNA repair methyltransferases (MTases) remove methyl or other alkyl groups from the O 6 position ... more DNA repair methyltransferases (MTases) remove methyl or other alkyl groups from the O 6 position of guanine or the O 4 position of thymine by transfering the group to an active site cysteine. In order to trap an MTase-DNA complex via a disulfide bond, 2′-deoxy-6-(cystamine)-2-aminopurine (d 6 Cys 2 AP) was synthesized and incorporated into oligonucleotides. d 6 Cys 2 AP has a disulfide bond within an alkyl chain linked to the 6 position of 2,6-diaminopurine, which disulfide can be reduced to form a free thiol. Addition of human MTase to reduced oligonucleotide resulted in a protein-DNA complex that was insensitive to denaturation by SDS and high salt, but which readily dissociated in the presence of dithiothreitol. Formation of this complex was prevented by methylation of the active site cysteine. Evidence that the active site cysteine is directly involved in disulfide bond formation was obtained by N-terminal sequencing of peptides that remained associated with DNA after proteolysis of the complex.
Transcription factors that belong to the same family typically have similar, but not identical, b... more Transcription factors that belong to the same family typically have similar, but not identical, binding specificities. As such, they can be expected to compete differentially for binding to different variants of their binding sites. Pho4 is a yeast factor whose nuclear concentration is up-regulated in low phosphate, while the related factor, Cbf1, is constitutively expressed. We constructed 16 GFPreporter genes containing all palindromic variants of the motif NNCACGTGNN, and determined their activities at a range of phosphate concentrations. Pho4 affinity did not explain expression data well except under fully induced conditions. However, reporter activity was quantitatively well explained under all conditions by a model in which Cbf1 itself has modest activating activity, and Pho4 and Cbf1 compete with one another. Chromatin immunoprecipitation and computational analyses of natural Pho4 target genes, along with the activities of the reporter constructs, indicates that genes differ in their sensitivity to intermediate induction signals in part because of differences in their affinity for Cbf1. The induction sensitivity of both natural Pho4 target genes and reporter genes was well explained only by a model that assumes a role for Cbf1 in remodeling chromatin. Our analyses highlight the importance of taking into account the activities of related transcription factors in explaining system-wide gene expression data.
Several recent papers show that differences in histone modification and the use of histone varian... more Several recent papers show that differences in histone modification and the use of histone variants at the 5′ and 3′ ends of genes influence the location and kinetics of transcriptional initiation. The ultimate target of most epigenetic mechanisms may be the regulation of nucleosome occupancy, which in turn controls access to DNA at specific genomic locations.
Homeodomains are 60 amino acid DNA binding domains found in numerous eukaryotic transcription fac... more Homeodomains are 60 amino acid DNA binding domains found in numerous eukaryotic transcription factors. The homeodomain family is a useful system for studying sequence-structure relationships because several hundred sequences are known and the structures of several homeodomains have been determined. Covariation of amino acid residues in the homeodomain family has been investigated to see whether strongly covariant residue pairs can be understood in terms of the structure and function of these domains. Among 16 strongly covariant pairs examined, 2 are explained by the ability to form salt bridges, and 9 appear related to the DNA binding function of the proteins. For the remaining 5 pairs, the rationale for covariance remains unclear and the likelihood of artifactual correlations is discussed in the context of experimental and evolutionary biases in the selection of sequences. No significant correlation was found between covariance and structural proximity in the hydrophobic core.
Introduction: Opsins are a large and sequence-diverse family of light-responsive G-protein couple... more Introduction: Opsins are a large and sequence-diverse family of light-responsive G-protein coupled receptors involved in vision, circadian rhythm, and other processes. Numerous subfamilies have been defined based on sequence similarity, cell-type localization, signal transduction mechanism, or biological function, but there is no consensus classification system. Methods: We used multiple hidden Markov models (HMMs) to identify opsins in the UniProt Reference Proteomes database. Opsin-specific HMMs were also used in an annotation procedure that represents sequences as a vector of HMM scores and assess the similarity of these vectors to those of annotated sequences. UniProt Reference Proteomes are built from genome sequences, allowing us to make meaningful comparisons of the number of opsins in each of the 260 species available at the time of the survey in absolute terms and relative to a larger superfamily of which opsins are a member. Results: More than 2,000 opsins were retrieved from 262 species (all metazoans). Discussion: Merging opsin counts into higher order taxa paints a broad view of the taxonomic distribution of opsins, and of opsin subfamilies, annotated according to three different schemes.
Hidden Markov models representing 167 protein sequence families were used to infer the presence o... more Hidden Markov models representing 167 protein sequence families were used to infer the presence or absence of homologs within the transcriptomes of 183 algal species/strains. Statistical analyses of the distribution of HMM hits across major clades of algae, or at branch points on the phylogenetic tree of 98 chlorophytes, confirmed and extended known cases of metabolic loss and gain, most notably the loss of the mevalonate pathway for terpenoid synthesis in green algae but not, as we show here, in the streptophyte algae. Evidence for novel events was found as well, most remarkably in the recurrent and coordinated gain or loss of enzymes for the glyoxylate shunt. We find, as well, a curious pattern of retention (or re-gain) of HMG-CoA synthase in chlorophytes that have otherwise lost the mevalonate pathway, suggesting a novel, co-opted function for this enzyme in select lineages. Finally, we find striking, phylogenetically linked distributions of coding sequences for three pathways that synthesize the major membrane lipid phosphatidylcholine, and a complementary phylogenetic distribution pattern for the non-phospholipid DGTS (diacyl-glyceryl-trimethylhomoserine). Mass spectrometric analysis of lipids from 25 species was used to validate the inference of DGTS synthesis from sequence data. Photosynthetic eukaryotes arose from the endosymbiosis of a photosynthetic bacterium by a eukaryotic cell 1 . The descendants of this endosymbiotic event, with the exception of the land plants (Embryophyta), are all referred to as 'algae' . The algae include the Glaucophyta (a small group of freshwater species that have features thought to resemble those of the common ancestor), the Rhodophyta ('red algae' , some of which are unicellular and others, complex seaweeds), and the 'green algae' . The 'green algae' are essentially the Viridiplantae ("green plants") that are not embryophytes (land plants). They are extraordinarily diverse organisms, and by definition consist of lineages that are as old or older than the land plants. They consist of multicellular seaweeds (Ulvophyceae) as well as numerous classes of mostly unicellular organisms, Chlorophyceae and Trebouxiophyceae being the two largest groups. The green algae also include prasinophytes, a paraphyletic group that includes several classes of mostly planktonic marine species. Together, all of these groups are considered to be members of the Chlorophyta. The green algae also include the streptophyte algae, a paraphyletic group, that, together with embryophytes, form the sister group to the Chlorophyta. Most algae are unicellular, but many are not; indeed, multicellularity evolved several times in this group 2 . Given this diversity of algal morphology, and the presence of species in nearly every imaginable ecological niche, it is reasonable to suppose that there have been interesting adaptations in metabolism. The potential for metabolic novelty among the algae is of particular interest because these organisms are poised to be of increasing importance to biotechnology 3 .
Background: Effects on gene expression due to environmental or genetic changes can be easily meas... more Background: Effects on gene expression due to environmental or genetic changes can be easily measured using microarrays. However, indirect effects on expression can be substantial. The indirect effects of a perturbation need to be distinguished from the direct effects if we are to understand the structure and behavior of regulatory networks. The most direct way to perturb a transcriptional network is to alter transcription factor activity. Here, for the first time, we compare expression changes and genomic binding in a simple regulon under conditions of both low and high transcription factor activity. Specifically, we assessed the effects on expression and binding due to deletion of the yeast LEU3 transcription factor gene and effects due to elevation of Leu3 activity. Leu3 activity was elevated through overexpression and the introduction of a mutation that renders the protein constitutively active. Genes that are bound and/or regulated by Leu3 under one or both conditions were characterized in terms of their functional annotations and their predicted potential to be bound by Leu3. We also assessed the evolutionary conservation of the predicted binding potential using a novel alignment-independent method. Both perturbations yield genes that are likely to be direct targets of Leu3, including most of the classically defined targets. Additional direct targets are identified by each of the methods. However, experimental and computational criteria suggest that most genes whose expression is affected by the Leu3 genotype are unlikely to be regulated by binding of the protein. Most genes that are differentially expressed by Leu3 are not direct targets despite the exceptional simplicity of the regulon, and the unusually direct nature of the perturbations investigated. These conclusions are reached through computational analyses that support and extend chromatin immunoprecipitation data on the identities of direct targets. These results have implications for the interpretation of expression experiments, especially in cases for which chromatin immunoprecipitation data are unavailable, incomplete, or ambiguous.
Carolina Digital Repository (University of North Carolina at Chapel Hill), 2008
The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recog... more The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNA-binding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes, one-third of which have not been previously reported. Several newly discovered binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially ~100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.
The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement an... more The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement and refined to an R-factor of 19.7% at a resolution of 2.1 A. This structure offers a high-resolution view of an important family of DNA-binding proteins and allows comparison to the structure of the same protein bound to DNA. The most significant difference between the current structure and that of the 2.8-A engrailed-DNA complex is the close packing of an extended strand against the rest of the protein in the unbound protein. Structural features of the protein not previously noted include a "herringbone" packing of 4 aromatic residues in the core of the protein and an extensive network of salt bridges that covers much of the helix 1-helix 2 surface. Other features that may play a role in stabilizing the native state include the interaction of buried carbonyl oxygen atoms with the edge of Phe 49 and a bias toward statistically preferred side-chain dihedral angles. There is substantial disorder at both ends of the 61 amino acid protein. A 51-amino acid variant of engrailed (residues 6-56) was synthesized and shown by CD and thermal denaturation studies to be structurally and thermodynamically similar to the full-length domain.
Carolina Digital Repository (University of North Carolina at Chapel Hill), 2006
Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in g... more Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in genomic DNA than do observed in vivo protein-DNA interactions. To determine how chromatin influences the utilization of particular DNA-binding sites, we compared the in vivo genome-wide binding location of the yeast transcription factor Leu3 to the binding location observed on the same genomic DNA in the absence of any protein cofactors. We found that the DNA-sequence motif recognized by Leu3 in vitro and in vivo was functionally indistinguishable, but Leu3 bound different genomic locations under the two conditions. Accounting for nucleosome occupancy in addition to DNA-sequence motifs significantly improved the prediction of protein-DNA interactions in vivo, but not the prediction of sites bound by purified Leu3 in vitro. Use of histone modification data does not further improve binding predictions, presumably because their effect is already manifest in the global histone distribution. Measurements of nucleosome occupancy in strains that differ in Leu3 genotype show that low nucleosome occupancy at loci bound by Leu3 is not a consequence of Leu3 binding. These results permit quantitation of the epigenetic influence that chromatin exerts on DNA binding-site selection, and provide evidence for an instructive, functionally important role for nucleosome occupancy in determining patterns of regulatory factor targeting genome-wide.
Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in g... more Sequence motifs that are potentially recognized by DNA-binding proteins occur far more often in genomic DNA than do observed in vivo protein-DNA interactions. To determine how chromatin influences the utilization of particular DNA-binding sites, we compared the in vivo genome-wide binding location of the yeast transcription factor Leu3 to the binding location observed on the same genomic DNA in the absence of any protein cofactors. We found that the DNA-sequence motif recognized by Leu3 in vitro and in vivo was functionally indistinguishable, but Leu3 bound different genomic locations under the two conditions. Accounting for nucleosome occupancy in addition to DNA-sequence motifs significantly improved the prediction of protein-DNA interactions in vivo, but not the prediction of sites bound by purified Leu3 in vitro. Use of histone modification data does not further improve binding predictions, presumably because their effect is already manifest in the global histone distribution. Measurements of nucleosome occupancy in strains that differ in Leu3 genotype show that low nucleosome occupancy at loci bound by Leu3 is not a consequence of Leu3 binding. These results permit quantitation of the epigenetic influence that chromatin exerts on DNA binding-site selection, and provide evidence for an instructive, functionally important role for nucleosome occupancy in determining patterns of regulatory factor targeting genome-wide.
Hidden Markov models representing 167 protein sequence families were used to infer the presence o... more Hidden Markov models representing 167 protein sequence families were used to infer the presence or absence of homologs within the transcriptomes of 183 algal species/strains. Statistical analyses of the distribution of HMM hits across major clades of algae, or at branch points on the phylogenetic tree of 98 chlorophytes, confirmed and extended known cases of metabolic loss and gain, most notably the loss of the mevalonate pathway for terpenoid synthesis in green algae but not, as we show here, in the streptophyte algae. Evidence for novel events was found as well, most remarkably in the recurrent and coordinated gain or loss of enzymes for the glyoxylate shunt. We find, as well, a curious pattern of retention (or re-gain) of HMG-CoA synthase in chlorophytes that have otherwise lost the mevalonate pathway, suggesting a novel, co-opted function for this enzyme in select lineages. Finally, we find striking, phylogenetically linked distributions of coding sequences for three pathways th...
More than 3 percent of the protein sequences inferred from the Caenorhabditis elegans genome cont... more More than 3 percent of the protein sequences inferred from the Caenorhabditis elegans genome contain sequence motifs characteristic of zinc-binding structural domains, and of these more than half are believed to be sequence-specific DNA-binding proteins. The distribution of these zinc-binding domains among the genomes of various organisms offers insights into the role of zinc-binding proteins in evolution. In addition, the complete genome sequence of C. elegans provides an opportunity to analyze, and perhaps predict, pathways of transcriptional regulation.
Using a physically principled method of scoring genomic sequences for the potential to be bound b... more Using a physically principled method of scoring genomic sequences for the potential to be bound by transcription factors, we have developed an algorithm for assessing the conservation of predicted binding occupancy that does not rely on sequence alignment of promoters. The method, which we call orthologweighting, assesses the degree to which the predicted binding occupancy of a transcription factor in a reference gene is also predicted in the promoters of orthologous genes. The analysis was performed separately for over 100 different transcription factors in S. cerevisiae. Statistical significance was evaluated by simulation using permuted versions of the position weight matrices. Orthologweighting produced about twice as many significantly high scoring genes as were obtained from the S. cerevisiae genome alone. Gene Ontology analysis found a similar twofold enrichment of genes. Both analyses suggest that orthologweighting improves the prediction of true regulatory targets. Interestingly, the method has only a marginal effect on the prediction of binding by chromatin immunoprecipitation (ChIP) assays. We suggest several possibilities for reconciling this result with the improved enrichment that we observe for functionally related promoters and for promoters that are under positive selection.
Journal of integrative bioinformatics, Jan 18, 2012
Advances in high throughput sequencing technology have enabled the identification of transcriptio... more Advances in high throughput sequencing technology have enabled the identification of transcription factor (TF) binding sites in genome scale. TF binding studies are important for medical applications and stem cell research. Somatic cells can be reprogrammed to a pluripotent state by the combined introduction of factors such as Oct4, Sox2, c-Myc, Klf4. These reprogrammed cells share many characteristics with embryonic stem cells (ESCs) and are known as induced pluripotent stem cells (iPSCs). The signaling requirements for maintenance of human and murine embryonic stem cells (ESCs) differ considerably. Genome wide ChIP-seq TF binding maps in mouse stem cells include Oct4, Sox2, Nanog, Tbx3, Smad2 as well as group of other factors. ChIP-seq allows study of new candidate transcription factors for reprogramming. It was shown that Nr5a2 could replace Oct4 for reprogramming. Epigenetic modifications play important role in regulation of gene expression adding additional complexity to transc...
The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recog... more The sequence specificity of DNA-binding proteins is the primary mechanism by which the cell recognizes genomic features. Here, we describe systematic determination of yeast transcription factor DNAbinding specificities. We obtained binding specificities for 112 DNA-binding proteins representing 19 distinct structural classes. One-third of the binding specificities have not been previously reported. Several binding sequences have striking genomic distributions relative to transcription start sites, supporting their biological relevance and suggesting a role in promoter architecture. Among these are Rsc3 binding sequences, containing the core CGCG, which are found preferentially $100 bp upstream of transcription start sites. Mutation of RSC3 results in a dramatic increase in nucleosome occupancy in hundreds of proximal promoters containing a Rsc3 binding element, but has little impact on promoters lacking Rsc3 binding sequences, indicating that Rsc3 plays a broad role in targeting nucleosome exclusion at yeast promoters.
The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement an... more The structure of the Drosophila engrailed homeodomain has been solved by molecular replacement and refined to an R-factor of 19.7% at a resolution of 2.1 A. This structure offers a high-resolution view of an important family of DNA-binding proteins and allows comparison to the structure of the same protein bound to DNA. The most significant difference between the current structure and that of the 2.8-A engrailed-DNA complex is the close packing of an extended strand against the rest of the protein in the unbound protein. Structural features of the protein not previously noted include a "herringbone" packing of 4 aromatic residues in the core of the protein and an extensive network of salt bridges that covers much of the helix 1-helix 2 surface. Other features that may play a role in stabilizing the native state include the interaction of buried carbonyl oxygen atoms with the edge of Phe 49 and a bias toward statistically preferred side-chain dihedral angles. There is substantial disorder at both ends of the 61 amino acid protein. A 51-amino acid variant of engrailed (residues 6-56) was synthesized and shown by CD and thermal denaturation studies to be structurally and thermodynamically similar to the full-length domain.
Homeodomains are 60 amino acid DNA binding domains found in numerous eukaryotic transcription fac... more Homeodomains are 60 amino acid DNA binding domains found in numerous eukaryotic transcription factors. The homeodomain family is a useful system for studying sequence-structure relationships because several hundred sequences are known and the structures of several homeodomains have been determined. Covariation of amino acid residues in the homeodomain family has been investigated to see whether strongly covariant residue pairs can be understood in terms of the structure and function of these domains. Among 16 strongly covariant pairs examined, 2 are explained by the ability to form salt bridges, and 9 appear related to the DNA binding function of the proteins. For the remaining 5 pairs, the rationale for covariance remains unclear and the likelihood of artifactual correlations is discussed in the context of experimental and evolutionary biases in the selection of sequences. No significant correlation was found between covariance and structural proximity in the hydrophobic core.
Uploads
Papers by Neil Clarke