Browsing by Subject "Saccharomyces cerevisiae Proteins"
Now showing 1 - 20 of 22
Results Per Page
Sort Options
Item Open Access A nucleosome-guided map of transcription factor binding sites in yeast.(PLoS Comput Biol, 2007-11) Narlikar, Leelavati; Gordân, Raluca; Hartemink, Alexander JFinding functional DNA binding sites of transcription factors (TFs) throughout the genome is a crucial step in understanding transcriptional regulation. Unfortunately, these binding sites are typically short and degenerate, posing a significant statistical challenge: many more matches to known TF motifs occur in the genome than are actually functional. However, information about chromatin structure may help to identify the functional sites. In particular, it has been shown that active regulatory regions are usually depleted of nucleosomes, thereby enabling TFs to bind DNA in those regions. Here, we describe a novel motif discovery algorithm that employs an informative prior over DNA sequence positions based on a discriminative view of nucleosome occupancy. When a Gibbs sampling algorithm is applied to yeast sequence-sets identified by ChIP-chip, the correct motif is found in 52% more cases with our informative prior than with the commonly used uniform prior. This is the first demonstration that nucleosome occupancy information can be used to improve motif discovery. The improvement is dramatic, even though we are using only a statistical model to predict nucleosome occupancy; we expect our results to improve further as high-resolution genome-wide experimental nucleosome occupancy data becomes increasingly available.Item Open Access Amino acid permeases require COPII components and the ER resident membrane protein Shr3p for packaging into transport vesicles in vitro.(J Cell Biol, 1996-11) Kuehn, MJ; Schekman, R; Ljungdahl, POIn S. cerevisiae lacking SHR3, amino acid permeases specifically accumulate in membranes of the endoplasmic reticulum (ER) and fail to be transported to the plasma membrane. We examined the requirements of transport of the permeases from the ER to the Golgi in vitro. Addition of soluble COPII components (Sec23/24p, Sec13/31p, and Sar1p) to yeast membrane preparations generated vesicles containing the general amino acid permease. Gap1p, and the histidine permease, Hip1p. Shr3p was required for the packaging of Gap1p and Hip1p but was not itself incorporated into transport vesicles. In contrast, the packaging of the plasma membrane ATPase, Pma1p, and the soluble yeast pheromone precursor, glycosylated pro alpha factor, was independent of Shr3p. In addition, we show that integral membrane and soluble cargo colocalize in transport vesicles, indicating that different types of cargo are not segregated at an early step in secretion. Our data suggest that specific ancillary proteins in the ER membrane recruit subsets of integral membrane protein cargo into COPII transport vesicles.Item Open Access Calmodulin dissociation regulates Myo5 recruitment and function at endocytic sites.(EMBO J, 2010-09-01) Grötsch, Helga; Giblin, Jonathan P; Idrissi, Fatima-Zahra; Fernández-Golbano, Isabel-María; Collette, John R; Newpher, Thomas M; Robles, Virginia; Lemmon, Sandra K; Geli, María-IsabelMyosins-I are conserved proteins that bear an N-terminal motor head followed by a Tail Homology 1 (TH1) lipid-binding domain. Some myosins-I have an additional C-terminal extension (C(ext)) that promotes Arp2/3 complex-dependent actin polymerization. The head and the tail are separated by a neck that binds calmodulin or calmodulin-related light chains. Myosins-I are known to participate in actin-dependent membrane remodelling. However, the molecular mechanisms controlling their recruitment and their biochemical activities in vivo are far from being understood. In this study, we provided evidence suggesting the existence of an inhibitory interaction between the TH1 domain of the yeast myosin-I Myo5 and its C(ext). The TH1 domain prevented binding of the Myo5 C(ext) to the yeast WIP homologue Vrp1, Myo5 C(ext)-induced actin polymerization and recruitment of the Myo5 C(ext) to endocytic sites. Our data also indicated that calmodulin dissociation from Myo5 weakened the interaction between the neck and TH1 domains and the C(ext). Concomitantly, calmodulin dissociation triggered Myo5 binding to Vrp1, extended the myosin-I lifespan at endocytic sites and activated Myo5-induced actin polymerization.Item Open Access Cell-cycle control of cell polarity in yeast.(The Journal of cell biology, 2019-01) Moran, Kyle D; Kang, Hui; Araujo, Ana V; Zyla, Trevin R; Saito, Koji; Tsygankov, Denis; Lew, Daniel JIn many cells, morphogenetic events are coordinated with the cell cycle by cyclin-dependent kinases (CDKs). For example, many mammalian cells display extended morphologies during interphase but round up into more spherical shapes during mitosis (high CDK activity) and constrict a furrow during cytokinesis (low CDK activity). In the budding yeast Saccharomyces cerevisiae, bud formation reproducibly initiates near the G1/S transition and requires activation of CDKs at a point called "start" in G1. Previous work suggested that CDKs acted by controlling the ability of cells to polarize Cdc42, a conserved Rho-family GTPase that regulates cell polarity and the actin cytoskeleton in many systems. However, we report that yeast daughter cells can polarize Cdc42 before CDK activation at start. This polarization operates via a positive feedback loop mediated by the Cdc42 effector Ste20. We further identify a major and novel locus of CDK action downstream of Cdc42 polarization, affecting the ability of several other Cdc42 effectors to localize to the polarity site.Item Open Access Chemotactic movement of a polarity site enables yeast cells to find their mates.(Proceedings of the National Academy of Sciences of the United States of America, 2021-06) Ghose, Debraj; Jacobs, Katherine; Ramirez, Samuel; Elston, Timothy; Lew, DanielHow small eukaryotic cells can interpret dynamic, noisy, and spatially complex chemical gradients to orient growth or movement is poorly understood. We address this question using Saccharomyces cerevisiae, where cells orient polarity up pheromone gradients during mating. Initial orientation is often incorrect, but polarity sites then move around the cortex in a search for partners. We find that this movement is biased by local pheromone gradients across the polarity site: that is, movement of the polarity site is chemotactic. A bottom-up computational model recapitulates this biased movement. The model reveals how even though pheromone-bound receptors do not mimic the shape of external pheromone gradients, nonlinear and stochastic effects combine to generate effective gradient tracking. This mechanism for gradient tracking may be applicable to any cell that searches for a target in a complex chemical landscape.Item Open Access Chromatin: bind at your own RSC.(Curr Biol, 2011-03-22) Buchler, Nicolas E; Bai, LuRecent work has identified a novel RSC-nucleosome complex that both strongly phases flanking nucleosomes and presents regulatory sites for ready access. These results challenge several widely held views.Item Open Access Clathrin is important for normal actin dynamics and progression of Sla2p-containing patches during endocytosis in yeast.(Traffic (Copenhagen, Denmark), 2006-05) Newpher, Thomas M; Lemmon, Sandra KClathrin is a major vesicle coat protein involved in receptor-mediated endocytosis. In yeast and higher eukaryotes, clathrin is recruited to the plasma membrane during the early stage of endocytosis along with clathrin-associated adaptors. As coated pits undergo maturation, a burst of actin polymerization accompanies and helps drive vesicle internalization. Here, we investigate the dynamics of clathrin relative to the early endocytic patch protein Sla2p. We find that clathrin is recruited to the cortex prior to Sla2p. In the absence of clathrin, normal numbers of Sla2p patches form, but many do not internalize or are dramatically delayed in completion of endocytosis. Patches that do internalize receive Sla1p late, which is followed by Abp1, which appears near the end of Sla2p lifetime. In addition, clathrin mutants develop actin comet tails, suggesting an important function in actin patch organization/dynamics. Similar to its mammalian counterparts, the light chain (LC) subunit of yeast clathrin interacts directly with the coiled-coil domain of Sla2p. A mutant of Sla2p that no longer interacts with LC (sla2Delta376-573) results in delayed progression of endocytic patches and aberrant actin dynamics. These data demonstrate an important role for clathrin in organization and progression of early endocytic patches to the late stages of endocytosis.Item Open Access Different Mechanisms Confer Gradual Control and Memory at Nutrient- and Stress-Regulated Genes in Yeast.(Mol Cell Biol, 2015-11) Rienzo, Alessandro; Poveda-Huertes, Daniel; Aydin, Selcan; Buchler, Nicolas E; Pascual-Ahuir, Amparo; Proft, MarkusCells respond to environmental stimuli by fine-tuned regulation of gene expression. Here we investigated the dose-dependent modulation of gene expression at high temporal resolution in response to nutrient and stress signals in yeast. The GAL1 activity in cell populations is modulated in a well-defined range of galactose concentrations, correlating with a dynamic change of histone remodeling and RNA polymerase II (RNAPII) association. This behavior is the result of a heterogeneous induction delay caused by decreasing inducer concentrations across the population. Chromatin remodeling appears to be the basis for the dynamic GAL1 expression, because mutants with impaired histone dynamics show severely truncated dose-response profiles. In contrast, the GRE2 promoter operates like a rapid off/on switch in response to increasing osmotic stress, with almost constant expression rates and exclusively temporal regulation of histone remodeling and RNAPII occupancy. The Gal3 inducer and the Hog1 mitogen-activated protein (MAP) kinase seem to determine the different dose-response strategies at the two promoters. Accordingly, GAL1 becomes highly sensitive and dose independent if previously stimulated because of residual Gal3 levels, whereas GRE2 expression diminishes upon repeated stimulation due to acquired stress resistance. Our analysis reveals important differences in the way dynamic signals create dose-sensitive gene expression outputs.Item Unknown DNA mismatches reveal conformational penalties in protein-DNA recognition.(Nature, 2020-11) Afek, Ariel; Shi, Honglue; Rangadurai, Atul; Sahay, Harshit; Senitzki, Alon; Xhani, Suela; Fang, Mimi; Salinas, Raul; Mielko, Zachery; Pufall, Miles A; Poon, Gregory MK; Haran, Tali E; Schumacher, Maria A; Al-Hashimi, Hashim M; Gordân, RalucaTranscription factors recognize specific genomic sequences to regulate complex gene-expression programs. Although it is well-established that transcription factors bind to specific DNA sequences using a combination of base readout and shape recognition, some fundamental aspects of protein-DNA binding remain poorly understood1,2. Many DNA-binding proteins induce changes in the structure of the DNA outside the intrinsic B-DNA envelope. However, how the energetic cost that is associated with distorting the DNA contributes to recognition has proven difficult to study, because the distorted DNA exists in low abundance in the unbound ensemble3-9. Here we use a high-throughput assay that we term SaMBA (saturation mismatch-binding assay) to investigate the role of DNA conformational penalties in transcription factor-DNA recognition. In SaMBA, mismatched base pairs are introduced to pre-induce structural distortions in the DNA that are much larger than those induced by changes in the Watson-Crick sequence. Notably, approximately 10% of mismatches increased transcription factor binding, and for each of the 22 transcription factors that were examined, at least one mismatch was found that increased the binding affinity. Mismatches also converted non-specific sites into high-affinity sites, and high-affinity sites into 'super sites' that exhibit stronger affinity than any known canonical binding site. Determination of high-resolution X-ray structures, combined with nuclear magnetic resonance measurements and structural analyses, showed that many of the DNA mismatches that increase binding induce distortions that are similar to those induced by protein binding-thus prepaying some of the energetic cost incurred from deforming the DNA. Our work indicates that conformational penalties are a major determinant of protein-DNA recognition, and reveals mechanisms by which mismatches can recruit transcription factors and thus modulate replication and repair activities in the cell10,11.Item Unknown Environmental and genetic determinants of colony morphology in yeast.(PLoS Genet, 2010-01-22) Granek, JA; Magwene, PMNutrient stresses trigger a variety of developmental switches in the budding yeast Saccharomyces cerevisiae. One of the least understood of such responses is the development of complex colony morphology, characterized by intricate, organized, and strain-specific patterns of colony growth and architecture. The genetic bases of this phenotype and the key environmental signals involved in its induction have heretofore remained poorly understood. By surveying multiple strain backgrounds and a large number of growth conditions, we show that limitation for fermentable carbon sources coupled with a rich nitrogen source is the primary trigger for the colony morphology response in budding yeast. Using knockout mutants and transposon-mediated mutagenesis, we demonstrate that two key signaling networks regulating this response are the filamentous growth MAP kinase cascade and the Ras-cAMP-PKA pathway. We further show synergistic epistasis between Rim15, a kinase involved in integration of nutrient signals, and other genes in these pathways. Ploidy, mating-type, and genotype-by-environment interactions also appear to play a role in the controlling colony morphology. Our study highlights the high degree of network reuse in this model eukaryote; yeast use the same core signaling pathways in multiple contexts to integrate information about environmental and physiological states and generate diverse developmental outputs.Item Unknown How cells determine the number of polarity sites.(eLife, 2021-04-26) Chiou, Jian-Geng; Moran, Kyle D; Lew, Daniel JThe diversity of cell morphologies arises, in part, through regulation of cell polarity by Rho-family GTPases. A poorly understood but fundamental question concerns the regulatory mechanisms by which different cells generate different numbers of polarity sites. Mass-conserved activator-substrate (MCAS) models that describe polarity circuits develop multiple initial polarity sites, but then those sites engage in competition, leaving a single winner. Theoretical analyses predicted that competition would slow dramatically as GTPase concentrations at different polarity sites increase toward a 'saturation point', allowing polarity sites to coexist. Here, we test this prediction using budding yeast cells, and confirm that increasing the amount of key polarity proteins results in multiple polarity sites and simultaneous budding. Further, we elucidate a novel design principle whereby cells can switch from competition to equalization among polarity sites. These findings provide insight into how cells with diverse morphologies may determine the number of polarity sites.Item Open Access How the kinetochore couples microtubule force and centromere stretch to move chromosomes.(Nature cell biology, 2016-04) Suzuki, Aussie; Badger, Benjamin L; Haase, Julian; Ohashi, Tomoo; Erickson, Harold P; Salmon, Edward D; Bloom, KerryThe Ndc80 complex (Ndc80, Nuf2, Spc24 and Spc25) is a highly conserved kinetochore protein essential for end-on anchorage to spindle microtubule plus ends and for force generation coupled to plus-end polymerization and depolymerization. Spc24/Spc25 at one end of the Ndc80 complex binds the kinetochore. The N-terminal tail and CH domains of Ndc80 bind microtubules, and an internal domain binds microtubule-associated proteins (MAPs) such as the Dam1 complex. To determine how the microtubule- and MAP-binding domains of Ndc80 contribute to force production at the kinetochore in budding yeast, we have inserted a FRET tension sensor into the Ndc80 protein about halfway between its microtubule-binding and internal loop domains. The data support a mechanical model of force generation at metaphase where the position of the kinetochore relative to the microtubule plus end reflects the relative strengths of microtubule depolymerization, centromere stretch and microtubule-binding interactions with the Ndc80 and Dam1 complexes.Item Open Access In vivo dynamics of clathrin and its adaptor-dependent recruitment to the actin-based endocytic machinery in yeast.(Developmental cell, 2005-07) Newpher, Thomas M; Smith, Robin P; Lemmon, Vance; Lemmon, Sandra KClathrin-mediated transport is a major pathway for endocytosis. However, in yeast, where cortical actin patches are essential for endocytosis, plasma membrane-associated clathrin has never been observed. Using live cell imaging, we demonstrate cortical clathrin in association with the actin-based endocytic machinery in yeast. Fluorescently tagged clathrin is found in highly mobile internal trans-Golgi/endosomal structures and in smaller cortical patches. Total internal reflection fluorescence microscopy showed that cortical patches are likely endocytic sites, as clathrin is recruited prior to a burst of intensity of the actin patch/endocytic marker, Abp1. Clathrin also accumulates at the cortex with internalizing alpha factor receptor, Ste2p. Cortical clathrin localizes with epsins Ent1/2p and AP180s, and its recruitment to the surface is dependent upon these adaptors. In contrast, Sla2p, End3p, Pan1p, and a dynamic actin cytoskeleton are not required for clathrin assembly or exchange but are required for the mobility, maturation, and/or turnover of clathrin-containing endocytic structures.Item Open Access Measuring fast gene dynamics in single cells with time-lapse luminescence microscopy.(Mol Biol Cell, 2014-11-05) Mazo-Vargas, Anyimilehidi; Park, Heungwon; Aydin, Mert; Buchler, Nicolas ETime-lapse fluorescence microscopy is an important tool for measuring in vivo gene dynamics in single cells. However, fluorescent proteins are limited by slow chromophore maturation times and the cellular autofluorescence or phototoxicity that arises from light excitation. An alternative is luciferase, an enzyme that emits photons and is active upon folding. The photon flux per luciferase is significantly lower than that for fluorescent proteins. Thus time-lapse luminescence microscopy has been successfully used to track gene dynamics only in larger organisms and for slower processes, for which more total photons can be collected in one exposure. Here we tested green, yellow, and red beetle luciferases and optimized substrate conditions for in vivo luminescence. By combining time-lapse luminescence microscopy with a microfluidic device, we tracked the dynamics of cell cycle genes in single yeast with subminute exposure times over many generations. Our method was faster and in cells with much smaller volumes than previous work. Fluorescence of an optimized reporter (Venus) lagged luminescence by 15-20 min, which is consistent with its known rate of chromophore maturation in yeast. Our work demonstrates that luciferases are better than fluorescent proteins at faithfully tracking the underlying gene expression.Item Open Access Mechanisms that ensure monogamous mating in Saccharomyces cerevisiae.(Molecular biology of the cell, 2021-04) Robertson, Corrina G; Clark-Cotton, Manuella R; Lew, Daniel JHaploid cells of the budding yeast Saccharomyces cerevisiae communicate using secreted pheromones and mate to form diploid zygotes. Mating is monogamous, resulting in the fusion of precisely one cell of each mating type. Monogamous mating in crowded conditions, where cells have access to more than one potential partner, raises the question of how multiple-mating outcomes are prevented. Here we identify mutants capable of mating with multiple partners, revealing the mechanisms that ensure monogamous mating. Before fusion, cells develop polarity foci oriented toward potential partners. Competition between these polarity foci within each cell leads to disassembly of all but one focus, thus favoring a single fusion event. Fusion promotes the formation of heterodimeric complexes between subunits that are uniquely expressed in each mating type. One complex shuts off haploid-specific gene expression, and the other shuts off the ability to respond to pheromone. Zygotes able to form either complex remain monogamous, but zygotes lacking both can re-mate.Item Open Access Novel function of clathrin light chain in promoting endocytic vesicle formation.(Molecular biology of the cell, 2006-10) Newpher, Thomas M; Idrissi, Fatima-Zahra; Geli, Maria Isabel; Lemmon, Sandra KClathrin-mediated endocytosis is a major pathway for uptake of lipid and protein cargo at the plasma membrane. The lattices of clathrin-coated pits and vesicles are comprised of triskelions, each consisting of three oligomerized heavy chains (HC) bound by a light chain (LC). In addition to binding HC, LC interacts with members of the Hip1/R family of endocytic proteins, including the budding yeast homologue, Sla2p. Here, using in vivo analysis in yeast, we provide novel insight into the role of this interaction. We find that overexpression of LC partially restores endocytosis to cells lacking clathrin HC. This suppression is dependent on the Sla2p binding region of LC. Using live cell imaging techniques to visualize endocytic vesicle formation, we find that the N-terminal Sla2p binding region of LC promotes the progression of arrested Sla2p patches that form in the absence of HC. We propose that LC binding to Sla2p positively regulates Sla2p for efficient endocytic vesicle formation.Item Open Access Ratiometric GPCR signaling enables directional sensing in yeast.(PLoS biology, 2019-10-17) Henderson, Nicholas T; Pablo, Michael; Ghose, Debraj; Clark-Cotton, Manuella R; Zyla, Trevin R; Nolen, James; Elston, Timothy C; Lew, Daniel JAccurate detection of extracellular chemical gradients is essential for many cellular behaviors. Gradient sensing is challenging for small cells, which can experience little difference in ligand concentrations on the up-gradient and down-gradient sides of the cell. Nevertheless, the tiny cells of the yeast Saccharomyces cerevisiae reliably decode gradients of extracellular pheromones to find their mates. By imaging the behavior of polarity factors and pheromone receptors, we quantified the accuracy of initial polarization during mating encounters. We found that cells bias the orientation of initial polarity up-gradient, even though they have unevenly distributed receptors. Uneven receptor density means that the gradient of ligand-bound receptors does not accurately reflect the external pheromone gradient. Nevertheless, yeast cells appear to avoid being misled by responding to the fraction of occupied receptors rather than simply the concentration of ligand-bound receptors. Such ratiometric sensing also serves to amplify the gradient of active G protein. However, this process is quite error-prone, and initial errors are corrected during a subsequent indecisive phase in which polarity clusters exhibit erratic mobile behavior.Item Open Access Recombinational Repair of Nuclease-Generated Mitotic Double-Strand Breaks with Different End Structures in Yeast.(G3 (Bethesda, Md.), 2020-10) Gamble, Dionna; Shaltz, Samantha; Jinks-Robertson, SueMitotic recombination is the predominant mechanism for repairing double-strand breaks in Saccharomyces cerevisiae Current recombination models are largely based on studies utilizing the enzyme I-SceI or HO to create a site-specific break, each of which generates broken ends with 3' overhangs. In this study sequence-diverged ectopic substrates were used to assess whether the frequent Pol δ-mediated removal of a mismatch 8 nucleotides from a 3' end affects recombination outcomes and whether the presence of a 3' vs. 5' overhang at the break site alters outcomes. Recombination outcomes monitored were the distributions of recombination products into crossovers vs. noncrossovers, and the position/length of transferred sequence (heteroduplex DNA) in noncrossover products. A terminal mismatch that was 22 nucleotides from the 3' end was rarely removed and the greater distance from the end did not affect recombination outcomes. To determine whether the recombinational repair of breaks with 3' vs. 5' overhangs differs, we compared the well-studied 3' overhang created by I-SceI to a 5' overhang created by a ZFN (Zinc Finger Nuclease). Initiation with the ZFN yielded more recombinants, consistent with more efficient cleavage and potentially faster repair rate relative to I-SceI. While there were proportionally more COs among ZFN- than I-SceI-initiated events, NCOs in the two systems were indistinguishable in terms of the extent of strand transfer. These data demonstrate that the method of DSB induction and the resulting differences in end polarity have little effect on mitotic recombination outcomes despite potential differences in repair rate.Item Open Access Role of Scd5, a protein phosphatase-1 targeting protein, in phosphoregulation of Sla1 during endocytosis.(J Cell Sci, 2012-10-15) Chi, Richard J; Torres, Onaidy T; Segarra, Verónica A; Lansley, Tanya; Chang, Ji Suk; Newpher, Thomas M; Lemmon, Sandra KPhosphorylation regulates assembly and disassembly of proteins during endocytosis. In yeast, Prk1 and Ark1 phosphorylate factors after vesicle internalization leading to coat disassembly. Scd5, a protein phosphatase-1 (PP1)-targeting subunit, is proposed to regulate dephosphorylation of Prk1/Ark1 substrates to promote new rounds of endocytosis. In this study we analyzed scd5-PP1Δ2, a mutation causing impaired PP1 binding. scd5-PP1Δ2 caused hyperphosphorylation of several Prk1 endocytic targets. Live-cell imaging of 15 endocytic components in scd5-PP1Δ2 revealed that most factors arriving before the invagination/actin phase of endocytosis had delayed lifetimes. Severely affected were early factors and Sla2 (Hip1R homolog), whose lifetime was extended nearly fourfold. In contrast, the lifetime of Sla1, a Prk1 target, was extended less than twofold, but its cortical recruitment was significantly reduced. Delayed Sla2 dynamics caused by scd5-PP1Δ2 were suppressed by SLA1 overexpression. This was dependent on the LxxQxTG repeats (SR) of Sla1, which are phosphorylated by Prk1 and bind Pan1, another Prk1 target, in the dephosphorylated state. Without the SR, Sla1ΔSR was still recruited to the cell surface, but was less concentrated in cortical patches than Pan1. sla1ΔSR severely impaired endocytic progression, but this was partially suppressed by overexpression of LAS17, suggesting that without the SR region the SH3 region of Sla1 causes constitutive negative regulation of Las17 (WASp). These results demonstrate that Scd5/PP1 is important for recycling Prk1 targets to initiate new rounds of endocytosis and provide new mechanistic information on the role of the Sla1 SR domain in regulating progression to the invagination/actin phase of endocytosis.Item Open Access Stability selection for regression-based models of transcription factor-DNA binding specificity.(Bioinformatics, 2013-07-01) Mordelet, Fantine; Horton, John; Hartemink, Alexander J; Engelhardt, Barbara E; Gordân, RalucaMOTIVATION: The DNA binding specificity of a transcription factor (TF) is typically represented using a position weight matrix model, which implicitly assumes that individual bases in a TF binding site contribute independently to the binding affinity, an assumption that does not always hold. For this reason, more complex models of binding specificity have been developed. However, these models have their own caveats: they typically have a large number of parameters, which makes them hard to learn and interpret. RESULTS: We propose novel regression-based models of TF-DNA binding specificity, trained using high resolution in vitro data from custom protein-binding microarray (PBM) experiments. Our PBMs are specifically designed to cover a large number of putative DNA binding sites for the TFs of interest (yeast TFs Cbf1 and Tye7, and human TFs c-Myc, Max and Mad2) in their native genomic context. These high-throughput quantitative data are well suited for training complex models that take into account not only independent contributions from individual bases, but also contributions from di- and trinucleotides at various positions within or near the binding sites. To ensure that our models remain interpretable, we use feature selection to identify a small number of sequence features that accurately predict TF-DNA binding specificity. To further illustrate the accuracy of our regression models, we show that even in the case of paralogous TF with highly similar position weight matrices, our new models can distinguish the specificities of individual factors. Thus, our work represents an important step toward better sequence-based models of individual TF-DNA binding specificity. AVAILABILITY: Our code is available at http://genome.duke.edu/labs/gordan/ISMB2013. The PBM data used in this article are available in the Gene Expression Omnibus under accession number GSE47026.