Determination of the nuclear magnetic resonance solution structure of an Antennapedia homeodomain-DNA complex

Transcriptional precision in photoreceptor development and diseases - Lessons from 25 years of CRX research

The vertebrate retina is made up of six specialized neuronal cell types and one glia that are generated from a common retinal progenitor. The development of these distinct cell types is programmed by transcription factors that regulate the expression of specific genes essential for cell fate specification and differentiation. Because of the complex nature of transcriptional regulation, understanding transcription factor functions in development and disease is challenging. Research on the Cone-rod homeobox transcription factor CRX provides an excellent model to address these challenges. In this review, we reflect on 25 years of mammalian CRX research and discuss recent progress in elucidating the distinct pathogenic mechanisms of four CRX coding variant classes. We highlight how in vitro biochemical studies of CRX protein functions facilitate understanding CRX regulatory principles in animal models. We conclude with a brief discussion of the emerging systems biology approaches that could accelerate precision medicine for CRX-linked diseases and beyond.

Mammal Reproductive Homeobox (Rhox) Genes: An Update of Their Involvement in Reproduction and Development

The reproductive homeobox on the X chromosome (RHOX) genes were first identified in the mouse during the 1990s and have a crucial role in reproduction. In various transcription factors with a key regulatory role, the homeobox sequence encodes a "homeodomain" DNA-binding motif. In the mouse, there are three clusters of Rhox genes (α, β, and γ) on the X chromosome. Each cluster shows temporal and/or quantitative collinearity, which regulates the progression of the embryonic development process. Although the RHOX family is conserved in mammals, the interspecies differences in the number of RHOX genes and pseudogenes testifies to a rich evolutionary history with several relatively recent events. In the mouse, Rhox genes are mainly expressed in reproductive tissues, and several have a role in the differentiation of primordial germ cells (Rhox1, Rhox6, and Rhox10) and in spermatogenesis (Rhox1, Rhox8, and Rhox13). Despite the lack of detailed data on human RHOX, these genes appear to be involved in the formation of germ cells because they are predominantly expressed during the early (RHOXF1) and late (RHOXF2/F2B) stages of germ cell development. Furthermore, the few variants identified to date are thought to induce or predispose to impaired spermatogenesis and severe oligozoospermia or azoospermia. In the future, research on the pathophysiology of the human RHOX genes is likely to confirm the essential role of this family in the reproductive process and might help us to better understand the various causes of infertility and characterize the associated human phenotypes.

KNOX Genes Were Involved in Regulating Axillary Bud Formation of Chrysanthemum × morifolium

Branching is an important agronomic and economic trait in cut chrysanthemums. The axillary meristem (AM) formation of the axillary buds of cut chrysanthemums has a decisive role in its branching characteristics. However, little is known about the regulation mechanism of axillary meristem formation in chrysanthemums at the molecular level. Members of the Homeobox gene family especially genes belonging to the class I KNOX branch play a key role in regulating the axillary bud growth and development processes of plants. In this study, three genes belonging to the class I KNOX branch, CmKNAT1, CmKNAT6, and CmSTM were cloned from chrysanthemums, and their functions in regulating axillary bud formation were examined. The subcellular localization test showed that these three KNOX genes were expressed in the nucleus, so all of them might function as transcription factors. The results of the expression profile analysis showed that these three KNOX genes were highly expressed in the AM formation stage of axillary buds. Overexpression of KNOX genes result in a wrinkled leaf phenotype in tobacco and Arabidopsis, which may be related to the excessive division of leaf cells, resulting in the proliferation of leaf tissue. Furthermore, overexpression of these three KNOX genes enhances the regeneration ability of tobacco leaves, indicating that these three KNOX genes may participate in the regulation of cell meristematic ability, thus promoting the formation of buds. In addition, the results of fluorescence quantitative testing showed that these three KNOX genes may promote the formation of chrysanthemum axillary buds by promoting the cytokinin pathway while inhibiting the auxin and gibberellin pathways. In conclusion, this study demonstrated that CmKNAT1, CmKNAT6, and CmSTM genes were involved in regulating axillary bud formation of Chrysanthemum × morifolium and preliminarily revealed the molecular mechanism of their regulation of AM formation. These findings may provide a theoretical basis and candidate gene resources for genetic engineering breeding of new varieties of cut chrysanthemums without lateral branches.

Genome-wide identification of KNOX transcription factors in cotton and the role of GhKNOX4-A and GhKNOX22-D in response to salt and drought stress

Cotton is one of the most important economic and fiber crops in the world. KNOX is one class of universal transcription factors, which plays important roles in plant growth and development as well as response to different stresses. Although there are many researches on KNOXs in other plant species, there are few reports on cotton. In this study, we systematically and comprehensively identified all KNOX genes in upland cotton and its two ancestral species; we also studied their functions by employing RNA-seq analysis and virus-induced gene silence (VIGS). A total of 89 KNOX genes were identified from three cotton species. Among them, 44 were from upland cotton, 22 and 23 were found in its ancestral species G. raimondii and G. arboreum, respectively. Plant polyploidization and domestication play a selective force driving KNOX gene evolution. Phylogenetic analysis displayed that KNOX genes were evolved into three Classes. The intron length and exon number differed in each Class. Transcriptome data showed that KNOX genes of Class II were widely expressed in multiple tissues, including fiber. The majority of KNOX genes were induced by different abiotic stresses. Additionally, we found multiple cis-elements related to stress in the promoter region of KNOX genes. VIGS silence of GhKNOX4-A and GhKNOX22-D genes showed significant growth and development effect in cotton seedlings under salt and drought treatments. Both GhKNOX4-A and GhKNOX22-D regulated plant tolerance; silencing both genes induced oxidative stresses, evidenced by reduced SOD activity and induced leave cell death, and also enhanced stomatal open and water loss. Thus, GhKNOX4-A and GhKNOX22-D may contribute to drought response by regulating stomata opening and oxidative stresses.

Biophysical insights into glucose-dependent transcriptional regulation by PDX1

