A single amino acid can determine the DNA binding specificity of homeodomain proteins

Missense mutations in CRX homeodomain cause dominant retinopathies through two distinct mechanisms

Homeodomain transcription factors (HD TFs) are instrumental to vertebrate development. Mutations in HD TFs have been linked to human diseases, but their pathogenic mechanisms remain elusive. Here, we use Cone-Rod Homeobox (CRX) as a model to decipher the disease-causing mechanisms of two HD mutations, p.E80A and p.K88N, that produce severe dominant retinopathies. Through integrated analysis of molecular and functional evidence in vitro and in knock-in mouse models, we uncover two novel gain-of-function mechanisms: p.E80A increases CRX-mediated transactivation of canonical CRX target genes in developing photoreceptors; p.K88N alters CRX DNA-binding specificity resulting in binding at ectopic sites and severe perturbation of CRX target gene expression. Both mechanisms produce novel retinal morphological defects and hinder photoreceptor maturation distinct from loss-of-function models. This study reveals the distinct roles of E80 and K88 residues in CRX HD regulatory functions and emphasizes the importance of transcriptional precision in normal development.

Missense mutations in CRX homeodomain cause dominant retinopathies through two distinct mechanisms

Homeodomain transcription factors (HD TFs) are instrumental to vertebrate development. Mutations in HD TFs have been linked to human diseases, but their pathogenic mechanisms remain elusive. Here we use Cone-Rod Homeobox (CRX) as a model to decipher the disease-causing mechanisms of two HD mutations, p.E80A and p.K88N, that produce severe dominant retinopathies. Through integrated analysis of molecular and functional evidence in vitro and in knock-in mouse models, we uncover two novel gain-of-function mechanisms: p.E80A increases CRX-mediated transactivation of canonical CRX target genes in developing photoreceptors; p.K88N alters CRX DNA-binding specificity resulting in binding at ectopic sites and severe perturbation of CRX target gene expression. Both mechanisms produce novel retinal morphological defects and hinder photoreceptor maturation distinct from loss-of-function models. This study reveals the distinct roles of E80 and K88 residues in CRX HD regulatory functions and emphasizes the importance of transcriptional precision in normal development.

Intra-Domain Residue Coevolution in Transcription Factors Contributes to DNA Binding Specificity

Understanding the basis of the DNA-binding specificity of transcription factors (TFs) has been of long-standing interest. Despite extensive efforts to map millions of putative TF binding sequences, identifying the critical determinants for DNA binding specificity remains a major challenge. The coevolution of residues in proteins occurs due to a shared evolutionary history. However, it is unclear how coevolving residues in TFs contribute to DNA binding specificity. Here, we systematically collected publicly available data sets from multiple large-scale high-throughput TF-DNA interaction screening experiments for the major TF families with large numbers of TF members. These families included the Homeobox, HLH, bZIP_1, Ets, HMG_box, ZF-C4, and Zn_clus TFs. We detected TF subclass-determining sites (TSDSs) and showed that the TSDSs were more likely to coevolve with other TSDSs than with non-TSDSs, particularly for the Homeobox, HLH, Ets, bZIP_1, and HMG_box TF families. By in silico modeling, we showed that mutation of the highly coevolving residues could significantly reduce the stability of the TF-DNA complex. The distant residues from the DNA interface also contributed to TF-DNA binding activity. Overall, our study gave evidence that coevolved residues relate to transcriptional regulation and provided insights into the potential application of engineered DNA-binding domains and proteins. IMPORTANCE While unraveling DNA-binding specificity of TFs is the key to understanding the basis and molecular mechanism of gene expression regulation, identifying the critical determinants that contribute to DNA binding specificity remains a major challenge. In this study, we provided evidence showing that coevolving residues in TF domains contributed to DNA binding specificity. We demonstrated that the TSDSs were more likely to coevolve with other TSDSs than with non-TSDSs. Mutation of the coevolving residue pairs (CRPs) could significantly reduce the stability of THE TF-DNA complex, and even the distant residues from the DNA interface contribute to TF-DNA binding activity. Collectively, our study expands our knowledge of the interactions among coevolved residues in TFs, tertiary contacting, and functional importance in refined transcriptional regulation. Understanding the impact of coevolving residues in TFs will help understand the details of transcription of gene regulation and advance the application of engineered DNA-binding domains and protein.

Hox genes collaborate with helix-loop-helix factor Grainyhead to promote neuroblast apoptosis along the anterior-posterior axis of the Drosophila larval central nervous system

Hox genes code for a family of a homeodomain (HD) containing transcription factors that use TALE-HD containing factors Pbx/Exd and Meis/Hth to specify the development of the anterior-posterior (AP) axis of an organism. However, the absence of TALE-HD containing factors from specific tissues emphasizes the need to identify and validate new Hox cofactors. In Drosophila central nervous system (CNS), Hox execute segment-specific apoptosis of neural stem cells (neuroblasts-NBs) and neurons. In abdominal segments of larval CNS, Hox gene Abdominal-A (AbdA) mediates NB apoptosis with the help of Exd and bHLH factor Grainyhead (Grh) using a 717 bp apoptotic enhancer. In this study, we show that this enhancer is critical for abdominal NB apoptosis and relies on two separable set of DNA binding motifs responsible for its initiation and maintenance. Our results also show that AbdA and Grh interact through their highly conserved DNA binding domains, and the DNA binding specificity of AbdA-HD is important for it to interact with Grh and essential for it to execute NB apoptosis in CNS. We also establish that Grh is required for Hox-dependent NB apoptosis in Labial and Sex Combs Reduced (Scr) expressing regions of the CNS, and it can physically interact with all the Hox proteins in vitro. Our biochemical and functional data collectively support the idea that Grh can function as a Hox cofactor and help them carry out their in vivo roles during development.

Synthetic reconstruction of the hunchback promoter specifies the role of Bicoid, Zelda and Hunchback in the dynamics of its transcription

For over 40 years, the Bicoid-hunchback (Bcd-hb) system in the fruit fly embryo has been used as a model to study how positional information in morphogen concentration gradients is robustly translated into step-like responses. A body of quantitative comparisons between theory and experiment have since questioned the initial paradigm that the sharp hb transcription pattern emerges solely from diffusive biochemical interactions between the Bicoid transcription factor and the gene promoter region. Several alternative mechanisms have been proposed, such as additional sources of positional information, positive feedback from Hb proteins or out-of-equilibrium transcription activation. By using the MS2-MCP RNA-tagging system and analysing in real time, the transcription dynamics of synthetic reporters for Bicoid and/or its two partners Zelda and Hunchback, we show that all the early hb expression pattern features and temporal dynamics are compatible with an equilibrium model with a short decay length Bicoid activity gradient as a sole source of positional information. Meanwhile, Bicoid’s partners speed-up the process by different means: Zelda lowers the Bicoid concentration threshold required for transcriptional activation while Hunchback reduces burstiness and increases the polymerase firing rate.

