DNA-guided transcription factor cooperativity shapes face and limb mesenchyme.

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

DNA-binding specificities of human transcription factors.

Although the proteins that read the gene regulatory code, transcription factors (TFs), have been largely identified, it is not well known which sequences TFs can recognize. We have analyzed the sequence-specific binding of human TFs using high-throughput SELEX and ChIP sequencing. A total of 830 binding profiles were obtained, describing 239 distinctly different binding specificities. The models represent the majority of human TFs, approximately doubling the coverage compared to existing systematic studies. Our results reveal additional specificity determinants for a large number of factors for which a partial specificity was known, including a commonly observed A- or T-rich stretch that flanks the core motifs. Global analysis of the data revealed that homodimer orientation and spacing preferences, and base-stacking interactions, have a larger role in TF-DNA binding than previously appreciated. We further describe a binding model incorporating these features that is required to understand binding of TFs to DNA.

Genome-Wide Association Study Identifies 4 Novel Risk Loci for Small Intestinal Neuroendocrine Tumors Including a Missense Mutation in LGR5.

BACKGROUND & AIMS: Small intestinal neuroendocrine tumor (SI-NET) is a rare disease, but its incidence has increased over the past 4 decades. Understanding the genetic risk factors underlying SI-NETs can help in disease prevention and may provide clinically beneficial markers for diagnosis. Here the results of the largest genome-wide association study of SI-NETs performed to date with 405 cases and 614,666 controls are reported. METHODS: Samples from 307 patients with SI-NETs and 287,137 controls in the FinnGen study were used for the identification of SI-NET risk-associated genetic variants. The results were also meta-analyzed with summary statistics from the UK Biobank (n = 98 patients with SI-NET and n = 327,529 controls). RESULTS: We identified 6 genome-wide significant (P < 5 × 10-8) loci associated with SI-NET risk, of which 4 (near SEMA6A, LGR5, CDKAL1, and FERMT2) are novel and 2 (near LTA4H-ELK and in KIF16B) have been reported previously. Interestingly, the top hit (rs200138614; P = 1.80 × 10-19) was a missense variant (p.Cys712Phe) in the LGR5 gene, a bona-fide marker of adult intestinal stem cells and a potentiator of canonical WNT signaling. The association was validated in an independent Finnish collection of 70 patients with SI-NETs, as well as in the UK Biobank exome sequence data (n = 92 cases and n = 392,814 controls). Overexpression of LGR5 p.Cys712Phe in intestinal organoids abolished the ability of R-Spondin1 to support organoid growth, indicating that the mutation perturbed R-Spondin-LGR5 signaling. CONCLUSIONS: Our study is the largest genome-wide association study to date on SI-NETs and reported 4 new associated genome-wide association study loci, including a novel missense mutation (rs200138614, p.Cys712Phe) in LGR5, a canonical marker of adult intestinal stem cells.

The interaction landscape between transcription factors and the nucleosome.

Nucleosomes cover most of the genome and are thought to be displaced by transcription factors in regions that direct gene expression. However, the modes of interaction between transcription factors and nucleosomal DNA remain largely unknown. Here we systematically explore interactions between the nucleosome and 220 transcription factors representing diverse structural families. Consistent with earlier observations, we find that the majority of the studied transcription factors have less access to nucleosomal DNA than to free DNA. The motifs recovered from transcription factors bound to nucleosomal and free DNA are generally similar. However, steric hindrance and scaffolding by the nucleosome result in specific positioning and orientation of the motifs. Many transcription factors preferentially bind close to the end of nucleosomal DNA, or to periodic positions on the solvent-exposed side of the DNA. In addition, several transcription factors usually bind to nucleosomal DNA in a particular orientation. Some transcription factors specifically interact with DNA located at the dyad position at which only one DNA gyre is wound, whereas other transcription factors prefer sites spanning two DNA gyres and bind specifically to each of them. Our work reveals notable differences in the binding of transcription factors to free and nucleosomal DNA, and uncovers a diverse interaction landscape between transcription factors and the nucleosome.

Binding specificities of human RNA-binding proteins toward structured and linear RNA sequences.