The pancreatic and duodenal homeobox 1 (PDX1) is a central regulator of glucose-dependent transcription of insulin in pancreatic β cells. PDX1 transcription factor activity is integral to the development and sustained health of the pancreas; accordingly, deciphering the complex network of cellular cues that lead to PDX1 activation or inactivation is an important step toward understanding the etiopathologies of pancreatic diseases and the development of novel therapeutics. Despite nearly 3 decades of research into PDX1 control of Insulin expression, the molecular mechanisms that dictate the function of PDX1 in response to glucose are still elusive. The transcriptional activation functions of PDX1 are regulated, in part, by its two intrinsically disordered regions, which pose a barrier to its structural and biophysical characterization. Indeed, many studies of PDX1 interactions, clinical mutations, and posttranslational modifications lack molecular level detail. Emerging methods for the quantitative study of intrinsically disordered regions and refined models for transactivation now enable us to validate and interrogate the biochemical and biophysical features of PDX1 that dictate its function. The goal of this review is to summarize existing PDX1 studies and, further, to generate a comprehensive resource for future studies of transcriptional control via PDX1.

Genome-wide identification of three-amino-acid-loop-extension gene family and their expression profile under hormone and abiotic stress treatments during stem development of Prunus mume

Transcription factors encoded by the three-amino-acid-loop-extension (TALE) gene family play a key role in regulating plant growth and development, and are involved in plant hormone regulatory pathways and responses to various environmental stresses. Researchers are currently studying TALE genes in different species, but Prunus mume TALE genes have not yet been studied. Therefore, based on the P. mume genome, we found a total of 23 TALE gene family members, which were distributed on eight chromosomes. TALE genes contained the characteristic domains of this family, and could be divided into KNOTTED-like homeobox (KNOX) subfamily and BEL1-like homeobox (BELL) subfamily. They can form heterodimers with each other. Fragment duplication and tandem duplication events were the main reasons for the expansion of P. mume TALE gene family members and the TALE genes were selected by different degrees of purification. The inter-species collinearity analysis showed that the relationship between P. mume and other four Prunus species was consistent with the distance of origin. Eleven members of P. mume TALE genes were specifically highly expressed in stem, mainly at the early stage of stem development. The cis-element analysis showed that the promoter of P. mume TALE genes contained a variety of hormone and abiotic stress response elements, and four TALE genes responded to two kinds of abiotic stresses and four kinds of hormones at the early stage of stem development. In conclusion, this study lays a foundation to explore the role of TALE gene family in P. mume growth and development.

Identification of TALE Transcription Factor Family and Expression Patterns Related to Fruit Chloroplast Development in Tomato (Solanum lycopersicum L.)

The TALE gene family is an important transcription factor family that regulates meristem formation, organ morphogenesis, signal transduction, and fruit development. A total of 24 genes of the TALE family were identified and analyzed in tomato. The 24 SlTALE family members could be classified into five BELL subfamilies and four KNOX subfamilies. SlTALE genes were unevenly distributed on every tomato chromosome, lacked syntenic gene pairs, and had conserved structures but diverse regulatory functions. Promoter activity analysis showed that cis-elements responsive to light, phytohormone, developmental regulation, and environmental stress were enriched in the promoter of SlTALE genes, and the light response elements were the most abundant. An abundance of TF binding sites was also enriched in the promoter of SlTALE genes. Phenotype identification revealed that the green shoulder (GS) mutant fruits showed significantly enhanced chloroplast development and chlorophyll accumulation, and a significant increase of chlorophyll fluorescence parameters in the fruit shoulder region. Analysis of gene expression patterns indicated that six SlTALE genes were highly expressed in the GS fruit shoulder region, and four SlTALE genes were highly expressed in the parts with less-developed chloroplasts. The protein-protein interaction networks predicted interaction combinations among these SlTALE genes, especially between the BELL subfamilies and the KNOX subfamilies, indicating a complex regulatory network of these SlTALE genes in chloroplast development and green fruit shoulder formation. In conclusion, our result provides detailed knowledge of the SlTALE gene for functional research and the utilization of the TALE gene family in fruit quality improvement.

Genome-Wide Identification, Expression, and Interaction Analysis of BEL-Like Homeodomain Gene Family in Peach

BEL1-like homeodomain (BLH) family genes as homeodomain transcription factors are found ubiquitously in plants to play important regulatory roles in reproductive development, morphological development, and stress response. Although BLH proteins have been reported in some species, there is little information about BLH genes in peach. In this study, we identified 11 peach PpBLH genes based on the conserved domain. Phylogenetic analysis suggested that the PpBLH proteins could be divided into five groups, which might be involved in different aspects of morphogenesis. Genomics structure analysis revealed that there were four exons in the PpBLH gene, and the length of the third exon was 61 bp. Chromosomal location analysis showed that the PpBLH genes were not distributed uniformly on six chromosomes. Promoter analysis showed that the promoter sequences of six PpBLH genes contained multiple cis-acting elements for hormones and stress. Six PpBLH genes were cloned by RT-PCR, and PpBLH1, PpBLH4, and PpBLH7 showed different expression patterns in the tested fruits under common temperature and high temperature. Y2H results indicated that PpBLH7 andPpBLH10 interacted with the PpOFP6 protein, and PpBLH1 interacted with the PpOFP1, PpOFP2, PpOFP4, and PpOFP13 proteins. These results provide new insight for further study of PpBLH genes, and construction of regulatory networks of PpBLH proteins in the growth, development, and stress response of peach.

Cognate DNA Recognition by Engrailed Homeodomain Involves a Conformational Change Controlled via an Electrostatic-Spring-Loaded Latch