Deciphering the Retinal Epigenome during Development, Disease and Reprogramming: Advancements, Challenges and Perspectives

Retinal neurogenesis is driven by concerted actions of transcription factors, some of which are expressed in a continuum and across several cell subtypes throughout development. While seemingly redundant, many factors diversify their regulatory outcome on gene expression, by coordinating variations in chromatin landscapes to drive divergent retinal specification programs. Recent studies have furthered the understanding of the epigenetic contribution to the progression of age-related macular degeneration, a leading cause of blindness in the elderly. The knowledge of the epigenomic mechanisms that control the acquisition and stabilization of retinal cell fates and are evoked upon damage, holds the potential for the treatment of retinal degeneration. Herein, this review presents the state-of-the-art approaches to investigate the retinal epigenome during development, disease, and reprogramming. A pipeline is then reviewed to functionally interrogate the epigenetic and transcriptional networks underlying cell fate specification, relying on a truly unbiased screening of open chromatin states. The related work proposes an inferential model to identify gene regulatory networks, features the first footprinting analysis and the first tentative, systematic query of candidate pioneer factors in the retina ever conducted in any model organism, leading to the identification of previously uncharacterized master regulators of retinal cell identity, such as the nuclear factor I, NFI. This pipeline is virtually applicable to the study of genetic programs and candidate pioneer factors in any developmental context. Finally, challenges and limitations intrinsic to the current next-generation sequencing techniques are discussed, as well as recent advances in super-resolution imaging, enabling spatio-temporal resolution of the genome.

The power of the (imperfect) palindrome: Sequence-specific roles of palindromic motifs in gene regulation

In human languages, a palindrome reads the same forward as backward (e.g., 'madam'). In regulatory DNA, a palindrome is an inverted sequence repeat that allows a transcription factor to bind as a homodimer or as a heterodimer with another type of transcription factor. Regulatory palindromes are typically imperfect, that is, the repeated sequences differ in at least one base pair, but the functional significance of this asymmetry remains poorly understood. Here, we review the use of imperfect palindromes in Drosophila photoreceptor differentiation and mammalian steroid receptor signaling. Moreover, we discuss mechanistic explanations for the predominance of imperfect palindromes over perfect palindromes in these two gene regulatory contexts. Lastly, we propose to elucidate whether specific imperfectly palindromic variants have specific regulatory functions in steroid receptor signaling and whether such variants can help predict transcriptional outcomes as well as the response of individual patients to drug treatments.

Suboptimal Intermediates Underlie Evolution of the Bicoid Homeodomain

Changes in regulatory networks generate materials for evolution to create phenotypic diversity. For transcription networks, multiple studies have shown that alterations in binding sites of cis-regulatory elements correlate well with the gain or loss of specific features of the body plan. Less is known about alterations in the amino acid sequences of the transcription factors (TFs) that bind these elements. Here we study the evolution of Bicoid (Bcd), a homeodomain (HD) protein that is critical for anterior embryo patterning in Drosophila. The ancestor of Bcd (AncBcd) emerged after a duplication of a Zerknullt (Zen)-like ancestral protein (AncZB) in a suborder of flies. AncBcd diverged from AncZB, gaining novel transcriptional and translational activities. We focus on the evolution of the HD of AncBcd, which binds to DNA and RNA, and is comprised of four subdomains: an N-terminal arm (NT) and three helices; H1, H2, and Recognition Helix (RH). Using chimeras of subdomains and gene rescue assays in Drosophila, we show that robust patterning activity of the Bcd HD (high frequency rescue to adulthood) is achieved only when amino acid substitutions in three separate subdomains (NT, H1, and RH) are combined. Other combinations of subdomains also yield full rescue, but with lower penetrance, suggesting alternative suboptimal activities. Our results suggest a multistep pathway for the evolution of the Bcd HD that involved intermediate HD sequences with suboptimal activities, which constrained and enabled further evolutionary changes. They also demonstrate critical epistatic forces that contribute to the robust function of a DNA-binding domain.

Molecular topography of an entire nervous system

We have produced gene expression profiles of all 302 neurons of the C. elegans nervous system that match the single-cell resolution of its anatomy and wiring diagram. Our results suggest that individual neuron classes can be solely identified by combinatorial expression of specific gene families. For example, each neuron class expresses distinct codes of ∼23 neuropeptide genes and ∼36 neuropeptide receptors, delineating a complex and expansive "wireless" signaling network. To demonstrate the utility of this comprehensive gene expression catalog, we used computational approaches to (1) identify cis-regulatory elements for neuron-specific gene expression and (2) reveal adhesion proteins with potential roles in process placement and synaptic specificity. Our expression data are available at https://cengen.org and can be interrogated at the web application CengenApp. We expect that this neuron-specific directory of gene expression will spur investigations of underlying mechanisms that define anatomy, connectivity, and function throughout the C. elegans nervous system.

A combinatorial cis-regulatory logic restricts color-sensing Rhodopsins to specific photoreceptor subsets in Drosophila

Color vision in Drosophila melanogaster is based on the expression of five different color-sensing Rhodopsin proteins in distinct subtypes of photoreceptor neurons. Promoter regions of less than 300 base pairs are sufficient to reproduce the unique, photoreceptor subtype-specific rhodopsin expression patterns. The underlying cis-regulatory logic remains poorly understood, but it has been proposed that the rhodopsin promoters have a bipartite structure: the distal promoter region directs the highly restricted expression in a specific photoreceptor subtype, while the proximal core promoter region provides general activation in all photoreceptors. Here, we investigate whether the rhodopsin promoters exhibit a strict specialization of their distal (subtype specificity) and proximal (general activation) promoter regions, or if both promoter regions contribute to generating the photoreceptor subtype-specific expression pattern. To distinguish between these two models, we analyze the expression patterns of a set of hybrid promoters that combine the distal promoter region of one rhodopsin with the proximal core promoter region of another rhodopsin. We find that the function of the proximal core promoter regions extends beyond providing general activation: these regions play a previously underappreciated role in generating the non-overlapping expression patterns of the different rhodopsins. Therefore, cis-regulatory motifs in both the distal and the proximal core promoter regions recruit transcription factors that generate the unique rhodopsin patterns in a combinatorial manner. We compare this combinatorial regulatory logic to the regulatory logic of olfactory receptor genes and discuss potential implications for the evolution of rhodopsins.

Each type of sensory receptor neuron in our body expresses a specific sensory receptor protein, which allows us to detect and discriminate a variety of environmental stimuli. The regulatory logic that controls this spatially precise and highly restricted expression of sensory receptor proteins remains poorly understood. As a model system, we study the mechanisms that control the expression of different color-sensing Rhodopsin proteins in distinct subtypes of Drosophila photoreceptors, which is the basis for color vision. Compact promoter regions of less than 300 base pairs are sufficient to reproduce the non-overlapping rhodopsin patterns. However, the regulatory logic that underlies the combination (sometimes called ‘grammar’) of the cis-regulatory motifs (sometimes called ‘vocabulary’) within the rhodopsin promoters remains poorly understood. Here, we find that specific combinations of cis-regulatory motifs in the distal and the proximal core promoter regions of each rhodopsin direct its unique expression pattern.