RNA-binding proteins (RBPs) regulate RNA metabolism at multiple levels by affecting splicing of nascent transcripts, RNA folding, base modification, transport, localization, translation, and stability. Despite their central role in RNA function, the RNA-binding specificities of most RBPs remain unknown or incompletely defined. To address this, we have assembled a genome-scale collection of RBPs and their RNA-binding domains (RBDs) and assessed their specificities using high-throughput RNA-SELEX (HTR-SELEX). Approximately 70% of RBPs for which we obtained a motif bound to short linear sequences, whereas ∼30% preferred structured motifs folding into stem-loops. We also found that many RBPs can bind to multiple distinctly different motifs. Analysis of the matches of the motifs in human genomic sequences suggested novel roles for many RBPs. We found that three cytoplasmic proteins-ZC3H12A, ZC3H12B, and ZC3H12C-bound to motifs resembling the splice donor sequence, suggesting that these proteins are involved in degradation of cytoplasmic viral and/or unspliced transcripts. Structural analysis revealed that the RNA motif was not bound by the conventional C3H1 RNA-binding domain of ZC3H12B. Instead, the RNA motif was bound by the ZC3H12B's PilT N terminus (PIN) RNase domain, revealing a potential mechanism by which unconventional RBDs containing active sites or molecule-binding pockets could interact with short, structured RNA molecules. Our collection containing 145 high-resolution binding specificity models for 86 RBPs is the largest systematic resource for the analysis of human RBPs and will greatly facilitate future analysis of the various biological roles of this important class of proteins.

DNA-dependent formation of transcription factor pairs alters their binding specificity.

Gene expression is regulated by transcription factors (TFs), proteins that recognize short DNA sequence motifs. Such sequences are very common in the human genome, and an important determinant of the specificity of gene expression is the cooperative binding of multiple TFs to closely located motifs. However, interactions between DNA-bound TFs have not been systematically characterized. To identify TF pairs that bind cooperatively to DNA, and to characterize their spacing and orientation preferences, we have performed consecutive affinity-purification systematic evolution of ligands by exponential enrichment (CAP-SELEX) analysis of 9,400 TF-TF-DNA interactions. This analysis revealed 315 TF-TF interactions recognizing 618 heterodimeric motifs, most of which have not been previously described. The observed cooperativity occurred promiscuously between TFs from diverse structural families. Structural analysis of the TF pairs, including a novel crystal structure of MEIS1 and DLX3 bound to their identified recognition site, revealed that the interactions between the TFs were predominantly mediated by DNA. Most TF pair sites identified involved a large overlap between individual TF recognition motifs, and resulted in recognition of composite sites that were markedly different from the individual TF's motifs. Together, our results indicate that the DNA molecule commonly plays an active role in cooperative interactions that define the gene regulatory lexicon.

Damaging heterozygous mutations in NFKB1 lead to diverse immunologic phenotypes.

BACKGROUND: The nuclear factor κ light-chain enhancer of activated B cells (NF-κB) signaling pathway is a key regulator of immune responses. Accordingly, mutations in several NF-κB pathway genes cause immunodeficiency. OBJECTIVE: We sought to identify the cause of disease in 3 unrelated Finnish kindreds with variable symptoms of immunodeficiency and autoinflammation. METHODS: We applied genetic linkage analysis and next-generation sequencing and functional analyses of NFKB1 and its mutated alleles. RESULTS: In all affected subjects we detected novel heterozygous variants in NFKB1, encoding for p50/p105. Symptoms in variant carriers differed depending on the mutation. Patients harboring a p.I553M variant presented with antibody deficiency, infection susceptibility, and multiorgan autoimmunity. Patients with a p.H67R substitution had antibody deficiency and experienced autoinflammatory episodes, including aphthae, gastrointestinal disease, febrile attacks, and small-vessel vasculitis characteristic of Behçet disease. Patients with a p.R157X stop-gain experienced hyperinflammatory responses to surgery and showed enhanced inflammasome activation. In functional analyses the p.R157X variant caused proteasome-dependent degradation of both the truncated and wild-type proteins, leading to a dramatic loss of p50/p105. The p.H67R variant reduced nuclear entry of p50 and showed decreased transcriptional activity in luciferase reporter assays. The p.I553M mutation in turn showed no change in p50 function but exhibited reduced p105 phosphorylation and stability. Affinity purification mass spectrometry also demonstrated that both missense variants led to altered protein-protein interactions. CONCLUSION: Our findings broaden the scope of phenotypes caused by mutations in NFKB1 and suggest that a subset of autoinflammatory diseases, such as Behçet disease, can be caused by rare monogenic variants in genes of the NF-κB pathway.

Whole-Genome Sequencing Identifies STAT4 as a Putative Susceptibility Gene in Classic Kaposi Sarcoma.