Transcription factors must scan genomic DNA, recognize the cognate sequence of their control element(s), and bind tightly to them. The DNA recognition process is primarily carried out by their DNA binding domains (DBD), which interact with the cognate site with high affinity and more weakly with any other DNA sequence. DBDs are generally thought to bind to their cognate DNA without changing conformation (lock-and-key). Here, we used nuclear magnetic resonance and circular dichroism to investigate the interplay between DNA recognition and DBD conformation in the engrailed homeodomain (enHD), as a model case for the homeodomain family of eukaryotic DBDs. We found that the conformational ensemble of enHD is rather flexible and becomes gradually more disordered as ionic strength decreases following a Debye–Hückel’s dependence. Our analysis indicates that enHD’s response to ionic strength is mediated by a built-in electrostatic spring-loaded latch that operates as a conformational transducer. We also found that, at moderate ionic strengths, enHD changes conformation upon binding to cognate DNA. This change is of larger amplitude and somewhat orthogonal to the response to ionic strength. As a consequence, very high ionic strengths (e.g., 700 mM) block the electrostatic-spring-loaded latch and binding to cognate DNA becomes lock-and-key. However, the interplay between enHD conformation and cognate DNA binding is robust across a range of ionic strengths (i.e., 45 to 300 mM) that covers the physiologically-relevant conditions. Therefore, our results demonstrate the presence of a mechanism for the conformational control of cognate DNA recognition on a eukaryotic DBD. This mechanism can function as a signal transducer that locks the DBD in place upon encountering the cognate site during active DNA scanning. The electrostatic-spring-loaded latch of enHD can also enable the fine control of DNA recognition in response to transient changes in local ionic strength induced by variate physiological processes.

Genome-wide analysis of TALE superfamily in Triticum aestivum reveals TaKNOX11-A is involved in abiotic stress response

BACKGROUND

Three-amino-loop-extension (TALE) superfamily genes are widely present in plants and function directly in plant growth and development and abiotic stress response. Although TALE genes have been studied in many plant species, members of the TALE family have not been identified in wheat.

RESULTS

In this study, we identified 70 wheat TALE protein candidate genes divided into two subfamilies, KNOX (KNOTTED-like homeodomain) and BEL1-like (BLH/BELL homeodomain). Genes in the same subfamily or branch in the phylogenetic tree are similar in structure, and their encoded proteins have similar motifs and conserved structures. Wheat TALE genes are unevenly distributed on 21 chromosomes and expanded on the fourth chromosome. Through gene duplication analysis, 53 pairs of wheat TALE genes were determined to result from segmental duplication events, and five pairs were caused by tandem duplication events. The Ka/Ks between TALE gene pairs indicates a strong purification and selection effect. There are multiple cis-elements in the 2000 bp promoter sequence that respond to hormones and abiotic stress, indicating that most wheat TALE genes are involved in the growth, development, and stress response of wheat. We also studied the expression profiles of wheat TALE genes in different developmental stages and tissues and under different stress treatments. We detected the expression levels of four TALE genes by qRT-PCR, and selected TaKNOX11-A for further downstream analysis. TaKNOX11-A enhanced the drought and salt tolerances of Arabidopsis thaliana. TaKNOX11-A overexpressing plants had decreased malondialdehyde content and increased proline content, allowing for more effective adaptation of plants to unfavorable environments.

CONCLUSIONS

We identified TALE superfamily members in wheat and conducted a comprehensive bioinformatics analysis. The discovery of the potential role of TaKNOX11-A in drought resistance and salt tolerance provides a basis for follow-up studies of wheat TALE family members, and also provides new genetic resources for improving the stress resistance of wheat.

Racemic phosphorothioate as a tool for NMR investigations of protein-DNA complexes

A major driving force for protein-nucleic acid association is electrostatic interactions via ion pairs of the positively charged basic side chains and negatively charged phosphates. For a better understanding of how proteins scan DNA and recognize particular signatures, it is important to gain atomic-level insight into the behavior of basic side chains at the protein-DNA interfaces. NMR spectroscopy is a powerful tool for investigating the structural, dynamic, and kinetic aspects of protein-DNA interactions. However, resonance assignment of basic side-chain cationic moieties at the molecular interfaces remains to be a major challenge. Here, we propose a fast, robust, and inexpensive approach that greatly facilitates resonance assignment of interfacial moieties and also allows for kinetic measurements of protein translocation between two DNA duplexes. This approach utilizes site-specific incorporation of racemic phosphorothioate at the position of a phosphate that interacts with a protein side chain. This modification retains the electric charge of phosphate and therefore is mild, but causes significant chemical shift perturbations for the proximal protein side chains, which facilitates resonance assignment. Due to the racemic nature of the modification, two different chemical shifts are observed for the species with different diastereomers R_P and S_P of the incorporated phosphorothioate group. Kinetic information on the exchange of the protein molecule between R_P and S_P DNA duplexes can be obtained by ¹⁵N_z exchange spectroscopy. We demonstrate the applications of this approach to the Antennapedia homeodomain-DNA complex and the CREB1 basic leucine-zipper (bZIP)-DNA complex.

Protein isoelectric point distribution in the interactomes across the domains of life

The distribution of the protein isoelectric point (pI) in the protein-protein interaction (PPI) networks across the domains of life has not been investigated yet. This work attempts to correlate the pI with the number of direct interacting partners in the experimentally supported networks involving 226.085 PPIs from 14 various organisms including human, mouse, yeast, bacteria, viruses and 53.606 virus-host interactions. The results showed that the acidic proteins (pI<3) have the highest average number of interactions in eukaryotes, while in bacteria more neutral proteins. On the contrary, the basic proteins (pI>11) have the lowest average number of interactions in human, mouse, yeast, bacteria and human-viral interactomes and the highest average in intraviral interactomes. We examined the correlation of the pI of the interacting partners by calculating the assortativity index of various PPI networks. We found that the interactions between the acidic, neutral and basic proteins have a fairly random mix, implying weak if any association between the acidic and basic proteins. Furthermore, protein features such as biological function, structurally order and disorder, subcellular localization, and homodimerization were classified according to pI in prokaryote and eukaryote proteomes.

Comprehensive analysis of the three-amino-acid-loop-extension gene family and its tissue-differential expression in response to salt stress in poplar