Systematic comparison and prediction of the effects of missense mutations on protein-DNA and protein-RNA interactions

The binding affinities of protein-nucleic acid interactions could be altered due to missense mutations occurring in DNA- or RNA-binding proteins, therefore resulting in various diseases. Unfortunately, a systematic comparison and prediction of the effects of mutations on protein-DNA and protein-RNA interactions (these two mutation classes are termed MPDs and MPRs, respectively) is still lacking. Here, we demonstrated that these two classes of mutations could generate similar or different tendencies for binding free energy changes in terms of the properties of mutated residues. We then developed regression algorithms separately for MPDs and MPRs by introducing novel geometric partition-based energy features and interface-based structural features. Through feature selection and ensemble learning, similar computational frameworks that integrated energy- and nonenergy-based models were established to estimate the binding affinity changes resulting from MPDs and MPRs, but the selected features for the final models were different and therefore reflected the specificity of these two mutation classes. Furthermore, the proposed methodology was extended to the identification of mutations that significantly decreased the binding affinities. Extensive validations indicated that our algorithm generally performed better than the state-of-the-art methods on both the regression and classification tasks. The webserver and software are freely available at http://liulab.hzau.edu.cn/PEMPNI and https://github.com/hzau-liulab/PEMPNI.

Protein-nucleic acid interactions play important roles in various cellular processes. Missense mutations occurring in DNA- or RNA-binding proteins (termed MPDs and MPRs, respectively) could change the binding affinities of these interactions. Previous studies have compared protein-DNA and protein-RNA interactions from multifaceted viewpoints, but less attention has been given to the similarities and specific differences between the effects of MPDs and MPRs and between the methodologies for predicting the affinity changes induced by the two mutation classes. Therefore, we systematically compared their impacts and demonstrated that MPDs and MPRs could have specific preferences for binding affinity changes. These observations motivated us to construct regression models separately for MPDs and MPRs by introducing novel energy and nonenergy descriptors. Although similar frameworks were developed to estimate these two categories of mutation effects, different descriptors were selected in the regression models and further revealed the specificity of mutation classes. The interplay between the energy and nonenergy modules effectively improved prediction performance. Our algorithm can also be adopted to disentangle mutations significantly decreasing binding affinities from other mutations.

A six-amino-acid motif is a major determinant in functional evolution of HOX1 proteins

Gene duplication and divergence is a major driver in the emergence of evolutionary novelties. How variations in amino acid sequences lead to loss of ancestral activity and functional diversification of proteins is poorly understood. We used cross-species functional analysis of Drosophila Labial and its mouse HOX1 orthologs (HOXA1, HOXB1, and HOXD1) as a paradigm to address this issue. Mouse HOX1 proteins display low (30%) sequence similarity with Drosophila Labial. However, substituting endogenous Labial with the mouse proteins revealed that HOXA1 has retained essential ancestral functions of Labial, while HOXB1 and HOXD1 have diverged. Genome-wide analysis demonstrated similar DNA-binding patterns of HOXA1 and Labial in mouse cells, while HOXB1 binds to distinct targets. Compared with HOXB1, HOXA1 shows an enrichment in co-occupancy with PBX proteins on target sites and exists in the same complex with PBX on chromatin. Functional analysis of HOXA1-HOXB1 chimeric proteins uncovered a novel six-amino-acid C-terminal motif (CTM) flanking the homeodomain that serves as a major determinant of ancestral activity. In vitro DNA-binding experiments and structural prediction show that CTM provides an important domain for interaction of HOXA1 proteins with PBX. Our findings show that small changes outside of highly conserved DNA-binding regions can lead to profound changes in protein function.

The transcription factor PITX1 drives astrocyte differentiation by regulating the SOX9 gene

Astrocytes perform multiple essential functions in the developing and mature brain, including regulation of synapse formation, control of neurotransmitter release and uptake, and maintenance of extracellular ion balance. As a result, astrocytes have been implicated in the progression of neurodegenerative disorders such as Alzheimer's disease, Huntington's disease, and Parkinson's disease. Despite these critical functions, the study of human astrocytes can be difficult because standard differentiation protocols are time-consuming and technically challenging, but a differentiation protocol recently developed in our laboratory enables the efficient derivation of astrocytes from human embryonic stem cells. We used this protocol along with microarrays, luciferase assays, electrophoretic mobility shift assays, and ChIP assays to explore the genes involved in astrocyte differentiation. We demonstrate that paired-like homeodomain transcription factor 1 (PITX1) is critical for astrocyte differentiation. PITX1 overexpression induced early differentiation of astrocytes, and its knockdown blocked astrocyte differentiation. PITX1 overexpression also increased and PITX1 knockdown decreased expression of sex-determining region Y box 9 (SOX9), known initiator of gliogenesis, during early astrocyte differentiation. Moreover, we determined that PITX1 activates the SOX9 promoter through a unique binding motif. Taken together, these findings indicate that PITX1 drives astrocyte differentiation by sustaining activation of the SOX9 promoter.

Structural basis for the complex DNA binding behavior of the plant stem cell regulator WUSCHEL

Stem cells are one of the foundational evolutionary novelties that allowed the independent emergence of multicellularity in the plant and animal lineages. In plants, the homeodomain (HD) transcription factor WUSCHEL (WUS) is essential for the maintenance of stem cells in the shoot apical meristem. WUS has been reported to bind to diverse DNA motifs and to act as transcriptional activator and repressor. However, the mechanisms underlying this remarkable behavior have remained unclear. Here, we quantitatively delineate WUS binding to three divergent DNA motifs and resolve the relevant structural underpinnings. We show that WUS exhibits a strong binding preference for TGAA repeat sequences, while retaining the ability to weakly bind to TAAT elements. This behavior is attributable to the formation of dimers through interactions of specific residues in the HD that stabilize WUS DNA interaction. Our results provide a mechanistic basis for dissecting WUS dependent regulatory networks in plant stem cell control.

Constraints and limitations on the transcriptional response downstream of the Bicoid morphogen gradient

The regulation of the hunchback promoter expression by the maternal Bicoid gradient has been studied as a model system in development for many years. Yet, at the level of quantitative agreement between data and theoretical models, even the first step of this regulation, transcription, continues to be challenging. This situation is slowly progressing, thanks to quantitative live-imaging techniques coupled to advanced statistical data analysis and modeling. Here, we outline the current state of our knowledge of this apparently "simple" step, highlighting the newly appreciated role of bursty transcription dynamics and its regulation.

One evolutionarily selected amino acid variation is sufficient to provide functional specificity in the cold shock protein paralogs of Staphylococcus aureus