BACKGROUND: Classic Kaposi sarcoma (cKS) is an inflammatory tumor caused by human herpesvirus 8 (HHV-8) commonly observed in elderly men of Mediterranean origin. We studied a Finnish family of 5 affected individuals in 2 generations. Except for atypical mycobacterial infection of the index case, the affected individuals did not have notable histories of infection. METHODS: We performed genome and exome sequencing and mapped shared chromosomal regions to identify genetic predisposition in the family. RESULTS: We identified 12 protein-coding candidate variants that segregated in the 3 affected cousins from whom we had samples. The affected mother of the index case was an obligatory carrier. Among the 12 candidates was a rare heterozygous substitution rs141331848 (c.1337C>T, p.Thr446Ile) in the DNA-binding domain of STAT4. The variant was not present in 242 Finnish control genomes or 180 additional regional controls. Activated T-helper cells from the HHV-8-negative variant carriers showed reduced interferon γ production, compared with age and sex matched wild-type individuals. We screened STAT4 in additional 18 familial KS cases and the variant site from 56 sporadic KS cases but detected no pathogenic mutations. CONCLUSIONS: Our data suggest that STAT4 is a potential cKS-predisposition gene, but further functional and genetic validation is needed.

Loss of SUFU function in familial multiple meningioma.

Meningiomas are the most common primary tumors of the CNS and account for up to 30% of all CNS tumors. An increased risk of meningiomas has been associated with certain tumor-susceptibility syndromes, especially neurofibromatosis type II, but no gene defects predisposing to isolated familial meningiomas have thus far been identified. Here, we report on a family of five meningioma-affected siblings, four of whom have multiple tumors. No NF2 mutations were identified in the germline or tumors. We combined genome-wide linkage analysis and exome sequencing, and we identified in suppressor of fused homolog (Drosophila), SUFU, a c.367C>T (p.Arg123Cys) mutation segregating with the meningiomas in the family. The variation was not present in healthy controls, and all seven meningiomas analyzed displayed loss of the wild-type allele according to the classic two-hit model for tumor-suppressor genes. In silico modeling predicted the variant to affect the tertiary structure of the protein, and functional analyses showed that the activity of the altered SUFU was significantly reduced and therefore led to dysregulated hedgehog (Hh) signaling. SUFU is a known tumor-suppressor gene previously associated with childhood medulloblastoma predisposition. Our genetic and functional analyses indicate that germline mutations in SUFU also predispose to meningiomas, particularly to multiple meningiomas. It is possible that other genic mutations resulting in aberrant activation of the Hh pathway might underlie meningioma predisposition in families with an unknown etiology.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme.

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how 'Coordinator', a long DNA motif comprised of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, while HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in shared regulation of genes involved in cell-type and positional identities, and ultimately shapes facial morphology and evolution.

Structural insights into the interaction between transcription factors and the nucleosome.

In eukaryotic cells, DNA interacts with two main types of binding proteins: transcription factors and histones. Histones form the core of nucleosomes and display weak sequence preference owing to differences in bendability of different DNA sequences. By contrast, the affinity of transcription factors is highly dependent on DNA sequence - all sequences are bound with moderate affinity, but only few specific sequences are bound more tightly via molecular recognition of the DNA bases. Transcription factors can interact with nucleosomes directly by recognizing nucleosome-associated DNA and also indirectly by recruiting histone-modifying enzymes and nucleosome remodelers. These interactions result in sequence-dependent formation of a pattern of open and closed chromatin, where specific positions are occupied by transcription factors, histone-modifying enzymes, and modified histones. These patterns are then recognized by large DNA-associated macromolecular complexes such as cohesin and RNA polymerase II, which are involved in regulation of higher-order chromatin structure and transcription, respectively. Here, we review recent work that has provided structural and mechanistic insight into the interactions between all these classes of DNA-associated proteins.

Structural study and thermodynamic characterization of inhibitor binding to lumazine synthase from Bacillus anthracis.

The crystal structure of lumazine synthase from Bacillus anthracis was solved by molecular replacement and refined to R(cryst) = 23.7% (R(free) = 28.4%) at a resolution of 3.5 A. The structure reveals the icosahedral symmetry of the enzyme and specific features of the active site that are unique in comparison with previously determined orthologues. The application of isothermal titration calorimetry in combination with enzyme kinetics showed that three designed pyrimidine derivatives bind to lumazine synthase with micromolar dissociation constants and competitively inhibit the catalytic reaction. Structure-based modelling suggested the binding modes of the inhibitors in the active site and allowed an estimation of the possible contacts formed upon binding. The results provide a structural framework for the design of antibiotics active against B. anthracis.