The three-amino-acid-loop-extension (TALE) transcription factor gene family is widely present in plants and plays an important role in its growth and development. However, studies on the gene family are limited in poplar. In this study, we investigated 35 TALE gene family members in terms of their evolutionary relationship, classification, physicochemical properties, gene structures, and protein motifs. We divided the genes into four classes, based on their protein sequences similarity. The members from each class share similar gene structures and motif compositions. Evidence from transcript profiling indicated that the majority of the TALE genes exhibited distinct expression patterns over leaf, stem, and root tissues. Out of the 35 genes, 17 genes are highly expressed in stems, suggesting that the TALE gene family may play an important role in secondary growth and wood formation. Furthermore, out of the 35 genes, 11 genes are responsive to salt stress, and the spatio-temporal expression patterns of these 11 genes under salt stress were analysed using RT-qPCR. Yeast two-hybridization analysis indicated that poplar TALE proteins from different classes can form heterodimers. These results lay the foundation for future studies on biological functions of poplar TALE genes.

Rapid and efficient purification of Drosophila homeodomain transcription factors for biophysical characterization

Homeodomain transcription factors (HD TFs) are a large class of evolutionarily conserved DNA binding proteins that contain a basic 60-amino acid region required for binding to specific DNA sites. In Drosophila melanogaster, many of these HD TFs are expressed in the early embryo and control transcription of target genes in development through their interaction with cis-regulatory modules. Previous studies where some of the Drosophila HD TFs were purified required the use of strong denaturants (i.e. 6 M urea) and multiple chromatography columns, making the downstream biochemical examination of the isolated protein difficult. To circumvent these obstacles, we have developed a streamlined expression and purification protocol to produce large yields of Drosophila HD TFs. Using the HD TFs FUSHI-TARAZU (FTZ), ANTENNAPEDIA (ANTP), ABDOMINAL-A (ABD-A), ABDOMINAL-B (ABD-B), and ULTRABITHORAX (UBX) as examples, we demonstrate that our 3-day protocol involving the overexpression of His₆-SUMO fusion constructs in E. coli followed by a Ni²⁺-IMAC, SUMO-tag cleavage with the SUMO protease Ulp1, and a heparin column purification produces pure, soluble protein in biological buffers around pH 7 in the absence of denaturants. Electrophoretic mobility shift assays (EMSA) confirm that the purified HD proteins are functional and nuclear magnetic resonance (NMR) spectra confirm that the purified HDs are well-folded. These purified HD TFs can be used in future biophysical experiments to structurally and biochemically characterize how and why these HD TFs bind to different DNA sequences and further probe how nucleotide differences contribute to TF-DNA specificity in the HD family.

Evolutionary Potential of Cis-Regulatory Mutations to Cause Rapid Changes in Transcription Factor Binding

Transcriptional regulation is crucial for all biological processes and well investigated at the molecular level for a wide range of organisms. However, it is quite unclear how innovations, such as the activity of a novel regulatory element, evolve. In the case of transcription factor (TF) binding, both a novel TF and a novel-binding site would need to evolve concertedly. Since promiscuous functions have recently been identified as important intermediate steps in creating novel specific functions in many areas such as enzyme evolution and protein-protein interactions, we ask here how promiscuous binding of TFs to TF-binding sites (TFBSs) affects the robustness and evolvability of this tightly regulated system. Specifically, we investigate the binding behavior of several hundred TFs from different species at unprecedented breadth. Our results illustrate multiple aspects of TF-binding interactions, ranging from correlations between the strength of the interaction bond and specificity, to preferences regarding TFBS nucleotide composition in relation to both domains and binding specificity. We identified a subset of high A/T binding motifs. Motifs in this subset had many functionally neutral one-error mutants, and were bound by multiple different binding domains. Our results indicate that, especially for some TF-TFBS associations, low binding specificity confers high degrees of evolvability, that is that few mutations facilitate rapid changes in transcriptional regulation, in particular for large and old TF families. In this study we identify binding motifs exhibiting behavior indicating high evolutionary potential for innovations in transcriptional regulation.

Common ancestry of heterodimerizing TALE homeobox transcription factors across Metazoa and Archaeplastida

Background

Complex multicellularity requires elaborate developmental mechanisms, often based on the versatility of heterodimeric transcription factor (TF) interactions. Homeobox TFs in the TALE superclass are deeply embedded in the gene regulatory networks that orchestrate embryogenesis. Knotted-like homeobox (KNOX) TFs, homologous to animal MEIS, have been found to drive the haploid-to-diploid transition in both unicellular green algae and land plants via heterodimerization with other TALE superclass TFs, demonstrating remarkable functional conservation of a developmental TF across lineages that diverged one billion years ago. Here, we sought to delineate whether TALE-TALE heterodimerization is ancestral to eukaryotes.

Results

We analyzed TALE endowment in the algal radiations of Archaeplastida, ancestral to land plants. Homeodomain phylogeny and bioinformatics analysis partitioned TALEs into two broad groups, KNOX and non-KNOX. Each group shares previously defined heterodimerization domains, plant KNOX-homology in the KNOX group and animal PBC-homology in the non-KNOX group, indicating their deep ancestry. Protein-protein interaction experiments showed that the TALEs in the two groups all participated in heterodimerization.

Conclusions

Our study indicates that the TF dyads consisting of KNOX/MEIS and PBC-containing TALEs must have evolved early in eukaryotic evolution. Based on our results, we hypothesize that in early eukaryotes, the TALE heterodimeric configuration provided transcription-on switches via dimerization-dependent subcellular localization, ensuring execution of the haploid-to-diploid transition only when the gamete fusion is correctly executed between appropriate partner gametes. The TALE switch then diversified in the several lineages that engage in a complex multicellular organization.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0605-5) contains supplementary material, which is available to authorized users.

Molecular dynamics simulation on the allosteric analysis of the c-di-GMP class I riboswitch induced by ligand binding