Bacterial genomes encode several families of protein paralogs. Discrimination between functional divergence and redundancy among paralogs is challenging due to their sequence conservation. Here, we investigated whether the amino acid differences present in the cold shock protein (CSP) paralogs of Staphylococcus aureus were responsible for functional specificity. Since deletion of cspA reduces the synthesis of staphyloxanthin (STX), we used it as an in vivo reporter of CSP functionality. Complementation of a ΔcspA strain with the different S. aureus CSP variants showed that only CspA could specifically restore STX production by controlling the activity of the stress‐associated sigma B factor (σ^B). To determine the amino acid residues responsible for CspA specificity, we created several chimeric CSPs that interchanged the amino acid differences between CspA and CspC, which shared the highest identity. We demonstrated that CspA Pro58 was responsible for the specific control of σ^B activity and its associated phenotypes. Interestingly, CspC gained the biological function of CspA when the E58P substitution was introduced. This study highlights how just one evolutionarily selected amino acid change may be sufficient to modify the specific functionality of CSP paralogs.

Identification of genes involved in glaucoma pathogenesis using combined network analysis and empirical studies

Glaucoma is a leading cause of blindness. We aimed in this study to identify genes that may make subtle and cumulative contributions to glaucoma pathogenesis. To this end, we identified molecular interactions and pathways that include transcription factors (TFs) FOXC1, PITX2, PAX6 and NFKB1 and various microRNAs including miR-204 known to have relevance to trabecular meshwork (TM) functions and/or glaucoma. TM tissue is involved in glaucoma pathogenesis. In-house microarray transcriptome results and data sources were used to identify target genes of the regulatory molecules. Bioinformatics analyses were done to filter TM and glaucoma relevant genes. These were submitted to network-creating softwares to define interactions, pathways and a network that would include the genes. The network was stringently scrutinized and minimized, then expanded by addition of microarray data and data on TF and microRNA-binding sites. Selected features of the network were confirmed by empirical studies such as dual luciferase assays, real-time PCR and western blot experiments and apoptosis assays. MYOC, WDR36, LTPBP2, RHOA, CYP1B1, OPA1, SPARC, MEIS2, PLEKHG5, RGS5, BBS5, ALDH1A1, NOMO2, CXCL6, FMNL2, ADAMTS5, CLOCK and DKK1 were among the genes included in the final network. Pathways identified included those that affect ECM properties, IOP, ciliary body functions, retinal ganglion cell viability, apoptosis, focal adhesion and oxidative stress response. The identification of many genes potentially involved in glaucoma pathology is consistent with its being a complex disease. The inclusion of several known glaucoma-related genes validates the approach used.

Cis-regulatory basis of sister cell type divergence in the vertebrate retina

Multicellular organisms evolved via repeated functional divergence of transcriptionally related sister cell types, but the mechanisms underlying sister cell type divergence are not well understood. Here, we study a canonical pair of sister cell types, retinal photoreceptors and bipolar cells, to identify the key cis-regulatory features that distinguish them. By comparing open chromatin maps and transcriptomic profiles, we found that while photoreceptor and bipolar cells have divergent transcriptomes, they share remarkably similar cis-regulatory grammars, marked by enrichment of K50 homeodomain binding sites. However, cell class-specific enhancers are distinguished by enrichment of E-box motifs in bipolar cells, and Q50 homeodomain motifs in photoreceptors. We show that converting K50 motifs to Q50 motifs represses reporter expression in bipolar cells, while photoreceptor expression is maintained. These findings suggest that partitioning of Q50 motifs within cell type-specific cis-regulatory elements was a critical step in the evolutionary divergence of the bipolar transcriptome from that of photoreceptors.

Brain-specific homeobox Bsx specifies identity of pineal gland between serially homologous photoreceptive organs in zebrafish

The pineal gland functioning as a photoreceptive organ in non-mammalian species is a serial homolog of the retina. Here we found that Brain-specific homeobox (Bsx) is a key regulator conferring individuality on the pineal gland between the two serially homologous photoreceptive organs in zebrafish. Bsx knock-down impaired the pineal development with reduced expression of exorh, the pineal-specific gene responsible for the photoreception, whereas it induced ectopic expression of rho, a retina-specific gene, in the pineal gland. Bsx remarkably transactivated the exorh promoter in combination with Otx5, but not with Crx, through its binding to distinct subtypes of PIRE, a DNA cis-element driving Crx/Otx-dependent pineal-specific gene expression. These results demonstrate that the identity of pineal photoreceptive neurons is determined by the combinatorial code of Bsx and Otx5, the former confers the pineal specificity at the tissue level and the latter determines the photoreceptor specificity at the cellular level.

Predicting protein-DNA binding free energy change upon missense mutations using modified MM/PBSA approach: SAMPDI webserver

Motivation

Protein-DNA interactions are essential for regulating many cellular processes, such as transcription, replication, recombination and translation. Amino acid mutations occurring in DNA-binding proteins have profound effects on protein-DNA binding and are linked with many diseases. Hence, accurate and fast predictions of the effects of mutations on protein-DNA binding affinity are essential for understanding disease-causing mechanisms and guiding plausible treatments.

Results

Here we report a new method Single Amino acid Mutation binding free energy change of Protein-DNA Interaction (SAMPDI). The method utilizes modified Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) approach along with an additional set of knowledge-based terms delivered from investigations of the physicochemical properties of protein-DNA complexes. The method is benchmarked against experimentally determined binding free energy changes caused by 105 mutations in 13 proteins (compiled ProNIT database and data from recent references), and results in correlation coefficient of 0.72.

Availability and implementation

http://compbio.clemson.edu/SAMPDI.

Contact

ealexov@clemson.edu.

Supplementary information

Supplementary data are available at Bioinformatics online.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

Diversification of transcription factor-DNA interactions and the evolution of gene regulatory networks

Sequence-specific transcription factors (TFs) bind short DNA sequences in the genome to regulate the expression of target genes. In the last decade, numerous technical advances have enabled the determination of the DNA-binding specificities of many of these factors. Large-scale screens of many TFs enabled the creation of databases of TF DNA-binding specificities, typically represented as position weight matrices (PWMs). Although great progress has been made in determining and predicting binding specificities systematically, there are still many surprises to be found when studying a particular TF's interactions with DNA in detail. Paralogous TFs' binding specificities can differ in subtle ways, in a manner that is not immediately apparent from looking at their PWMs. These differences affect gene regulatory outputs and enable TFs to rewire transcriptional networks over evolutionary time. This review discusses recent observations made in the study of TF-DNA interactions that highlight the importance of continued in-depth analysis of TF-DNA interactions and their inherent complexity. This article is categorized under: Biological Mechanisms > Regulatory Biology.

A feed-forward relay integrates the regulatory activities of Bicoid and Orthodenticle via sequential binding to suboptimal sites

Datta et al. define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bicoid (Bcd) and Orthodenticle (Otd). The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins.