Virtual screening, selection and development of a benzindolone structural scaffold for inhibition of lumazine synthase.

Virtual screening of a library of commercially available compounds versus the structure of Mycobacterium tuberculosis lumazine synthase identified 2-(2-oxo-1,2-dihydrobenzo[cd]indole-6-sulfonamido)acetic acid (9) as a possible lead compound. Compound 9 proved to be an effective inhibitor of M. tuberculosis lumazine synthase with a K(i) of 70microM. Lead optimization through replacement of the carboxymethylsulfonamide sidechain with sulfonamides substituted with alkyl phosphates led to a four-carbon phosphate 38 that displayed a moderate increase in enzyme inhibitory activity (K(i) 38microM). Molecular modeling based on known lumazine synthase/inhibitor crystal structures suggests that the main forces stabilizing the present benzindolone/enzyme complexes involve pi-pi stacking interactions with Trp27 and hydrogen bonding of the phosphates with Arg128, the backbone nitrogens of Gly85 and Gln86, and the side chain hydroxyl of Thr87.

Second career of a biosynthetic enzyme: Lumazine synthase as a virus-like nanoparticle in vaccine development.

Naturally occurring and computationally ab initio designed protein cages can now be considered as extremely suitable materials for new developments in nanotechnology. Via self-assembly from single identical or non-identical protomers large oligomeric particles can be formed. Virus-like particles have today found a number of quite successful applications in the development of new vaccines. Complex chimeric nanoparticles can serve as suitable platforms for the presentation of natural or designed antigens to the immune system of the host. The scaffolds can be cage forming highly symmetric biological macromolecules like lumazine synthase or symmetric self-assembling virus-like particles generated by computational ab initio design. Symmetric nanoparticle carriers display a structurally ordered array of immunogens. This feature can lead to a more favorable interaction with B-cell receptors, in comparison to the administration of single recombinant immunogens. Several pre-clinical animal studies and clinical studies have recently pointed out the efficiency of nanoparticle antigens produced recombinantly in creating strong immune responses against infectious diseases like HIV, Malaria, Borrelia, Influenza.

Structural insights into the adaptation of proliferating cell nuclear antigen (PCNA) from Haloferax volcanii to a high-salt environment.

The sliding clamp proliferating cell nuclear antigen (PCNA) plays vital roles in many aspects of DNA replication and repair in eukaryotic cells and in archaea. Realising the full potential of archaea as a model for PCNA function requires a combination of biochemical and genetic approaches. In order to provide a platform for subsequent reverse genetic analysis, PCNA from the halophilic archaeon Haloferax volcanii was subjected to crystallographic analysis. The gene was cloned and expressed in Escherichia coli and the protein was purified by affinity chromatography and crystallized by the vapour-diffusion technique. The structure was determined by molecular replacement and refined at 3.5 A resolution to a final R factor of 23.7% (R(free) = 25%). PCNA from H. volcanii was found to be homotrimeric and to resemble other homotrimeric PCNA clamps but with several differences that appear to be associated with adaptation of the protein to the high intracellular salt concentrations found in H. volcanii cells.

Novel TMEM173 Mutation and the Role of Disease Modifying Alleles.

Upon binding to pathogen or self-derived cytosolic nucleic acids cyclic GMP-AMP synthase (cGAS) triggers the production of cGAMP that further activates transmembrane protein STING. Upon activation STING translocates from ER via Golgi to vesicles. Monogenic STING gain-of-function mutations cause early-onset type I interferonopathy, with disease presentation ranging from fatal vasculopathy to mild chilblain lupus. Molecular mechanisms underlying the variable phenotype-genotype correlation are presently unclear. Here, we report a novel gain-of-function G207E STING mutation causing a distinct phenotype with alopecia, photosensitivity, thyroid dysfunction, and features of STING-associated vasculopathy with onset in infancy (SAVI), such as livedo reticularis, skin vasculitis, nasal septum perforation, facial erythema, and bacterial infections. Polymorphism in TMEM173 and IFIH1 showed variable penetrance in the affected family, implying contribution to varying phenotype spectrum. The G207E mutation constitutively activates inflammation-related pathways in vitro, and causes aberrant interferon signature and inflammasome activation in patient PBMCs. Treatment with Janus kinase 1 and 2 (JAK1/2) inhibitor baricitinib was beneficiary for a vasculitic ulcer, induced hair regrowth and improved overall well-being in one patient. Protein-protein interactions propose impaired cellular trafficking of G207E mutant. These findings reveal the molecular landscape of STING and propose common polymorphisms in TMEM173 and IFIH1 as likely modifiers of the phenotype.