Riboswitches are RNA molecules that regulate gene expression using conformation change, affected by binding of small molecule ligands. Although a number of ligand-bound aptamer complex structures have been solved, it is important to know ligand-free conformations of the aptamers in order to understand the mechanism of specific binding by ligands. In this paper, we use dynamics simulations on a series of models to characterize the ligand-free and ligand-bound aptamer domain of the c-di-GMP class I (GEMM-I) riboswitch. The results revealed that the ligand-free aptamer has a stable state with a folded P2 and P3 helix, an unfolded P1 helix and open binding pocket. The first Mg ions binding to the aptamer is structurally favorable for the successive c-di-GMP binding. The P1 helix forms when c-di-GMP is successive bound. Three key junctions J1/2, J2/3 and J1/3 in the GEMM-I riboswitch contributing to the formation of P1 helix have been found. The binding of the c-di-GMP ligand to the GEMM-I riboswitch induces the riboswitch's regulation through the direct allosteric communication network in GEMM-I riboswitch from the c-di-GMP binding sites in the J1/2 and J1/3 junctions to the P1 helix, the indirect ones from those in the J2/3 and P2 communicating to P1 helix via the J1/2 and J1/3 media.

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.

Internal Motions of Basic Side Chains of the Antennapedia Homeodomain in the Free and DNA-Bound States

Basic side chains play crucial roles in protein-DNA interactions. In this study, using NMR spectroscopy, we investigated the dynamics of Arg and Lys side chains of the fruit fly Antennapedia homeodomain in the free state and in the complex with target DNA. We measured 15N relaxation for Arg and Lys side chains at two magnetic fields, from which generalized order parameters for the cationic groups were determined. Mobility of the R5 side chain, which makes hydrogen bonds with a thymine base in the DNA minor groove, was greatly dampened. Several Lys and Arg side chains that form intermolecular ion pairs with DNA phosphates were found to retain high mobility with the order parameter being <0.6 in the DNA-bound state. Interestingly, some of the interfacial cationic groups in the complex were more mobile than in the free protein. The retained or enhanced mobility of the Arg and Lys side chains in the complex should mitigate the overall loss of conformational entropy in the protein-DNA association and allow dynamic molecular recognition.

Selective Distance Measurements Using Triple Spin Labeling with Gd3+, Mn2+, and a Nitroxide

Distance measurements by pulse electron paramagnetic resonance techniques, such as double electron-electron resonance (DEER, also called PELDOR), have become an established tool to explore structural properties of biomacromolecules and their assemblies. In such measurements a pair of spin labels provides a single distance constraint. Here we show that by employing three different types of spin labels that differ in their spectroscopic and spin dynamics properties it is possible to extract three independent distances from a single sample. We demonstrate this using the Antennapedia homeodomain orthogonally labeled with Gd³⁺ and Mn²⁺ tags in complex with its cognate DNA binding site labeled with a nitroxide.

Landscape of Pleiotropic Proteins Causing Human Disease: Structural and System Biology Insights

Pleiotropy is the phenomenon by which the same gene can result in multiple phenotypes. Pleiotropic proteins are emerging as important contributors to rare and common disorders. Nevertheless, little is known on the mechanisms underlying pleiotropy and the characteristic of pleiotropic proteins. We analyzed disease-causing proteins reported in UniProt and observed that 12% are pleiotropic (variants in the same protein cause more than one disease). Pleiotropic proteins were enriched in deleterious and rare variants, but not in common variants. Pleiotropic proteins were more likely to be involved in the pathogenesis of neoplasms, neurological, and circulatory diseases and congenital malformations, whereas non-pleiotropic proteins in endocrine and metabolic disorders. Pleiotropic proteins were more essential and had a higher number of interacting partners compared with non-pleiotropic proteins. Significantly more pleiotropic than non-pleiotropic proteins contained at least one intrinsically long disordered region (P < 0.001). Deleterious variants occurring in structurally disordered regions were more commonly found in pleiotropic, rather than non-pleiotropic proteins. In conclusion, pleiotropic proteins are an important contributor to human disease. They represent a biologically different class of proteins compared with non-pleiotropic proteins and a better understanding of their characteristics and genetic variants can greatly aid in the interpretation of genetic studies and drug design.

Elastic Scattering Spectroscopy (ESS): an Instrument-Concept for Dynamics of Complex (Bio-) Systems From Elastic Neutron Scattering

A new type of neutron-scattering spectroscopy is presented that is designed specifically to measure dynamics in bio-systems that are difficult to obtain in any other way. The temporal information is largely model-free and is analogous to relaxation processes measured with dielectric spectroscopy, but provides additional spacial and geometric aspects of the underlying dynamics. Numerical simulations of the basic instrument design show the neutron beam can be highly focussed, giving efficiency gains that enable the use of small samples. Although we concentrate on continuous neutron sources, the extension to pulsed neutron sources is proposed, both requiring minimal data-treatment and being broadly analogous with dielectric spectroscopy, they will open the study of dynamics to new areas of biophysics.

Dynamic Local Polymorphisms in the Gbx1 Homeodomain Induced by DNA Binding

The Gastrulation Brain Homeobox 1 (Gbx1) gene encodes the Gbx1 homeodomain that targets TAATTA motifs in double-stranded DNA (dsDNA). Residues Glu17 and Arg52 in Gbx1 form a salt bridge, which is preserved in crystal structures and molecular dynamics simulations of homologous homeodomain-DNA complexes. In contrast, our nuclear magnetic resonance (NMR) studies show that DNA binding to Gbx1 induces dynamic local polymorphisms, which include breaking of the Glu17-Arg52 salt bridge. To study this interaction, we produced a variant with Glu17Arg and Arg52Glu mutations, which exhibited the same fold as the wild-type protein, but a 2-fold reduction in affinity for dsDNA. Analysis of the NMR structures of the Gbx1 homeodomain in the free form, the Gbx1[E17R,R52E] variant, and a Gbx1 homeodomain-DNA complex showed that stabilizing interactions of the Arg52 side chain with the DNA backbone are facilitated by transient breakage of the Glu17-Arg52 salt bridge in the DNA-bound Gbx1.

The genetic regulatory network centered on Pto-Wuschela and its targets involved in wood formation revealed by association studies