The K50 (lysine at amino acid position 50) homeodomain (HD) protein Orthodenticle (Otd) is critical for anterior patterning and brain and eye development in most metazoans. In Drosophila melanogaster, another K50HD protein, Bicoid (Bcd), has evolved to replace Otd's ancestral function in embryo patterning. Bcd is distributed as a long-range maternal gradient and activates transcription of a large number of target genes, including otd. Otd and Bcd bind similar DNA sequences in vitro, but how their transcriptional activities are integrated to pattern anterior regions of the embryo is unknown. Here we define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bcd and Otd. Class 1 enhancers are initially activated by Bcd, and activation is transferred to Otd via a feed-forward relay (FFR) that involves sequential binding of the two proteins to the same DNA motif. Class 2 enhancers are activated by Bcd and maintained by an Otd-independent mechanism. Class 3 enhancers are never bound by Bcd, but Otd binds and activates them in a second wave of zygotic transcription. The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins. Our results define specific patterning roles for Bcd and Otd and provide mechanisms for coordinating the precise timing of gene expression patterns during embryonic development.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

Combinatorial regulation of a Blimp1 (Prdm1) enhancer in the mouse retina

The mouse retina comprises seven major cell types that exist in differing proportions. They are generated from multipotent progenitors in a stochastic manner, such that the relative frequency of any given type generated changes over time. The mechanisms determining the proportions of each cell type are only partially understood. Photoreceptors and bipolar interneurons are derived from cells that express Otx2. Within this population, Blimp1 (Prdm1) helps set the balance between photoreceptors and bipolar cells by suppressing bipolar identity in most of the cells. How only a subset of these Otx2+ cells decides to upregulate Blimp1 and adopt photoreceptor fate is unknown. To understand this, we investigated how Blimp1 transcription is regulated. We identified several potential Blimp1 retinal enhancer elements using DNase hypersensitivity sequencing. Only one of the elements recapitulated Blimp1 spatial and temporal expression in cultured explant assays and within the retinas of transgenic mice. Mutagenesis of this retinal Blimp1 enhancer element revealed four discrete sequences that were each required for its activity. These included highly conserved Otx2 and ROR (retinoic acid receptor related orphan receptor) binding sites. The other required sequences do not appear to be controlled by Otx2 or ROR factors, increasing the complexity of the Blimp1 gene regulatory network. Our results show that the intersection of three or more transcription factors is required to correctly regulate the spatial and temporal features of Blimp1 enhancer expression. This explains how Blimp1 expression can diverge from Otx2 and set the balance between photoreceptor and bipolar fates.

Hierarchical mutational events compensate for glutamate auxotrophy of a Bacillus subtilis gltC mutant

Glutamate is the major donor of nitrogen for anabolic reactions. The Gram-positive soil bacterium Bacillus subtilis either utilizes exogenously provided glutamate or synthesizes it using the gltAB-encoded glutamate synthase (GOGAT). In the absence of glutamate, the transcription factor GltC activates expression of the GOGAT genes for glutamate production. Consequently, a gltC mutant strain is auxotrophic for glutamate. Using a genetic selection and screening system, we could isolate and differentiate between gltC suppressor mutants in one step. All mutants had acquired the ability to synthesize glutamate, independent of GltC. We identified (i) gain-of-function mutations in the gltR gene, encoding the transcription factor GltR, (ii) mutations in the promoter of the gltAB operon and (iii) massive amplification of the genomic locus containing the gltAB operon. The mutants belonging to the first two classes constitutively expressed the gltAB genes and produced sufficient glutamate for growth. By contrast, mutants that belong to the third class appeared most frequently and solved glutamate limitation by increasing the copy number of the poorly expressed gltAB genes. Thus, glutamate auxotrophy of a B. subtilis gltC mutant can be relieved in multiple ways. Moreover, recombination-dependent amplification of the gltAB genes is the predominant mutational event indicating a hierarchy of mutations.

Expression of CXCL6 and BBS5 that may be glaucoma relevant genes is regulated by PITX2

The transcription factor PITX2 is implicated in glaucoma pathology. In an earlier study we had used microarray analysis to identify genes in the trabecular meshwork (TM) that are affected by knock down of PITX2. Here, those studies were pursued to identify genes that are direct targets of PITX2 and that may be relevant to glaucoma. Initially, bioinformatics tools were used to select among the genes that had been affected by PITX2 knock down those that have PITX2 binding sites and that may be involved in glaucoma related functions. Subsequently, the effect of PITX2 was tested using the dual luciferase assay in four cell cultures including two primary TM cultures co-transfected with vectors containing promoter fragments of six candidate genes upstream of a luciferase gene and a vector that expressed PITX2. Finally, the effect of PITX2 on endogenous expression of two genes was assessed by over expression and knock down of PITX2 in TM cells. Thirty four genes were found to contain PITX2 binding sites in their putative promoter regions, and 16 were found to be associated with TM-specific and/or glaucoma associated functions. Results of dual luciferase assays confirmed that two of six genes tested were directly targeted by PITX2. The two genes were CXCL6 (chemokine (C-X-C motif) ligand 6) and BBS5 (Bardet-Biedl syndrome 5). Over expression and knock down of PITX2 showed that this transcription factor affects endogenous expression of these two genes in TM cells. CXCL6 encodes a pro-inflammatory cytokine, and many studies have suggested that cytokines and other immune system functions are involved in glaucoma pathogenesis. BBS5 is a member of the BBS family of genes that affect ciliary functions, and ciliary bodies in the anterior chamber of the eye produce the aqueous fluid that affects intraocular pressure. Immune related functions and intraocular pressure are both important components of glaucoma pathology. The role of PITX2 in glaucoma may be mediated partly by regulating the expression of CXCL6 and BBS5 and thus affecting immune functions and intraocular pressure.

60 YEARS OF POMC: Transcriptional and epigenetic regulation of POMC gene expression

Expression of the pro-opiomelanocortin (POMC) gene integrates numerous inputs that reflect the developmental history of POMC-expressing cells of the pituitary and hypothalamus, as well as their critical role in the endocrine system. These inputs are integrated at specific regulatory sequences within the promoter and pituitary or hypothalamic enhancers of the POMC locus. Investigations of developmental mechanisms and transcription factors (TFs) responsible for pituitary activation of POMC transcription led to the discovery of the Pitx factors that have critical roles in pituitary development and striking patterning functions in embryonic development. Terminal differentiation of the two pituitary POMC lineages, the corticotrophs and melanotrophs, is controlled by Tpit; mutations of the human TPIT gene cause isolated adrenocorticotrophic hormone deficiency. Intermediate lobe and melanotroph identity is provided by the pioneer TF Pax7 that remodels chromatin to reveal a new repertoire of enhancers for Tpit action. Many signaling pathways regulate POMC transcription including activation by hypothalamic corticotrophin-releasing hormone acting through the orphan nuclear receptors of the Nur family and feedback repression by glucocorticoids and their glucocorticoid receptor. TFs of the basic helix-loop-helix, Smad, Stat, Etv, and nuclear factor-B families also mediate signals for control of POMC transcription. Whereas most of these regulatory processes are conserved in different species, there are also notable differences between specific targets for regulation of the human compared with mouse POMC genes.