Impact of constitutional TET2 haploinsufficiency on molecular and clinical phenotype in humans.

Clonal hematopoiesis driven by somatic heterozygous TET2 loss is linked to malignant degeneration via consequent aberrant DNA methylation, and possibly to cardiovascular disease via increased cytokine and chemokine expression as reported in mice. Here, we discover a germline TET2 mutation in a lymphoma family. We observe neither unusual predisposition to atherosclerosis nor abnormal pro-inflammatory cytokine or chemokine expression. The latter finding is confirmed in cells from three additional unrelated TET2 germline mutation carriers. The TET2 defect elevates blood DNA methylation levels, especially at active enhancers and cell-type specific regulatory regions with binding sequences of master transcription factors involved in hematopoiesis. The regions display reduced methylation relative to all open chromatin regions in four DNMT3A germline mutation carriers, potentially due to TET2-mediated oxidation. Our findings provide insight into the interplay between epigenetic modulators and transcription factor activity in hematological neoplasia, but do not confirm the putative role of TET2 in atherosclerosis.

Two distinct DNA sequences recognized by transcription factors represent enthalpy and entropy optima.

Most transcription factors (TFs) can bind to a population of sequences closely related to a single optimal site. However, some TFs can bind to two distinct sequences that represent two local optima in the Gibbs free energy of binding (ΔG). To determine the molecular mechanism behind this effect, we solved the structures of human HOXB13 and CDX2 bound to their two optimal DNA sequences, CAATAAA and TCGTAAA. Thermodynamic analyses by isothermal titration calorimetry revealed that both sites were bound with similar ΔG. However, the interaction with the CAA sequence was driven by change in enthalpy (ΔH), whereas the TCG site was bound with similar affinity due to smaller loss of entropy (ΔS). This thermodynamic mechanism that leads to at least two local optima likely affects many macromolecular interactions, as ΔG depends on two partially independent variables ΔH and ΔS according to the central equation of thermodynamics, ΔG = ΔH - TΔS.

Isolation, crystallization and preliminary crystallographic analysis of Salmonella typhimurium uridine phosphorylase crystallized with 2,2'-anhydrouridine.

Uridine phosphorylase (UPh; EC 2.4.2.3) is a member of the pyrimidine nucleoside phosphorylase family of enzymes which catalyzes the phosphorolytic cleavage of the C-N glycoside bond of uridine, with the formation of ribose 1-phosphate and uracil. This enzyme has been shown to be important in the activation and catabolism of fluoropyrimidines. Modulation of its enzymatic activity may affect the therapeutic efficacy of chemotherapeutic agents. The structural investigation of the bacterial uridine phosphorylases, both unliganded and complexed with substrate/product analogues and inhibitors, may help in understanding the catalytic mechanism of the phosphorolytic cleavage of uridine. Salmonella typhimurium uridine phosphorylase has been crystallized with 2,2'-anhydrouridine. X-ray diffraction data were collected to 2.15 A. Preliminary analysis of the diffraction data indicates that the crystal belongs to space group P2(1)2(1)2(1), with unit-cell parameters a = 88.52, b = 123.98, c = 133.52 A. The solvent content is 45.51%, assuming the presence of one hexamer molecule per asymmetric unit.

Structural perspective of cooperative transcription factor binding.

In prokaryotes, individual transcription factors (TFs) can recognize long DNA motifs that are alone sufficient to define the genes that they induce or repress. In contrast, in higher organisms that have larger genomes, TFs recognize sequences that are too short to define unique genomic positions. In addition, development of multicellular organisms requires molecular systems that are capable of executing combinatorial logical operations. Co-operative recognition of DNA by multiple TFs allows both definition of unique genomic positions in large genomes, and complex information processing at the level of individual regulatory elements. The TFs can co-operate in multiple different ways, and the precise mechanism used for co-operation determines important features of the regulatory interactions. Here, we present an overview of the structural basis of the different mechanisms by which TFs can cooperate, focusing on insight from recent functional studies and structural analyses of specific TF-TF-DNA complexes.