Transcription factors (TFs) regulate gene expression and can strongly affect phenotypes. However, few studies have examined TF variants and TF interactions with their targets in plants. Here, we used genetic association in 435 unrelated individuals of Populus tomentosa to explore the variants in Pto-Wuschela and its targets to decipher the genetic regulatory network of Pto-Wuschela. Our bioinformatics and co-expression analysis identified 53 genes with the motif TCACGTGA as putative targets of Pto-Wuschela. Single-marker association analysis showed that Pto-Wuschela was associated with wood properties, which is in agreement with the observation that it has higher expression in stem vascular tissues in Populus. Also, SNPs in the 53 targets were associated with growth or wood properties under additive or dominance effects, suggesting these genes and Pto-Wuschela may act in the same genetic pathways that affect variation in these quantitative traits. Epistasis analysis indicated that 75.5% of these genes directly or indirectly interacted Pto-Wuschela, revealing the coordinated genetic regulatory network formed by Pto-Wuschela and its targets. Thus, our study provides an alternative method for dissection of the interactions between a TF and its targets, which will strength our understanding of the regulatory roles of TFs in complex traits in plants.

Fuzzy complexes: Specific binding without complete folding

Specific molecular recognition is assumed to require a well-defined set of contacts and devoid of conformational and interaction ambiguities. Growing experimental evidence demonstrates however, that structural multiplicity or dynamic disorder can be retained in protein complexes, termed as fuzziness. Fuzzy regions establish alternative contacts between specific partners usually via transient interactions. Nature often tailors the dynamic properties of these segments via post-translational modifications or alternative splicing to fine-tune affinity. Most experimentally characterized fuzzy complexes are involved in regulation of gene-expression, signal transduction and cell-cycle regulation. Fuzziness is also characteristic to viral protein complexes, cytoskeleton structure, and surprisingly in a few metabolic enzymes. A plausible role of fuzzy complexes in increasing half-life of intrinsically disordered proteins is also discussed.

Interaction investigations of HipA binding to HipB dimer and HipB dimer + DNA complex: a molecular dynamics simulation study

We carried out molecular dynamics simulations and free energy calculations for a series of ternary and diplex models for the HipA protein, HipB dimer, and DNA molecule to address the mechanism of HipA sequestration and the binding order of events from apo HipB/HipA to 2HipA + HipB dimer + DNA complex. The results revealed that the combination of DNA with the HipB dimer is energetically favorable for the combination of HipB dimer with HipA protein. The binding of DNA to HipB dimer induces a long-range allosteric communication from the HipB2 -DNA interface to the HipA-HipB2 interface, which involves the closeness of α1 helices of HipB dimer to HipA protein and formations of extra hydrogen bonds in the HipA-HipB2 interface through the extension of α2/3 helices in the HipB dimer. These simulated results suggested that the DNA molecule, as a regulative media, modulates the HipB dimer conformation, consequently increasing the interactions of HipB dimer with the HipA proteins, which explains the mechanism of HipA sequestration reported by the previous experiment. Simultaneously, these simulations also explored that the thermodynamic binding order in a simulated physiological environment, that is, the HipB dimer first bind to DNA to form HipB dimer + DNA complex, then capturing strongly the HipA proteins to form a ternary complex, 2HipA + HipB dimer + DNA, for sequestrating HipA in the nucleoid.

Hox gene dysregulation in acute myeloid leukemia

ABSTRACT: 

In humans, class I homeobox genes (HOX genes) are distributed in four clusters. Upstream regulators include transcriptional activators and members of the CDX family of transcription factors. HOX genes encode proteins and need cofactor interactions, to increase their specificity and selectivity. HOX genes contribute to the organization and regulation of hematopoiesis by controlling the balance between proliferation and differentiation. Changes in HOX gene expression can be associated with chromosomal rearrangements generating fusion genes, such as those involving MLL and NUP98, or molecular defects, such as mutations in NPM1 and CEBPA for example. Several miRNAs are involved in the control of HOX gene expression and their expression correlates with HOX gene dysregulation. HOX genes dysregulation is a dominant mechanism of leukemic transformation. A better knowledge of their target genes and the mechanisms by which their dysregulated expression contributes to leukemogenesis could lead to the development of new drugs.

Binding polymorphism in the DNA bound state of the Pdx1 homeodomain

The subtle effects of DNA-protein recognition are illustrated in the homeodomain fold. This is one of several small DNA binding motifs that, in spite of limited DNA binding specificity, adopts crucial, specific roles when incorporated in a transcription factor. The homeodomain is composed of a 3-helix domain and a mobile N-terminal arm. Helix 3 (the recognition helix) interacts with the DNA bases through the major groove, while the N-terminal arm becomes ordered upon binding a specific sequence through the minor groove. Although many structural studies have characterized the DNA binding properties of homeodomains, the factors behind the binding specificity are still difficult to elucidate. A crystal structure of the Pdx1 homeodomain bound to DNA (PDB 2H1K) obtained previously in our lab shows two complexes with differences in the conformation of the N-terminal arm, major groove contacts, and backbone contacts, raising new questions about the DNA recognition process by homeodomains. Here, we carry out fully atomistic Molecular Dynamics simulations both in crystal and aqueous environments in order to elucidate the nature of the difference in binding contacts. The crystal simulations reproduce the X-ray experimental structures well. In the absence of crystal packing constraints, the differences between the two complexes increase during the solution simulations. Thus, the conformational differences are not an artifact of crystal packing. In solution, the homeodomain with a disordered N-terminal arm repositions to a partially specific orientation. Both the crystal and aqueous simulations support the existence of different stable binding conformers identified in the original crystallographic data with different degrees of specificity. We propose that protein-protein and protein-DNA interactions favor a subset of the possible conformations. This flexibility in DNA binding may facilitate multiple functions for the same transcription factor.

Thermodynamic and structural determinants of differential Pdx1 binding to elements from the insulin and IAPP promoters