HOPX: The Unusual Homeodomain-Containing Protein

The homeodomain-only protein homeobox (HOPX) is the smallest known member of the homeodomain-containing protein family, atypically unable to bind DNA. HOPX is widely expressed in diverse tissues, where it is critically involved in the regulation of proliferation and differentiation. In human skin, HOPX controls epidermal formation through the regulation of late differentiation markers, and HOPX expression correlates with the level of differentiation in cutaneous pathologies. In mouse skin, Hopx was additionally identified as a lineage tracing marker of quiescent hair follicle stem cells. This review discusses current knowledge of HOPX structure and function in normal and pathological conditions.

Structure of human heat-shock transcription factor 1 in complex with DNA

Heat-shock transcription factor 1 (HSF1) has a central role in mediating the protective response to protein conformational stresses in eukaryotes. HSF1 consists of an N-terminal DNA-binding domain (DBD), a coiled-coil oligomerization domain, a regulatory domain and a transactivation domain. Upon stress, HSF1 trimerizes via its coiled-coil domain and binds to the promoters of heat shock protein-encoding genes. Here, we present cocrystal structures of the human HSF1 DBD in complex with cognate DNA. A comparative analysis of the HSF1 paralog Skn7 from Chaetomium thermophilum showed that single amino acid changes in the DBD can switch DNA binding specificity, thus revealing the structural basis for the interaction of HSF1 with cognate DNA. We used a crystal structure of the coiled-coil domain of C. thermophilum Skn7 to develop a model of the active human HSF1 trimer in which HSF1 embraces the heat-shock-element DNA.

Dermestes maculatus: an intermediate-germ beetle model system for evo-devo

BACKGROUND

Understanding how genes change during evolution to direct the development of diverse body plans is a major goal of the evo-devo field. Achieving this will require the establishment of new model systems that represent key points in phylogeny. These new model systems must be amenable to laboratory culture, and molecular and functional approaches should be feasible. To date, studies of insects have been best represented by the model system Drosophila melanogaster. Given the enormous diversity represented by insect taxa, comparative studies within this clade will provide a wealth of information about the evolutionary potential and trajectories of alternative developmental strategies.

RESULTS

Here we established the beetle Dermestes maculatus, a member of the speciose clade Coleoptera, as a new insect model system. We have maintained a continuously breeding culture in the lab and documented Dermestes maculatus embryogenesis using nuclear and phalloidin staining. Anterior segments are specified during the blastoderm stage before gastrulation, and posterior segments are added sequentially during germ band elongation. We isolated and studied the expression and function of the pair-rule segmentation gene paired in Dermestes maculatus. In this species, paired is expressed in stripes during both blastoderm and germ band stages: four primary stripes arise prior to gastrulation, confirming an intermediate-germ mode of development for this species. As in other insects, these primary stripes then split into secondary stripes. To study gene function, we established both embryonic and parental RNAi. Knockdown of Dmac-paired with either method resulted in pair-rule-like segmentation defects, including loss of Engrailed expression in alternate stripes.

CONCLUSIONS

These studies establish basic approaches necessary to use Dermestes maculatus as a model system. Methods are now available for use of this intermediate-germ insect for future studies of the evolution of regulatory networks controlling insect segmentation, as well as of other processes in development and homeostasis. Consistent with the role of paired in long-germ Drosophila and shorter-germ Tribolium, paired functions as a pair-rule segmentation gene in Dermestes maculatus. Thus, paired retains pair-rule function in insects with different modes of segment addition.

Role of DNA conformation & energetic insights in Msx-1-DNA recognition as revealed by molecular dynamics studies on specific and nonspecific complexes

In most of homeodomain-DNA complexes, glutamine or lysine is present at 50th position and interacts with 5th and 6th nucleotide of core recognition region. Molecular dynamics simulations of Msx-1-DNA complex (Q50-TG) and its variant complexes, that is specific (Q50K-CC), nonspecific (Q50-CC) having mutation in DNA and (Q50K-TG) in protein, have been carried out. Analysis of protein-DNA interactions and structure of DNA in specific and nonspecific complexes show that amino acid residues use sequence-dependent shape of DNA to interact. The binding free energies of all four complexes were analysed to define role of amino acid residue at 50th position in terms of binding strength considering the variation in DNA on stability of protein-DNA complexes. The order of stability of protein-DNA complexes shows that specific complexes are more stable than nonspecific ones. Decomposition analysis shows that N-terminal amino acid residues have been found to contribute maximally in binding free energy of protein-DNA complexes. Among specific protein-DNA complexes, K50 contributes more as compared to Q50 towards binding free energy in respective complexes. The sequence dependence of local conformation of DNA enables Q50/Q50K to make hydrogen bond with nucleotide(s) of DNA. The changes in amino acid sequence of protein are accommodated and stabilized around TAAT core region of DNA having variation in nucleotides.

Impacts of the ubiquitous factor Zelda on Bicoid-dependent DNA binding and transcription in Drosophila

The Drosophila transcription factor Bicoid (Bcd) binds thousands of genomic sites during early embryogenesis, but it is unclear how many of these binding events are functionally important. Here, Small and colleagues test the role of the maternal factor Zelda (Zld) in Bcd-mediated binding and transcription. Embryos lacking Zld show enhanced Bcd binding to a subset of genomic locations, causing early activation of target genes normally silent until later stages. This study demonstrates a critical role for Zld in controlling Bcd binding and target gene activation in the early embryo.

In vivo cross-linking studies suggest that the Drosophila transcription factor Bicoid (Bcd) binds to several thousand sites during early embryogenesis, but it is not clear how many of these binding events are functionally important. In contrast, reporter gene studies have identified >60 Bcd-dependent enhancers, all of which contain clusters of the consensus binding sequence TAATCC. These studies also identified clusters of TAATCC motifs (inactive fragments) that failed to drive Bcd-dependent activation. In general, active fragments showed higher levels of Bcd binding in vivo and were enriched in predicted binding sites for the ubiquitous maternal protein Zelda (Zld). Here we tested the role of Zld in Bcd-mediated binding and transcription. Removal of Zld function and mutations in Zld sites caused significant reductions in Bcd binding to known enhancers and variable effects on the activation and spatial positioning of Bcd-dependent expression patterns. Also, insertion of Zld sites converted one of six inactive fragments into a Bcd-responsive enhancer. Genome-wide binding experiments in zld mutants showed variable effects on Bcd-binding peaks, ranging from strong reductions to significantly enhanced levels of binding. Increases in Bcd binding caused the precocious Bcd-dependent activation of genes that are normally not expressed in early embryos, suggesting that Zld controls the genome-wide binding profile of Bcd at the qualitative level and is critical for selecting target genes for activation in the early embryo. These results underscore the importance of combinatorial binding in enhancer function and provide data that will help predict regulatory activities based on DNA sequence.

Distinct molecular strategies for Hox-mediated limb suppression in Drosophila: from cooperativity to dispensability/antagonism in TALE partnership

The emergence following gene duplication of a large repertoire of Hox paralogue proteins underlies the importance taken by Hox proteins in controlling animal body plans in development and evolution. Sequence divergence of paralogous proteins accounts for functional specialization, promoting axial morphological diversification in bilaterian animals. Yet functionally specialized paralogous Hox proteins also continue performing ancient common functions. In this study, we investigate how highly divergent Hox proteins perform an identical function. This was achieved by comparing in Drosophila the mode of limb suppression by the central (Ultrabithorax and AbdominalA) and posterior class (AbdominalB) Hox proteins. Results highlight that Hox-mediated limb suppression relies on distinct modes of DNA binding and a distinct use of TALE cofactors. Control of common functions by divergent Hox proteins, at least in the case studied, relies on evolving novel molecular properties. Thus, changes in protein sequences not only provide the driving force for functional specialization of Hox paralogue proteins, but also provide means to perform common ancient functions in distinct ways.

BmTGIF, a Bombyx mori homolog of Drosophila DmTGIF, regulates progression of spermatogenesis

TG-interacting factor (TGIF) in Drosophila consists of two tandemly-repeated genes, achintya (Dmachi) and vismay (Dmvis), which act as transcriptional activators in Drosophila spermatogenesis. In contrast, TGIF in humans is a transcriptional repressor that binds directly to DNA or interacts with corepressors to repress the transcription of target genes. In this study, we investigated the characteristics and functions of BmTGIF, a Bombyx mori homolog of DmTGIF. Like DmTGIF, BmTGIF is predominantly expressed in the testes and ovaries. Four alternatively spliced isoforms could be isolated from testes, and two isoforms from ovaries. Quantitative polymerase chain reaction indicated BmTGIF was abundantly expressed in the testis of 3rd instar larvae, when the testis is almost full of primary spermatocytes. The results of luciferase assays indicated that BmTGIF contains two adjacent acidic domains that activate the transcription of reporter genes. Immunofluorescence assay in BmN cells showed that the BmTGIF protein was located mainly in the nucleus, and paraffin sections of testis showed BmTGIF was grossly expressed in primary spermatocytes and mature sperms. Consistent with the role of DmVis in Drosophila development, BmTGIF significantly affected spermatid differentiation, as indicated by hematoxylin-eosin staining of paraffin sections of testis from BmTGIF-small interfering RNA (siRNA)-injected male silkworms. Co-immunoprecipitation experiments suggested that BmTGIF interacted with BmAly, and that they may recruit other factors to form a complex to regulate the genes required for meiotic divisions and spermatid differentiation. The results of this analysis of BmTGIF will improve our understanding of the mechanism of spermatid differentiation in B. mori, with potential applications for pest control.

The N-terminus of the human RecQL4 helicase is a homeodomain-like DNA interaction motif

The RecQL4 helicase is involved in the maintenance of genome integrity and DNA replication. Mutations in the human RecQL4 gene cause the Rothmund–Thomson, RAPADILINO and Baller–Gerold syndromes. Mouse models and experiments in human and Xenopus have proven the N-terminal part of RecQL4 to be vital for cell growth. We have identified the first 54 amino acids of RecQL4 (RecQL4_N54) as the minimum interaction region with human TopBP1. The solution structure of RecQL4_N54 was determined by heteronuclear liquid–state nuclear magnetic resonance (NMR) spectroscopy (PDB 2KMU; backbone root-mean-square deviation 0.73 Å). Despite low-sequence homology, the well-defined structure carries an overall helical fold similar to homeodomain DNA-binding proteins but lacks their archetypical, minor groove-binding N-terminal extension. Sequence comparison indicates that this N-terminal homeodomain-like fold is a common hallmark of metazoan RecQL4 and yeast Sld2 DNA replication initiation factors. RecQL4_N54 binds DNA without noticeable sequence specificity yet with apparent preference for branched over double-stranded (ds) or single-stranded (ss) DNA. NMR chemical shift perturbation observed upon titration with Y-shaped, ssDNA and dsDNA shows a major contribution of helix α3 to DNA binding, and additional arginine side chain interactions for the ss and Y-shaped DNA.

Exploring the DNA-recognition potential of homeodomains

The recognition potential of most families of DNA-binding domains (DBDs) remains relatively unexplored. Homeodomains (HDs), like many other families of DBDs, display limited diversity in their preferred recognition sequences. To explore the recognition potential of HDs, we utilized a bacterial selection system to isolate HD variants, from a randomized library, that are compatible with each of the 64 possible 3' triplet sites (i.e., TAANNN). The majority of these selections yielded sets of HDs with overrepresented residues at specific recognition positions, implying the selection of specific binders. The DNA-binding specificity of 151 representative HD variants was subsequently characterized, identifying HDs that preferentially recognize 44 of these target sites. Many of these variants contain novel combinations of specificity determinants that are uncommon or absent in extant HDs. These novel determinants, when grafted into different HD backbones, produce a corresponding alteration in specificity. This information was used to create more explicit HD recognition models, which can inform the prediction of transcriptional regulatory networks for extant HDs or the engineering of HDs with novel DNA-recognition potential. The diversity of recovered HD recognition sequences raises important questions about the fitness barrier that restricts the evolution of alternate recognition modalities in natural systems.

Cofactor binding evokes latent differences in DNA binding specificity between Hox proteins

Members of transcription factor families typically have similar DNA binding specificities yet execute unique functions in vivo. Transcription factors often bind DNA as multiprotein complexes, raising the possibility that complex formation might modify their DNA binding specificities. To test this hypothesis, we developed an experimental and computational platform, SELEX-seq, that can be used to determine the relative affinities to any DNA sequence for any transcription factor complex. Applying this method to all eight Drosophila Hox proteins, we show that they obtain novel recognition properties when they bind DNA with the dimeric cofactor Extradenticle-Homothorax (Exd). Exd-Hox specificities group into three main classes that obey Hox gene collinearity rules and DNA structure predictions suggest that anterior and posterior Hox proteins prefer DNA sequences with distinct minor groove topographies. Together, these data suggest that emergent DNA recognition properties revealed by interactions with cofactors contribute to transcription factor specificities in vivo.

Comparisons of the embryonic development of Drosophila, Nasonia, and Tribolium