In adult mammals, the production of insulin and other peptide hormones, such as the islet amyloid polypeptide (IAPP), is limited to β-cells due to tissue-specific expression of a set of transcription factors, the best known of which is pancreatic duodenal homeobox protein 1 (Pdx1). Like many homeodomain transcription factors, Pdx1 binds to a core DNA recognition sequence containing the tetranucleotide 5'-TAAT-3'; its consensus recognition element is 5'-CTCTAAT(T/G)AG-3'. Currently, a complete thermodynamic profile of Pdx1 binding to near-consensus and native promoter sequences has not been established, obscuring the mechanism of target site selection by this critical transcription factor. Strikingly, while Pdx1 responsive elements in the human insulin promoter conform to the pentanucleotide 5'-CTAAT-3' sequence, the Pdx1 responsive elements in the human iapp promoter all contain a substitution to 5'-TTAAT-3'. The crystal structure of Pdx1 bound to the consensus nucleotide sequence does not explain how Pdx1 identifies this natural variation, if it does at all. Here we report a combination of isothermal calorimetric titrations, NMR spectroscopy, and extensive multi-microsecond molecular dynamics calculations of Pdx1 that define its interactions with a panel of natural promoter elements and consensus-derived sequences. Our results show a small preference of Pdx1 for a C base 5' relative to the core TAAT promoter element. Molecular mechanics calculations, corroborated by experimental NMR data, lead to a rational explanation for sequence discrimination at this position. Taken together, our results suggest a molecular mechanism for differential Pdx1 affinity to elements from the insulin and iapp promoter sequences.

Crystal structure of the human NKX2.5 homeodomain in complex with DNA target

NKX2.5 is a homeodomain containing transcription factor regulating cardiac formation and function, and its mutations are linked to congenital heart disease. Here we provide the first report of the crystal structure of the NKX2.5 homeodomain in complex with double-stranded DNA of its endogenous target, locating within the proximal promoter -242 site of the atrial natriuretic factor gene. The crystal structure, determined at 1.8 Å resolution, demonstrates that NKX2.5 homeodomains occupy both DNA binding sites separated by five nucleotides without physical interaction between themselves. The two homeodomains show identical conformation despite the differences in the DNA sequences they bind, and no significant bending of the DNA was observed. Tyr54, absolutely conserved in NK2 family proteins, mediates sequence-specific interaction with the TAAG motif. This high resolution crystal structure of NKX2.5 protein provides a detailed picture of protein and DNA interactions, which allows us to predict DNA binding of mutants identified in human patients.

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.

Structural basis for homeodomain recognition by the cell-cycle regulator Geminin

Homeodomain-containing transcription factors play a fundamental role in the regulation of numerous developmental and cellular processes. Their multiple regulatory functions are accomplished through context-dependent inputs of target DNA sequences and collaborating protein partners. Previous studies have well established the sequence-specific DNA binding to homeodomains; however, little is known about how protein partners regulate their functions through targeting homeodomains. Here we report the solution structure of the Hox homeodomain in complex with the cell-cycle regulator, Geminin, which inhibits Hox transcriptional activity and enrolls Hox in cell proliferative control. Side-chain carboxylates of glutamates and aspartates in the C terminus of Geminin generate an overall charge pattern resembling the DNA phosphate backbone. These residues provide electrostatic interactions with homeodomain, which combine with the van der Waals contacts to form the stereospecific complex. We further showed that the interaction with Geminin is homeodomain subclass-selective and Hox paralog-specific, which relies on the stapling role of residues R43 and M54 in helix III and the basic amino acid cluster in the N terminus. Interestingly, we found that the C-terminal residue Ser184 of Geminin could be phosphorylated by Casein kinase II, resulting in the enhanced binding to Hox and more potent inhibitory effect on Hox transcriptional activity, indicating an additional layer of regulation. This structure provides insight into the molecular mechanism underlying homeodomain-protein recognition and may serve as a paradigm for interactions between homeodomains and DNA-competitive peptide inhibitors.

The use of information theory in evolutionary biology

Information is a key concept in evolutionary biology. Information stored in a biological organism's genome is used to generate the organism and to maintain and control it. Information is also that which evolves. When a population adapts to a local environment, information about this environment is fixed in a representative genome. However, when an environment changes, information can be lost. At the same time, information is processed by animal brains to survive in complex environments, and the capacity for information processing also evolves. Here, I review applications of information theory to the evolution of proteins and to the evolution of information processing in simulated agents that adapt to perform a complex task.

Structural and biophysical insights into the ligand-free Pitx2 homeodomain and a ring dermoid of the cornea inducing homeodomain mutant

The homeodomain-containing transcription factor Pitx2 (pituitary homeobox protein 2) is present in many developing embryonic tissues, including the heart. Its homeodomain is responsible for the recognition and binding to target DNA sequences and thus constitutes a major functional unit in the Pitx2 protein. Nuclear magnetic resonance techniques were employed to determine the solution structure of the native Pitx2 homeodomain and a R24H mutant that causes autosomal dominantly inherited ring dermoid of the cornea syndrome. The structures reveal that both isoforms possess the canonical homeodomain fold. However, the R24H mutation results in a 2-fold increase in DNA binding affinity and a 5 °C decrease in thermal stability, while changing the dynamic environment of the homeodomain only locally. When introduced into full-length Pitx2c, the mutation results in an only 25% loss of transactivation activity. Our data correlate well with clinical observations suggesting a milder deficiency for the R24H mutation compared to those of other Pitx2 homeodomain mutations.

Preparation and optimization of protein-DNA complexes suitable for detailed NMR studies

This chapter describes the methods to form and optimize samples of protein-DNA complexes that are suitable for detailed structure and dynamics studies by NMR spectroscopy.

Exploring translocation of proteins on DNA by NMR

While an extensive body of knowledge has accumulated on the structures of transcription factors, DNA and their complexes from both NMR and crystallography, much less is known at a molecular level regarding the mechanisms whereby transcription factors locate their specific DNA target site within an overwhelming sea of non-specific DNA sites. Indirect kinetic data suggested that three processes are involved in the search procedure: jumping by dissociation of the protein from the DNA followed by re-association at another site, direct transfer from one DNA molecule or segment to another, and one-dimensional sliding. In this brief perspective I summarize recent NMR developments from our laboratory that have permitted direct characterization of the species and molecular mechanisms involved in the target search process, including the detection of highly transient sparsely-populated states. The main tool in these studies involves the application of paramagnetic relaxation enhancement, supplemented by z-exchange spectroscopy, lineshape analysis and residual dipolar couplings. These studies led to the first direct demonstration of rotation-coupled sliding of a protein along the DNA and the direct transfer of a protein from one DNA molecule to another without dissociating into free solution.

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.

Dynamic protein-DNA recognition: beyond what can be seen

Traditionally, specific DNA recognition is thought to rely on static contacts with the bases or phosphates. Recent results, however, indicate that residues far outside the binding context can crucially influence selectivity or binding affinity via transient, dynamic interactions with the DNA binding interface. These regions usually do not adopt a well-defined structure, even when bound to DNA, and thus form a fuzzy complex. Here, we propose the existence of a dynamic DNA readout mechanism, wherein distant segments modulate conformational preferences, flexibility or spacing of the DNA binding motifs or serve as competitive partners. Despite their low sequence similarity, these intrinsically disordered regions are often conserved at the structural level, and exploited for regulation of the transcription machinery via protein-protein interactions, post-translational modifications or alternative splicing.

Nuclear magnetic resonance analysis of protein-DNA interactions

Recent methodological and instrumental advances in solution-state nuclear magnetic resonance have opened up the way to investigating challenging problems in structural biology such as large macromolecular complexes. This review focuses on the experimental strategies currently employed to solve structures of protein-DNA complexes and to analyse their dynamics. It highlights how these approaches can help in understanding detailed molecular mechanisms of target recognition.

Homeodomain subtypes and functional diversity

The homeodomain is a protein domain of about 60 amino acids that is encoded by homeobox genes. The homeodomain is a DNA binding domain, and hence homeodomain proteins are essentially transcription factors (TFs). They have been shown to play major roles in many developmental processes of animals, as well as fungi and plants. A primary function of homeodomain proteins is to regulate the expression of other genes in development and differentiation. Thousands of homeobox genes have been identified, and they can be grouped into many different classes. Often other conserved protein domains are found linked to a homeodomain. Several particular types of homeobox genes are organized into chromosomal clusters. The best-known cluster, the HOX cluster, is found in all bilaterian animals. Tetrapods contain four HOX clusters that arose through duplication in early vertebrate evolution. The genes in these clusters are called Hox genes. Lower chordates, insects and nematodes tend to have only one HOX cluster. Of particular interest is that many of the HOX cluster genes function in the process of pattern formation along the anterior-posterior body axis. Many other types of homeodomain proteins play roles in the determination of cell fates and cell differentiation. Homeobox genes thus perform key roles for all aspects of the development of an organism.

DNA search efficiency is modulated by charge composition and distribution in the intrinsically disordered tail

Intrinsically disordered tails are common in DNA-binding proteins and can affect their search efficiency on nonspecific DNA by promoting the brachiation dynamics of intersegment transfer. During brachiation, the protein jumps between distant DNA regions via an intermediate state in which the tail and globular moieties are bound to different DNA segments. While the disordered tail must be long and positively charged to facilitate DNA search, the effect of its residue sequence on brachiation is unknown. We explored this issue using the NK-2 and Antp homeodomain transcription factors. We designed 566 NK-2 tail-variants and 55 Antp tail-variants having different net charges and positive charge distributions and studied their dynamics and DNA search efficiencies using coarse-grained molecular dynamics simulations. More intersegment transfers occur when the tail is moderately positively charged and the positive charges are clustered together in the middle of the tail or towards its N terminus. The presence of a negatively charged residue does not significantly affect protein brachiation, although it is likely that the presence of many negatively charged residues will complicate the DNA search mechanism. A bioinformatic analysis of 1,384 wild-type homeodomains illustrates that the charge composition and distribution in their N-tail sequences are consistent with an optimal charge pattern to promote intersegment transfer. Our study thus indicates that the residue sequence of the disordered tails of DNA-binding proteins has unique characteristics that were evolutionarily selected to achieve optimized function and suggests that the sequence-structure-function paradigm known for structured proteins is valid for intrinsically disordered proteins as well.

Interpreting protein/DNA interactions: distinguishing specific from non-specific and electrostatic from non-electrostatic components

We discuss the effectiveness of existing methods for understanding the forces driving the formation of specific protein-DNA complexes. Theoretical approaches using the Poisson-Boltzmann (PB) equation to analyse interactions between these highly charged macromolecules to form known structures are contrasted with an empirical approach that analyses the effects of salt on the stability of these complexes and assumes that release of counter-ions associated with the free DNA plays the dominant role in their formation. According to this counter-ion condensation (CC) concept, the salt-dependent part of the Gibbs energy of binding, which is defined as the electrostatic component, is fully entropic and its dependence on the salt concentration represents the number of ionic contacts present in the complex. It is shown that although this electrostatic component provides the majority of the Gibbs energy of complex formation and does not depend on the DNA sequence, the salt-independent part of the Gibbs energy--usually regarded as non-electrostatic--is sequence specific. The CC approach thus has considerable practical value for studying protein/DNA complexes, while practical applications of PB analysis have yet to demonstrate their merit.

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.

Cooperative DNA-binding and sequence-recognition mechanism of aristaless and clawless

To achieve accurate gene regulation, some homeodomain proteins bind cooperatively to DNA to increase those site specificities. We report a ternary complex structure containing two homeodomain proteins, aristaless (Al) and clawless (Cll), bound to DNA. Our results show that the extended conserved sequences of the Cll homeodomain are indispensable to cooperative DNA binding. In the Al-Cll-DNA complex structure, the residues in the extended regions are used not only for the intermolecular contacts between the two homeodomain proteins but also for the sequence-recognition mechanism of DNA by direct interactions. The residues in the extended N-terminal arm lie within the minor groove of DNA to form direct interactions with bases, whereas the extended conserved region of the C-terminus of the homeodomain interacts with Al to stabilize and localize the third alpha helix of the Cll homeodomain. This structure suggests a novel mode for the cooperativity of homeodomain proteins.