Studying the embryogenesis of diverse insect species is crucial to understanding insect evolution. Here, we review current advances in understanding the development of two emerging model organisms: the wasp Nasonia vitripennis and the beetle Tribolium castaneum in comparison with the well-studied fruit fly Drosophila melanogaster. Although Nasonia represents the most basally branching order of holometabolous insects, it employs a derived long germband mode of embryogenesis, more like that of Drosophila, whereas Tribolium undergoes an intermediate germband mode of embryogenesis, which is more similar to the ancestral mechanism. Comparing the embryonic development and genetic regulation of early patterning events in these three insects has given invaluable insights into insect evolution. The similar mode of embryogenesis of Drosophila and Nasonia is reflected in their reliance on maternal morphogenetic gradients. However, they employ different genes as maternal factors, reflecting the evolutionary distance separating them. Tribolium, on the other hand, relies heavily on self-regulatory mechanisms other than maternal cues, reflecting its sequential nature of segmentation and the need for reiterated patterning.

Interlocked feedforward loops control cell-type-specific Rhodopsin expression in the Drosophila eye

How complex networks of activators and repressors lead to exquisitely specific cell-type determination during development is poorly understood. In the Drosophila eye, expression patterns of Rhodopsins define at least eight functionally distinct though related subtypes of photoreceptors. Here, we describe a role for the transcription factor gene defective proventriculus (dve) as a critical node in the network regulating Rhodopsin expression. dve is a shared component of two opposing, interlocked feedforward loops (FFLs). Orthodenticle and Dve interact in an incoherent FFL to repress Rhodopsin expression throughout the eye. In R7 and R8 photoreceptors, a coherent FFL relieves repression by Dve while activating Rhodopsin expression. Therefore, this network uses repression to restrict and combinatorial activation to induce cell-type-specific expression. Furthermore, Dve levels are finely tuned to yield cell-type- and region-specific repression or activation outcomes. This interlocked FFL motif may be a general mechanism to control terminal cell-fate specification.

Structural basis for sequence specific DNA binding and protein dimerization of HOXA13

The homeobox gene (HOXA13) codes for a transcription factor protein that binds to AT-rich DNA sequences and controls expression of genes during embryonic morphogenesis. Here we present the NMR structure of HOXA13 homeodomain (A13DBD) bound to an 11-mer DNA duplex. A13DBD forms a dimer that binds to DNA with a dissociation constant of 7.5 nM. The A13DBD/DNA complex has a molar mass of 35 kDa consistent with two molecules of DNA bound at both ends of the A13DBD dimer. A13DBD contains an N-terminal arm (residues 324 – 329) that binds in the DNA minor groove, and a C-terminal helix (residues 362 – 382) that contacts the ATAA nucleotide sequence in the major groove. The N370 side-chain forms hydrogen bonds with the purine base of A5* (base paired with T5). Side-chain methyl groups of V373 form hydrophobic contacts with the pyrimidine methyl groups of T5, T6* and T7*, responsible for recognition of TAA in the DNA core. I366 makes similar methyl contacts with T3* and T4*. Mutants (I366A, N370A and V373G) all have decreased DNA binding and transcriptional activity. Exposed protein residues (R337, K343, and F344) make intermolecular contacts at the protein dimer interface. The mutation F344A weakens protein dimerization and lowers transcriptional activity by 76%. We conclude that the non-conserved residue, V373 is critical for structurally recognizing TAA in the major groove, and that HOXA13 dimerization is required to activate transcription of target genes.

Systematic characterization of protein-DNA interactions

Sequence-specific protein-DNA interactions (PDIs) are critical for regulating many cellular processes, including transcription, DNA replication, repair, and rearrangement. We review recent experimental advances in high-throughput technologies designed to characterize PDIs and discuss recent studies that use these tools, including ChIP-chip/seq, SELEX-based approaches, yeast one-hybrid, bacterial one-hybrid, protein binding microarray, and protein microarray. The results of these studies have challenged some long-standing concepts of PDI and provide valuable insights into the complex transcriptional regulatory networks.

Using the recognition code to swap homeodomain target specificity in cell culture

The homeodomain (HD) is a 60 amino acid-long DNA-binding domain. A large fraction of HDs binds with high affinity sequences containing the 5'-TAAT-3' core motif. However, NK-2 class HDs recognizes sequences containing the 5'-CAAG-3' core motif. By using a cell transfection approach, here we show that modification of residues located in the N-terminal arm (at positions 6, 7 and 8) and in the recognition helix (at position 54) is enough to swap the "in vivo" binding specificity of TTF-1 HD (which is a member of the NK-2 class HD) from 5'-CAAG-3' to 5'-TAAT-3'-containing targets. The role of residue at position 54 is also supported by data obtained with the HD of the Drosophila engrailed protein. These data support the notion that DNA-binding specificity "in vivo" is dictated by few critical residues.

Origins of specificity in protein-DNA recognition

Specific interactions between proteins and DNA are fundamental to many biological processes. In this review, we provide a revised view of protein-DNA interactions that emphasizes the importance of the three-dimensional structures of both macromolecules. We divide protein-DNA interactions into two categories: those when the protein recognizes the unique chemical signatures of the DNA bases (base readout) and those when the protein recognizes a sequence-dependent DNA shape (shape readout). We further divide base readout into those interactions that occur in the major groove from those that occur in the minor groove. Analogously, the readout of the DNA shape is subdivided into global shape recognition (for example, when the DNA helix exhibits an overall bend) and local shape recognition (for example, when a base pair step is kinked or a region of the minor groove is narrow). Based on the >1500 structures of protein-DNA complexes now available in the Protein Data Bank, we argue that individual DNA-binding proteins combine multiple readout mechanisms to achieve DNA-binding specificity. Specificity that distinguishes between families frequently involves base readout in the major groove, whereas shape readout is often exploited for higher resolution specificity, to distinguish between members within the same DNA-binding protein family.

Autoregulatory binding sites in the zebrafish six3a promoter region define a new recognition sequence for Six3 proteins

The homeodomain (HD) transcription factor Six3, which is a member of the Six/Sine oculis family, is essential for development of the eyes and forebrain in vertebrates. It has recently been claimed that the HDs of Six3 and other members of the Six family have a common recognition sequence, TGATAC. However, a different recognition sequence including the typical TAAT core motif, which has not yet been fully defined, has also been proposed for the Six3 HD in mice. Our study of the zebrafish orthologue six3a, which has an identical HD, shows that it binds in vitro to multiple TAAT-containing sites within its promoter region. Comparison of the different binding affinities for these sequences identifies three high-affinity sites with a common TAATGTC motif. Notably, this new recognition sequence, which is supported by our analysis of the influence of single-nucleotide substitutions on the DNA-binding affinity, is distinct from all of the DNA-binding specificities previously described in surveys of HDs. In addition, our comparison of Six3a HD binding to the novel TAATGTC motif and the common recognition sequence of Six family HDs (TGATAC) shows very similar affinities, suggesting two distinct DNA-binding modes. Transient reporter assays of the six3a promoter in zebrafish embryos also indicate that the three high-affinity sites are involved in autoregulation. In support of this, chromatin immunoprecipitation experiments show enrichment of Six3a binding to a six3a promoter fragment containing two clustered high-affinity sites. These findings provide strong evidence that the TAATGTC motif is an important target sequence for vertebrate Six3 proteins in vivo.