Crystal structure of the yeast MATalpha2/MCM1/DNA ternary complex

Experimental approaches to investigate biophysical interactions between homeodomain transcription factors and DNA

Homeodomain transcription factors (TFs) bind to specific DNA sequences to regulate the expression of target genes. Structural work has provided insight into molecular identities and aided in unraveling structural features of these TFs. However, the detailed affinity and specificity by which these TFs bind to DNA sequences is still largely unknown. Qualitative methods, such as DNA footprinting, Electrophoretic Mobility Shift Assays (EMSAs), Systematic Evolution of Ligands by Exponential Enrichment (SELEX), Bacterial One Hybrid (B1H) systems, Surface Plasmon Resonance (SPR), and Protein Binding Microarrays (PBMs) have been widely used to investigate the biochemical characteristics of TF-DNA binding events. In addition to these qualitative methods, bioinformatic approaches have also assisted in TF binding site discovery. Here we discuss the advantages and limitations of these different approaches, as well as the benefits of utilizing more quantitative approaches, such as Mechanically Induced Trapping of Molecular Interactions (MITOMI), Microscale Thermophoresis (MST) and Isothermal Titration Calorimetry (ITC), in determining the biophysical basis of binding specificity of TF-DNA complexes and improving upon existing computational approaches aimed at affinity predictions.

SEPALLATA-driven MADS transcription factor tetramerization is required for inner whorl floral organ development

MADS transcription factors are master regulators of plant reproduction and flower development. The SEPALLATA (SEP) subfamily of MADS transcription factors is required for the development of floral organs and plays roles in inflorescence architecture and development of the floral meristem. SEPALLATAs act as organizers of MADS complexes, forming both heterodimers and heterotetramers in vitro. To date, the MADS complexes characterized in angiosperm floral organ development contain at least 1 SEPALLATA protein. Whether DNA binding by SEPALLATA-containing dimeric MADS complexes is sufficient for launching floral organ identity programs, however, is not clear as only defects in floral meristem determinacy were observed in tetramerization-impaired SEPALLATA mutant proteins. Here, we used a combination of genome-wide-binding studies, high-resolution structural studies of the SEP3/AGAMOUS (AG) tetramerization domain, structure-based mutagenesis and complementation experiments in Arabidopsis (Arabidopsis thaliana) sep1 sep2 sep3 and sep1 sep2 sep3 ag-4 plants transformed with versions of SEP3 encoding tetramerization mutants. We demonstrate that while SEP3 heterodimers can bind DNA both in vitro and in vivo and recognize the majority of SEP3 wild-type-binding sites genome-wide, tetramerization is required not only for floral meristem determinacy but also for floral organ identity in the second, third, and fourth whorls.

Accurate prediction of protein function using statistics-informed graph networks

Understanding protein function is pivotal in comprehending the intricate mechanisms that underlie many crucial biological activities, with far-reaching implications in the fields of medicine, biotechnology, and drug development. However, more than 200 million proteins remain uncharacterized, and computational efforts heavily rely on protein structural information to predict annotations of varying quality. Here, we present a method that utilizes statistics-informed graph networks to predict protein functions solely from its sequence. Our method inherently characterizes evolutionary signatures, allowing for a quantitative assessment of the significance of residues that carry out specific functions. PhiGnet not only demonstrates superior performance compared to alternative approaches but also narrows the sequence-function gap, even in the absence of structural information. Our findings indicate that applying deep learning to evolutionary data can highlight functional sites at the residue level, providing valuable support for interpreting both existing properties and new functionalities of proteins in research and biomedicine.

Evolution of a new form of haploid-specific gene regulation appearing in a limited clade of ascomycete yeast species

Over evolutionary timescales, the logic and pattern of cell-type specific gene expression can remain constant, yet the molecular mechanisms underlying such regulation can drift between alternative forms. Here, we document a new example of this principle in the regulation of the haploid-specific genes in a small clade of fungal species. For most ascomycete fungal species, transcription of these genes is repressed in the a/α cell type by a heterodimer of two homeodomain proteins, Mata1 and Matα2. We show that in the species Lachancea kluyveri, most of the haploid-specific genes are regulated in this way, but repression of one haploid-specific gene (GPA1) requires, in addition to Mata1 and Matα2, a third regulatory protein, Mcm1. Model building, based on x-ray crystal structures of the three proteins, rationalizes the requirement for all three proteins: no single pair of the proteins is optimally arranged, and we show that no single pair can bring about repression. This case study exemplifies the idea that the energy of DNA binding can be "shared out" in different ways and can result in different DNA-binding solutions across different genes-while maintaining the same overall pattern of gene expression.

Cracking the Floral Quartet Code: How Do Multimers of MIKCC-Type MADS-Domain Transcription Factors Recognize Their Target Genes?

MADS-domain transcription factors (MTFs) are involved in the control of many important processes in eukaryotes. They are defined by the presence of a unique and highly conserved DNA-binding domain, the MADS domain. MTFs bind to double-stranded DNA as dimers and recognize specific sequences termed CArG boxes (such as 5'-CC(A/T)6GG-3') and similar sequences that occur hundreds of thousands of times in a typical flowering plant genome. The number of MTF-encoding genes increased by around two orders of magnitude during land plant evolution, resulting in roughly 100 genes in flowering plant genomes. This raises the question as to how dozens of different but highly similar MTFs accurately recognize the cis-regulatory elements of diverse target genes when the core binding sequence (CArG box) occurs at such a high frequency. Besides the usual processes, such as the base and shape readout of individual DNA sequences by dimers of MTFs, an important sublineage of MTFs in plants, termed MIKCC-type MTFs (MC-MTFs), has evolved an additional mechanism to increase the accurate recognition of target genes: the formation of heterotetramers of closely related proteins that bind to two CArG boxes on the same DNA strand involving DNA looping. MC-MTFs control important developmental processes in flowering plants, ranging from root and shoot to flower, fruit and seed development. The way in which MC-MTFs bind to DNA and select their target genes is hence not only of high biological interest, but also of great agronomic and economic importance. In this article, we review the interplay of the different mechanisms of target gene recognition, from the ordinary (base readout) via the extravagant (shape readout) to the idiosyncratic (recognition of the distance and orientation of two CArG boxes by heterotetramers of MC-MTFs). A special focus of our review is on the structural prerequisites of MC-MTFs that enable the specific recognition of target genes.

ZHX2 promotes HIF1α oncogenic signaling in triple-negative breast cancer

Triple-negative breast cancer (TNBC) is an aggressive and highly lethal disease, which warrants the critical need to identify new therapeutic targets. We show that Zinc Fingers and Homeoboxes 2 (ZHX2) is amplified or overexpressed in TNBC cell lines and patients. Functionally, depletion of ZHX2 inhibited TNBC cell growth and invasion in vitro, orthotopic tumor growth, and spontaneous lung metastasis in vivo. Mechanistically, ZHX2 bound with hypoxia-inducible factor (HIF) family members and positively regulated HIF1α activity in TNBC. Integrated ChIP-seq and gene expression profiling demonstrated that ZHX2 co-occupied with HIF1α on transcriptionally active promoters marked by H3K4me3 and H3K27ac, thereby promoting gene expression. Among the identified ZHX2 and HIF1α coregulated genes, overexpression of AP2B1, COX20, KDM3A, or PTGES3L could partially rescue TNBC cell growth defect by ZHX2 depletion, suggested that these downstream targets contribute to the oncogenic role of ZHX2 in an accumulative fashion. Furthermore, multiple residues (R491, R581, and R674) on ZHX2 are important in regulating its phenotype, which correspond with their roles on controlling ZHX2 transcriptional activity in TNBC cells. These studies establish that ZHX2 activates oncogenic HIF1α signaling, therefore serving as a potential therapeutic target for TNBC.

How structural biology transformed studies of transcription regulation

The past 4 decades have seen remarkable advances in our understanding of the structural basis of gene regulation. Technological advances in protein expression, nucleic acid synthesis, and structural biology made it possible to study the proteins that regulate transcription in the context of ever larger complexes containing proteins bound to DNA. This review, written on the occasion of the 50th anniversary of the founding of the Protein Data Bank focuses on the insights gained from structural studies of protein-DNA complexes and the role the PDB has played in driving this research. I cover highlights in the field, beginning with X-ray crystal structures of the first DNA-binding domains to be studied, through recent cryo-EM structures of transcription factor binding to nucleosomal DNA.

Speed-Specificity Trade-Offs in the Transcription Factors Search for Their Genomic Binding Sites

Transcription factors (TFs) regulate gene expression by binding DNA sequences recognized by their DNA-binding domains (DBDs). DBD-recognized motifs are short and highly abundant in genomes. The ability of TFs to bind a specific subset of motif-containing sites, and to do so rapidly upon activation, is fundamental for gene expression in all eukaryotes. Despite extensive interest, our understanding of the TF-target search process is fragmented; although binding specificity and detection speed are two facets of this same process, trade-offs between them are rarely addressed. In this opinion article, we discuss potential speed-specificity trade-offs in the context of existing models. We further discuss the recently described 'distributed specificity' paradigm, suggesting that intrinsically disordered regions (IDRs) promote specificity while reducing the TF-target search time.

Functional and Regulatory Roles of Fold-Switching Proteins

Fold-switching proteins respond to cellular stimuli by remodeling their secondary structures and changing their functions. Whereas several previous reviews have focused on various structural, physical-chemical, and evolutionary aspects of this newly emerging class of proteins, this minireview focuses on how fold switching modulates protein function and regulates biological processes. It first compares and contrasts fold switchers with other known types of proteins. Second, it presents examples of how various proteins can change their functions through fold switching. Third, it demonstrates that fold switchers can regulate biological processes by discussing two proteins, RfaH and KaiB, whose dramatic secondary structure remodeling events directly affect gene expression and a circadian clock, respectively. Finally, this minireview discusses how the field of protein fold switching might advance.

Crystal Structures of Ternary Complexes of MEF2 and NKX2-5 Bound to DNA Reveal a Disease Related Protein-Protein Interaction Interface

MEF2 and NKX2-5 transcription factors interact with each other in cardiogenesis and are necessary for normal heart formation. Despite evidence suggesting that these two transcription factors function synergistically and possibly through direct physical interactions, molecular mechanisms by which they interact are not clear. Here we determined the crystal structures of ternary complexes of MEF2 and NKX2-5 bound to myocardin enhancer DNA in two crystal forms. These crystal structures are the first example of human MADS-box/homeobox ternary complex structures involved in cardiogenesis. Our structures reveal two possible modes of interactions between MEF2 and NKX2-5: MEF2 and NKX bind to adjacent DNA sites to recognize DNA in cis; and MEF2 and NKX bind to different DNA strands to interact with each other in trans via a conserved protein-protein interface observed in both crystal forms. Disease-related mutations are mapped to the observed protein-protein interface. Our structural studies provide a starting point to understand and further study the molecular mechanisms of the interactions between MEF2 and NKX2.5 and their roles in cardiogenesis.

Design of orthogonal regulatory systems for modulating gene expression in plants

Agricultural biotechnology strategies often require the precise regulation of multiple genes to effectively modify complex plant traits. However, most efforts are hindered by a lack of characterized tools that allow for reliable and targeted expression of transgenes. We have successfully engineered a library of synthetic transcriptional regulators that modulate expression strength in planta. By leveraging orthogonal regulatory systems from Saccharomyces spp., we have developed a strategy for the design of synthetic activators, synthetic repressors, and synthetic promoters and have validated their use in Nicotiana benthamiana and Arabidopsis thaliana. This characterization of contributing genetic elements that dictate gene expression represents a foundation for the rational design of refined synthetic regulators. Our findings demonstrate that these tools provide variation in transcriptional output while enabling the concerted expression of multiple genes in a tissue-specific and environmentally responsive manner, providing a basis for generating complex genetic circuits that process endogenous and environmental stimuli.

Protein-coding changes preceded cis-regulatory gains in a newly evolved transcription circuit

Changes in both the coding sequence of transcriptional regulators and in the cis-regulatory sequences recognized by these regulators have been implicated in the evolution of transcriptional circuits. However, little is known about how they evolved in concert. We describe an evolutionary pathway in fungi where a new transcriptional circuit (a-specific gene repression by the homeodomain protein Matα2) evolved by coding changes in this ancient regulator, followed millions of years later by cis-regulatory sequence changes in the genes of its future regulon. By analyzing a group of species that has acquired the coding changes but not the cis-regulatory sites, we show that the coding changes became necessary for the regulator's deeply conserved function, thereby poising the regulator to jump-start formation of the new circuit.

Experimental maps of DNA structure at nucleotide resolution distinguish intrinsic from protein-induced DNA deformations

Recognition of DNA by proteins depends on DNA sequence and structure. Often unanswered is whether the structure of naked DNA persists in a protein–DNA complex, or whether protein binding changes DNA shape. While X-ray structures of protein–DNA complexes are numerous, the structure of naked cognate DNA is seldom available experimentally. We present here an experimental and computational analysis pipeline that uses hydroxyl radical cleavage to map, at single-nucleotide resolution, DNA minor groove width, a recognition feature widely exploited by proteins. For 11 protein–DNA complexes, we compared experimental maps of naked DNA minor groove width with minor groove width measured from X-ray co-crystal structures. Seven sites had similar minor groove widths as naked DNA and when bound to protein. For four sites, part of the DNA in the complex had the same structure as naked DNA, and part changed structure upon protein binding. We compared the experimental map with minor groove patterns of DNA predicted by two computational approaches, DNAshape and ORChID2, and found good but not perfect concordance with both. This experimental approach will be useful in mapping structures of DNA sequences for which high-resolution structural data are unavailable. This approach allows probing of protein family-dependent readout mechanisms.

Structural Basis for Plant MADS Transcription Factor Oligomerization

MADS transcription factors (TFs) are DNA binding proteins found in almost all eukaryotes that play essential roles in diverse biological processes. While present in animals and fungi as a small TF family, the family has dramatically expanded in plants over the course of evolution, with the model flowering plant, Arabidopsis thaliana, possessing over 100 type I and type II MADS TFs. All MADS TFs contain a core and highly conserved DNA binding domain called the MADS or M domain. Plant MADS TFs have diversified this domain with plant-specific auxiliary domains. Plant type I MADS TFs have a highly diverse and largely unstructured Carboxy-terminal (C domain), whereas type II MADS have added oligomerization domains, called Intervening (I domain) and Keratin-like (K domain), in addition to the C domain. In this mini review, we describe the overall structure of the type II "MIKC" type MADS TFs in plants, with a focus on the K domain, a critical oligomerization module. We summarize the determining factors for oligomerization and provide mechanistic insights on how secondary structural elements are required for oligomerization capability and specificity. Using MADS TFs that are involved in flower organ specification as an example, we provide case studies and homology modeling of MADS TFs complex formation. Finally, we highlight outstanding questions in the field.

Chemical cross-linking in the structural analysis of protein assemblies

For decades, chemical cross-linking of proteins has been an established method to study protein interaction partners. The chemical cross-linking approach has recently been revived by mass spectrometric analysis of the cross-linking reaction products. Chemical cross-linking and mass spectrometric analysis (CXMS) enables the identification of residues that are close in three-dimensional (3D) space but not necessarily close in primary sequence. Therefore, this approach provides medium resolution information to guide de novo structure prediction, protein interface mapping and protein complex model building. The robustness and compatibility of the CXMS approach with multiple biochemical methods have made it especially appealing for challenging systems with multiple biochemical compositions and conformation states. This review provides an overview of the CXMS approach, describing general procedures in sample processing, data acquisition and analysis. Selection of proper chemical cross-linking reagents, strategies for cross-linked peptide identification, and successful application of CXMS in structural characterization of proteins and protein complexes are discussed.

Genome-wide determinants of sequence-specific DNA binding of general regulatory factors

General regulatory factors (GRFs), such as Reb1, Abf1, Rap1, Mcm1, and Cbf1, positionally organize yeast chromatin through interactions with a core consensus DNA sequence. It is assumed that sequence recognition via direct base readout suffices for specificity and that spurious nonfunctional sites are rendered inaccessible by chromatin. We tested these assumptions through genome-wide mapping of GRFs in vivo and in purified biochemical systems at near–base pair (bp) resolution using several ChIP-exo–based assays. We find that computationally predicted DNA shape features (e.g., minor groove width, helix twist, base roll, and propeller twist) that are not defined by a unique consensus sequence are embedded in the nonunique portions of GRF motifs and contribute critically to sequence-specific binding. This dual source specificity occurs at GRF sites in promoter regions where chromatin organization starts. Outside of promoter regions, strong consensus sites lack the shape component and consequently lack an intrinsic ability to bind cognate GRFs, without regard to influences from chromatin. However, sites having a weak consensus and low intrinsic affinity do exist in these regions but are rendered inaccessible in a chromatin environment. Thus, GRF site-specificity is achieved through integration of favorable DNA sequence and shape readouts in promoter regions and by chromatin-based exclusion from fortuitous weak sites within gene bodies. This study further revealed a severe G/C nucleotide cross-linking selectivity inherent in all formaldehyde-based ChIP assays, which includes ChIP-seq. However, for most tested proteins, G/C selectivity did not appreciably affect binding site detection, although it does place limits on the quantitativeness of occupancy levels.

The Cancer Mutation D83V Induces an α-Helix to β-Strand Conformation Switch in MEF2B

MEF2B is a major target of somatic mutations in non-Hodgkin lymphoma. Most of these mutations are non-synonymous substitutions of surface residues in the MADS-box/MEF2 domain. Among them, D83V is the most frequent mutation found in tumor cells. The link between this hotspot mutation and cancer is not well understood. Here we show that the D83V mutation induces a dramatic α-helix to β-strand switch in the MEF2 domain. Located in an α-helix region rich in β-branched residues, the D83V mutation not only removes the extensive helix stabilization interactions but also introduces an additional β-branched residue that further shifts the conformation equilibrium from α-helix to β-strand. Cross-database analyses of cancer mutations and chameleon sequences revealed a number of well-known cancer targets harboring β-strand favoring mutations in chameleon α-helices, suggesting a commonality of such conformational switch in certain cancers and a new factor to consider when stratifying the rapidly expanding cancer mutation data.

Genome-wide prediction of minor-groove electrostatic potential enables biophysical modeling of protein-DNA binding

Protein–DNA binding is a fundamental component of gene regulatory processes, but it is still not completely understood how proteins recognize their target sites in the genome. Besides hydrogen bonding in the major groove (base readout), proteins recognize minor-groove geometry using positively charged amino acids (shape readout). The underlying mechanism of DNA shape readout involves the correlation between minor-groove width and electrostatic potential (EP). To probe this biophysical effect directly, rather than using minor-groove width as an indirect measure for shape readout, we developed a methodology, DNAphi, for predicting EP in the minor groove and confirmed the direct role of EP in protein–DNA binding using massive sequencing data. The DNAphi method uses a sliding-window approach to mine results from non-linear Poisson–Boltzmann (NLPB) calculations on DNA structures derived from all-atom Monte Carlo simulations. We validated this approach, which only requires nucleotide sequence as input, based on direct comparison with NLPB calculations for available crystal structures. Using statistical machine-learning approaches, we showed that adding EP as a biophysical feature can improve the predictive power of quantitative binding specificity models across 27 transcription factor families. High-throughput prediction of EP offers a novel way to integrate biophysical and genomic studies of protein–DNA binding.

Cross-Family Transcription Factor Interactions: An Additional Layer of Gene Regulation

Specific and dynamic gene expression strongly depends on transcription factor (TF) activity and most plant TFs function in a combinatorial fashion. They can bind to DNA and control the expression of the corresponding gene in an additive fashion or cooperate by physical interactions, forming larger protein complexes. The importance of protein-protein interactions between members of a particular plant TF family has long been recognised; however, a significant number of interfamily TF interactions has recently been reported. The biological implications and the molecular mechanisms involved in cross-family interactions have now started to be elucidated and the examples illustrate potential roles in the bridging of biological processes. Hence, cross-family TF interactions expand the molecular toolbox for plants with additional mechanisms to control and fine-tune robust gene expression patterns and to adapt to their continuously changing environment.

A Model for Dimerization of the SOX Group E Transcription Factor Family

Group E members of the SOX transcription factor family include SOX8, SOX9, and SOX10. Preceding the high mobility group (HMG) domain in each of these proteins is a thirty-eight amino acid region that supports the formation of dimers on promoters containing tandemly inverted sites. The purpose of this study was to obtain new structural insights into how the dimerization region functions with the HMG domain. From a mutagenic scan of the dimerization region, the most essential amino acids of the dimerization region were clustered on the hydrophobic face of a single, predicted amphipathic helix. Consistent with our hypothesis that the dimerization region directly contacts the HMG domain, a peptide corresponding to the dimerization region bound a preassembled HMG-DNA complex. Sequence conservation among Group E members served as a basis to identify two surface exposed amino acids in the HMG domain of SOX9 that were necessary for dimerization. These data were combined to make a molecular model that places the dimerization region of one SOX9 protein onto the HMG domain of another SOX9 protein situated at the opposing site of a tandem promoter. The model provides a detailed foundation for assessing the impact of mutations on SOX Group E transcription factors.

Mating-type Gene Switching in Saccharomyces cerevisiae

The budding yeast Saccharomyces cerevisiae has two alternative mating types designated MATa and MATα. These are distinguished by about 700 bp of unique sequences, Ya or Yα, including divergent promoter sequences and part of the open reading frames of genes that regulate mating phenotype. Homothallic budding yeast, carrying an active HO endonuclease gene, HO, can switch mating type through a recombination process known as gene conversion, in which a site-specific double-strand break (DSB) created immediately adjacent to the Y region results in replacement of the Y sequences with a copy of the opposite mating type information, which is harbored in one of two heterochromatic donor loci, HMLα or HMRa. HO gene expression is tightly regulated to ensure that only half of the cells in a lineage switch to the opposite MAT allele, thus promoting conjugation and diploid formation. Study of the silencing of these loci has provided a great deal of information about the role of the Sir2 histone deacetylase and its associated Sir3 and Sir4 proteins in creating heterochromatic regions. MAT switching has been examined in great detail to learn about the steps in homologous recombination. MAT switching is remarkably directional, with MATa recombining preferentially with HMLα and MATα using HMRa. Donor preference is controlled by a cis-acting recombination enhancer located near HML. RE is turned off in MATα cells but in MATa binds multiple copies of the Fkh1 transcription factor whose forkhead-associated phosphothreonine binding domain localizes at the DSB, bringing HML into conjunction with MATa.

Degradation elements coincide with cofactor binding sites in a short-lived transcription factor

Elaborate control of gene expression by transcription factors is common to all kingdoms of life. In eukaryotes, transcription factor abundance and activity are often regulated by targeted proteolysis via the ubiquitin-proteasome system (UPS). The yeast MATα2 (α2) cell type regulator has long served as a model for UPS-dependent transcription factor degradation. Proteolysis of α2 is complex: it involves at least 2 ubiquitylation pathways and multiple regions of α2 affect its degradation. Such complexity also exists for the degradation of other UPS substrates. Here I review α2 degradation, most notably our recent identification of 2 novel degradation elements within α2 that overlap corepressor binding sites. I discuss possible implications of these findings and consider how principles of α2 proteolysis may be relevant to the degradation of other UPS substrates.

Evolutionary Dynamics of Floral Homeotic Transcription Factor Protein-Protein Interactions

Protein–protein interactions (PPIs) have widely acknowledged roles in the regulation of development, but few studies have addressed the timing and mechanism of shifting PPIs over evolutionary history. The B-class MADS-box transcription factors, PISTILLATA (PI) and APETALA3 (AP3) are key regulators of floral development. PI-like (PI^L) and AP3-like (AP3^L) proteins from a number of plants, including Arabidopsis thaliana (Arabidopsis) and the grass Zea mays (maize), bind DNA as obligate heterodimers. However, a PI^L protein from the grass relative Joinvillea can bind DNA as a homodimer. To ascertain whether Joinvillea PI^L homodimerization is an anomaly or indicative of broader trends, we characterized PI^L dimerization across the Poales and uncovered unexpected evolutionary lability. Both obligate B-class heterodimerization and PI^L homodimerization have evolved multiple times in the order, by distinct molecular mechanisms. For example, obligate B-class heterodimerization in maize evolved very recently from PI^L homodimerization. A single amino acid change, fixed during domestication, is sufficient to toggle one maize PI^L protein between homodimerization and obligate heterodimerization. We detected a signature of positive selection acting on residues preferentially clustered in predicted sites of contact between MADS-box monomers and dimers, and in motifs that mediate MADS PPI specificity in Arabidopsis. Changing one positively selected residue can alter PI^L dimerization activity. Furthermore, ectopic expression of a Joinvillea PI^L homodimer in Arabidopsis can homeotically transform sepals into petals. Our results provide a window into the evolutionary remodeling of PPIs, and show that novel interactions have the potential to alter plant form in a context-dependent manner.

Grid-based backbone correction to the ff12SB protein force field for implicit-solvent simulations

Force fields, such as Amber's ff12SB, can be fairly accurate models of the physical forces in proteins and other biomolecules. When coupled with accurate solvation models, force fields are able to bring insight into the conformational preferences, transitions, pathways, and free energies for these biomolecules. When computational speed/cost matters, implicit solvent is often used but at the cost of accuracy. We present an empirical grid-like correction term, in the spirit of cMAPs, to the combination of the ff12SB protein force field and the GBneck2 implicit-solvent model. Ff12SB-cMAP is parametrized on experimental helicity data. We provide validation on a set of peptides and proteins. Ff12SB-cMAP successfully improves the secondary structure biases observed in ff12SB + Gbneck2. Ff12SB-cMAP can be downloaded ( https://github.com/laufercenter/Amap.git ) and used within the Amber package. It can improve the agreement of force fields + implicit solvent with experiments.

STUbL-mediated degradation of the transcription factor MATα2 requires degradation elements that coincide with corepressor binding sites

The yeast cell type regulator MATα (α2) is degraded through two ubiquitylation pathways, one of which has been minimally characterized. We identify two regions in α2 important for this pathway and show that these regions overlap specific binding sites for α2 corepressors, suggesting that α2 degradation is coordinated with its functional status.

The yeast transcription factor MATα2 (α2) is a short-lived protein known to be ubiquitylated by two distinct pathways, one involving the ubiquitin-conjugating enzymes (E2s) Ubc6 and Ubc7 and the ubiquitin ligase (E3) Doa10 and the other operating with the E2 Ubc4 and the heterodimeric E3 Slx5/Slx8. Although Slx5/Slx8 is a small ubiquitin-like modifier (SUMO)-targeted ubiquitin ligase (STUbL), it does not require SUMO to target α2 but instead directly recognizes α2. Little is known about the α2 determinants required for its Ubc4- and STUbL-mediated degradation or how these determinants substitute for SUMO in recognition by the STUbL pathway. We describe two distinct degradation elements within α2, both of which are necessary for α2 recognition specifically by the Ubc4 pathway. Slx5/Slx8 can directly ubiquitylate a C-terminal fragment of α2, and mutating one of the degradation elements impairs this ubiquitylation. Surprisingly, both degradation elements identified here overlap specific interaction sites for α2 corepressors: the Mcm1 interaction site in the central α2 linker and the Ssn6 (Cyc8) binding site in the α2 homeodomain. We propose that competitive binding to α2 by the ubiquitylation machinery and α2 cofactors is balanced so that α2 can function in transcription repression yet be short lived enough to allow cell-type switching.

Decoupling of divergent gene regulation by sequence-specific DNA binding factors

Divergent gene pairs (DGPs) are abundant in eukaryotic genomes. Since two genes in a DGP potentially share the same regulatory sequence, one might expect that they should be co-regulated. However, an inspection of yeast DGPs containing cell-cycle or stress response genes revealed that most DGPs are differentially-regulated. The mechanism underlying DGP differential regulation is not understood. Here, we showed that co- versus differential regulation cannot be explained by genetic features including promoter length, binding site orientation, TATA elements, nucleosome distribution, or presence of non-coding RNAs. Using time-lapse fluorescence microscopy, we carried out an in-depth study of a differentially regulated DGP, PFK26-MOB1. We found that their differential regulation is mainly achieved through two DNA-binding factors, Tbf1 and Mcm1. Similar to 'enhancer-blocking insulators' in higher eukaryotes, these factors shield the proximal promoter from the action of more distant transcription regulators. We confirmed the blockage function of Tbf1 using synthetic promoters. We further presented evidence that the blockage mechanism is widely used among genome-wide DGPs. Besides elucidating the DGP regulatory mechanism, our work revealed a novel class of insulators in yeast.

Organization of the MADS box from human SRF revealed by tyrosine perturbation

MADS box family transcription factors are involved in signal transduction and development control through DNA specific sequence recognition. The DNA binding domain of these proteins contains a conservative 55-60 amino acid sequence which defines the membership of this large family. Here we present a thorough study of the MADS segment of serum response factor (MADS(SRF)). Fluorescence, UV-absorption, and Raman spectroscopy studies were performed in order to disclose its behavior and basic functional properties in an aqueous environment. The secondary structure of MADS(SRF) estimated by analysis of Raman spectra and supported by CD has revealed only the C-terminal part as homologous with those of free core-SRF, while the N-terminal part has lost the stable α-helical structure found in both the free core-SRF and its specific complex with DNA. The three tyrosine residues of the MADS(SRF) were used as spectroscopic inner probes. The effect of environmental conditions, especially pH variations and addition of variously charged quenchers, on their spectra was examined. Two-component fluorescence quenching was revealed using factor analysis and corresponding Stern-Volmer constants determined. Factor analysis of absorbance and fluorescence pH titration led to determination of three dissociation constants pKa1 = 6.4 ± 0.2, pKa2 = 7.3 ± 0.2, and pKa3 = 9.6 ± 0.6. Critical comparison of all experiments identified the deprotonation of His193 hydrogen bonded to Tyr195 as a candidate for pKa1 (and that of Tyr158 as a candidate for pKa2). Within MADS(SRF), His193 is a key intermediary between the N-terminal primary DNA binding element and the hydrophobic C-terminal protein dimerization element.

Structural basis for the oligomerization of the MADS domain transcription factor SEPALLATA3 in Arabidopsis

In plants, MADS domain transcription factors act as central regulators of diverse developmental pathways. In Arabidopsis thaliana, one of the most central members of this family is SEPALLATA3 (SEP3), which is involved in many aspects of plant reproduction, including floral meristem and floral organ development. SEP3 has been shown to form homo and heterooligomeric complexes with other MADS domain transcription factors through its intervening (I) and keratin-like (K) domains. SEP3 function depends on its ability to form specific protein-protein complexes; however, the atomic level determinants of oligomerization are poorly understood. Here, we report the 2.5-Å crystal structure of a small portion of the intervening and the complete keratin-like domain of SEP3. The domains form two amphipathic alpha helices separated by a rigid kink, which prevents intramolecular association and presents separate dimerization and tetramerization interfaces comprising predominantly hydrophobic patches. Mutations to the tetramerization interface demonstrate the importance of highly conserved hydrophobic residues for tetramer stability. Atomic force microscopy was used to show SEP3-DNA interactions and the role of oligomerization in DNA binding and conformation. Based on these data, the oligomerization patterns of the larger family of MADS domain transcription factors can be predicted and manipulated based on the primary sequence.

Virtual screening on an α-helix to β-strand switchable region of the FGFR2 extracellular domain revealed positive and negative modulators

The secondary structure of some protein segments may vary between α-helix and β-strand. To predict these switchable segments, we have developed an algorithm, Switch-P, based solely on the protein sequence. This algorithm was used on the extracellular parts of FGF receptors. For FGFR2, it predicted that β4 and β5 strands of the third Ig-like domain were highly switchable. These two strands possess a high number of somatic mutations associated with cancer. Analysis of PDB structures of FGF receptors confirmed the switchability prediction for β5. We thus evaluated if compound-driven α-helix/β-strand switching of β5 could modulate FGFR2 signaling. We performed the virtual screening of a library containing 1.4 million of chemical compounds with two models of the third Ig-like domain of FGFR2 showing different secondary structures for β5, and we selected 32 compounds. Experimental testing using proliferation assays with FGF7-stimulated SNU-16 cells and a FGFR2-dependent Erk1/2 phosphorylation assay with FGFR2-transfected L6 cells, revealed activators and inhibitors of FGFR2. Our method for the identification of switchable proteinic regions, associated with our virtual screening approach, provides an opportunity to discover new generation of drugs with under-explored mechanism of action.

Mapping yeast transcriptional networks

The term “transcriptional network” refers to the mechanism(s) that underlies coordinated expression of genes, typically involving transcription factors (TFs) binding to the promoters of multiple genes, and individual genes controlled by multiple TFs. A multitude of studies in the last two decades have aimed to map and characterize transcriptional networks in the yeast Saccharomyces cerevisiae. We review the methodologies and accomplishments of these studies, as well as challenges we now face. For most yeast TFs, data have been collected on their sequence preferences, in vivo promoter occupancy, and gene expression profiles in deletion mutants. These systematic studies have led to the identification of new regulators of numerous cellular functions and shed light on the overall organization of yeast gene regulation. However, many yeast TFs appear to be inactive under standard laboratory growth conditions, and many of the available data were collected using techniques that have since been improved. Perhaps as a consequence, comprehensive and accurate mapping among TF sequence preferences, promoter binding, and gene expression remains an open challenge. We propose that the time is ripe for renewed systematic efforts toward a complete mapping of yeast transcriptional regulatory mechanisms.

Structural determinants of DNA recognition by plant MADS-domain transcription factors

Plant MADS-domain transcription factors act as key regulators of many developmental processes. Despite the wealth of information that exists about these factors, the mechanisms by which they recognize their cognate DNA-binding site, called CArG-box (consensus CCW₆GG), and how different MADS-domain proteins achieve DNA-binding specificity, are still largely unknown. We used information from in vivo ChIP-seq experiments, in vitro DNA-binding data and evolutionary conservation to address these important questions. We found that structural characteristics of the DNA play an important role in the DNA binding of plant MADS-domain proteins. The central region of the CArG-box largely resembles a structural motif called ‘A-tract’, which is characterized by a narrow minor groove and may assist bending of the DNA by MADS-domain proteins. Periodically spaced A-tracts outside the CArG-box suggest additional roles for this structure in the process of DNA binding of these transcription factors. Structural characteristics of the CArG-box not only play an important role in DNA-binding site recognition of MADS-domain proteins, but also partly explain differences in DNA-binding specificity of different members of this transcription factor family and their heteromeric complexes.

Solvent-accessible surface area: How well can be applied to hot-spot detection?

A detailed comprehension of protein-based interfaces is essential for the rational drug development. One of the key features of these interfaces is their solvent accessible surface area profile. With that in mind, we tested a group of 12 SASA-based features for their ability to correlate and differentiate hot- and null-spots. These were tested in three different data sets, explicit water MD, implicit water MD, and static PDB structure. We found no discernible improvement with the use of more comprehensive data sets obtained from molecular dynamics. The features tested were shown to be capable of discerning between hot- and null-spots, while presenting low correlations. Residue standardization such as rel SASAi or rel/res SASAi , improved the features as a tool to predict ΔΔGbinding values. A new method using support machine learning algorithms was developed: SBHD (Sasa-Based Hot-spot Detection). This method presents a precision, recall, and F1 score of 0.72, 0.81, and 0.76 for the training set and 0.91, 0.73, and 0.81 for an independent test set.

Following gene duplication, paralog interference constrains transcriptional circuit evolution

Most models of gene duplication assume that the ancestral functions of the preduplication gene are independent and can therefore be neatly partitioned between descendant paralogs. However, many gene products, such as transcriptional regulators, are components within cooperative assemblies; here, we show that a natural consequence of duplication and divergence of such proteins can be competitive interference between the paralogs. Our example is based on the duplication of the essential MADS-box transcriptional regulator Mcm1, which is found in all fungi and regulates a large set of genes. We show that a set of historical amino acid sequence substitutions minimized paralog interference in contemporary species and, in doing so, increased the molecular complexity of this gene regulatory network. We propose that paralog interference is a common constraint on gene duplicate evolution, and its resolution, which can generate additional regulatory complexity, is needed to stabilize duplicated genes in the genome.

Crystallization studies of the keratin-like domain from Arabidopsis thaliana SEPALLATA 3

In higher plants, the MADS-box genes encode a large family of transcription factors (TFs) involved in key developmental processes, most notably plant reproduction, flowering and floral organ development. SEPALLATA 3 (SEP3) is a member of the MADS TF family and plays a role in the development of the floral organs through the formation of multiprotein complexes with other MADS-family TFs. SEP3 is divided into four domains: the M (MADS) domain, involved in DNA binding and dimerization, the I (intervening) domain, a short domain involved in dimerization, the K (keratin-like) domain important for multimeric MADS complex formation and the C (C-terminal) domain, a largely unstructured region putatively important for higher-order complex formation. The entire K domain along with a portion of the I and C domains of SEP3 was crystallized using high-throughput robotic screening followed by optimization. The crystals belonged to space group P2(1)2(1)2, with unit-cell parameters a = 123.44, b = 143.07, c = 49.83 Å, and a complete data set was collected to 2.53 Å resolution.

Computational Alanine Scanning Mutagenesis-An Improved Methodological Approach for Protein-DNA Complexes

Proteins and protein-based complexes are the basis of many key systems in nature and have been the subject of intense research in the last decades, in an attempt to acquire comprehensive knowledge of reactions that take place in nature. Computational Alanine Scanning Mutagenesis approaches have been extensively used in the study of protein interfaces and in the determination of the most important residues for complex formation, the Hot-spots. However, as it is usually applied to the study of protein-protein interfaces, we tried to modify and apply it to the study of protein-DNA interfaces, which are also crucial in nature but have not been the subject of as much research. In this work, we carry out MD simulations of seven protein-DNA complexes and tested the influence of the variation of different parameters on the determination of the binding free energy terms (ΔΔGbinding) of 78 mutations: solvent representation, internal dielectric constant, Linear and Nonlinear Poisson-Boltzmann equation, Generalized Born model, simulation time, number of structures analyzed, number of MD trajectories, force field used, and energetic terms involved. Overall, this new approach gave an average error of 1.55 kcal/mol, and P, R, F1, accuracy, and specificity values of 0.78, 0.50, 0.61, 0.77, and 0.92, respectively. This improved computational alanine scanning mutagenesis approach may serve as a tool to explore the behavior of this important class of complexes.

Cloning and functional analysis of the flowering gene GmSOC1-like, a putative SUPPRESSOR OF OVEREXPRESSION CO1/AGAMOUS-LIKE 20 (SOC1/AGL20) ortholog in soybean

KEY MESSAGE

The major insight in this manuscript is that we identified a new flowering regulator, GmSOC1-like, which may participate in the initiation and maintenance of flowering in soybean. Flowering is pivotal for the reproductive behavior of plants, and it is regulated by complex and coordinated genetic networks that are fine-tuned by endogenous cues and environmental signals. To better understand the molecular basis of flowering regulation in soybean, we isolated GmSOC1 and GmSOC1-like, two putative soybean orthologs for the Arabidopsis SUPPRESSOR OF OVEREXPRESSION OF CO1/AGAMOUS-LIKE 20 (SOC1/AGL20). The expression pattern of GmSOC1-like was analyzed by qRT-PCR in Zigongdongdou, a photoperiod-sensitive soybean cultivar. GmSOC1-like was widely expressed at different levels in most organs of the soybean, with the highest expression in the shoot apex during the early stage of floral transition. In addition, its expression showed a circadian rhythm pattern, with the highest expression at midnight under short-day (SD) condition. Intriguingly, GmSOC1-like was induced 4 days earlier than GmSOC1 during flowering transition in SD, suggesting that GmSOC1 and GmSOC1-like expression might be differentially regulated. However, under long-day (LD) condition, the expression of GmSOC1 and GmSOC1-like decreased gradually in the shoot apex of Zigongdongdou, which is in accordance with the fact that Zigongdongdou maintains vegetative growth in LD. In addition, overexpression of GmSOC1-like stimulated the flowering of Lotus corniculatus cv. supperroot plants. In conclusion, the results of this study indicate that GmSOC1-like may act as a flowering inducer in soybean.

Bck2 acts through the MADS box protein Mcm1 to activate cell-cycle-regulated genes in budding yeast

The Bck2 protein is a potent genetic regulator of cell-cycle-dependent gene expression in budding yeast. To date, most experiments have focused on assessing a potential role for Bck2 in activation of the G1/S-specific transcription factors SBF (Swi4, Swi6) and MBF (Mbp1, Swi6), yet the mechanism of gene activation by Bck2 has remained obscure. We performed a yeast two-hybrid screen using a truncated version of Bck2 and discovered six novel Bck2-binding partners including Mcm1, an essential protein that binds to and activates M/G1 promoters through Early Cell cycle Box (ECB) elements as well as to G2/M promoters. At M/G1 promoters Mcm1 is inhibited by association with two repressors, Yox1 or Yhp1, and gene activation ensues once repression is relieved by an unknown activating signal. Here, we show that Bck2 interacts physically with Mcm1 to activate genes during G1 phase. We used chromatin immunoprecipitation (ChIP) experiments to show that Bck2 localizes to the promoters of M/G1-specific genes, in a manner dependent on functional ECB elements, as well as to the promoters of G1/S and G2/M genes. The Bck2-Mcm1 interaction requires valine 69 on Mcm1, a residue known to be required for interaction with Yox1. Overexpression of BCK2 decreases Yox1 localization to the early G1-specific CLN3 promoter and rescues the lethality caused by overexpression of YOX1. Our data suggest that Yox1 and Bck2 may compete for access to the Mcm1-ECB scaffold to ensure appropriate activation of the initial suite of genes required for cell cycle commitment.

Local conformational changes in the DNA interfaces of proteins

When a protein binds to DNA, a conformational change is often induced so that the protein will fit into the DNA structure. Therefore, quantitative analyses were conducted to understand the conformational changes in proteins. The results showed that conformational changes in DNA interfaces are more frequent than in non-interfaces, and DNA interfaces have more conformational variations in the DNA-free form. As expected, the former indicates that interaction with DNA has some influence on protein structure. The latter suggests that the intrinsic conformational flexibility of DNA interfaces is important for adjusting their conformation for DNA. The amino acid propensities of the conformationally changed regions in DNA interfaces indicate that hydrophilic residues are preferred over the amino acids that appear in the conformationally unchanged regions. This trend is true for disordered regions, suggesting again that intrinsic flexibility is of importance not only for DNA binding but also for interactions with other molecules. These results demonstrate that fragments destined to be DNA interfaces have an intrinsic flexibility and are composed of amino acids with the capability of binding to DNA. This information suggests that the prediction of DNA binding sites may be improved by the integration of amino acid preference for DNA and one for disordered regions.

Protein modularity, cooperative binding, and hybrid regulatory states underlie transcriptional network diversification

We examine how different transcriptional network structures can evolve from an ancestral network. By characterizing how the ancestral mode of gene regulation for genes specific to a-type cells in yeast species evolved from an activating paradigm to a repressing one, we show that regulatory protein modularity, conversion of one cis-regulatory sequence to another, distribution of binding energy among protein-protein and protein-DNA interactions, and exploitation of ancestral network features all contribute to the evolution of a novel regulatory mode. The formation of this derived mode of regulation did not disrupt the ancestral mode and thereby created a hybrid regulatory state where both means of transcription regulation (ancestral and derived) contribute to the conserved expression pattern of the network. Finally, we show how this hybrid regulatory state has resolved in different ways in different lineages to generate the diversity of regulatory network structures observed in modern species.

Mating-type genes and MAT switching in Saccharomyces cerevisiae

Mating type in Saccharomyces cerevisiae is determined by two nonhomologous alleles, MATa and MATα. These sequences encode regulators of the two different haploid mating types and of the diploids formed by their conjugation. Analysis of the MATa1, MATα1, and MATα2 alleles provided one of the earliest models of cell-type specification by transcriptional activators and repressors. Remarkably, homothallic yeast cells can switch their mating type as often as every generation by a highly choreographed, site-specific homologous recombination event that replaces one MAT allele with different DNA sequences encoding the opposite MAT allele. This replacement process involves the participation of two intact but unexpressed copies of mating-type information at the heterochromatic loci, HMLα and HMRa, which are located at opposite ends of the same chromosome-encoding MAT. The study of MAT switching has yielded important insights into the control of cell lineage, the silencing of gene expression, the formation of heterochromatin, and the regulation of accessibility of the donor sequences. Real-time analysis of MAT switching has provided the most detailed description of the molecular events that occur during the homologous recombinational repair of a programmed double-strand chromosome break.

Transcriptional regulation in Saccharomyces cerevisiae: transcription factor regulation and function, mechanisms of initiation, and roles of activators and coactivators

Here we review recent advances in understanding the regulation of mRNA synthesis in Saccharomyces cerevisiae. Many fundamental gene regulatory mechanisms have been conserved in all eukaryotes, and budding yeast has been at the forefront in the discovery and dissection of these conserved mechanisms. Topics covered include upstream activation sequence and promoter structure, transcription factor classification, and examples of regulated transcription factor activity. We also examine advances in understanding the RNA polymerase II transcription machinery, conserved coactivator complexes, transcription activation domains, and the cooperation of these factors in gene regulatory mechanisms.

Classification of proteins: available structural space for molecular modeling

The wealth of available protein structural data provides unprecedented opportunity to study and better understand the underlying principles of protein folding and protein structure evolution. A key to achieving this lies in the ability to analyse these data and to organize them in a coherent classification scheme. Over the past years several protein classifications have been developed that aim to group proteins based on their structural relationships. Some of these classification schemes explore the concept of structural neighbourhood (structural continuum), whereas other utilize the notion of protein evolution and thus provide a discrete rather than continuum view of protein structure space. This chapter presents a strategy for classification of proteins with known three-dimensional structure. Steps in the classification process along with basic definitions are introduced. Examples illustrating some fundamental concepts of protein folding and evolution with a special focus on the exceptions to them are presented.

A DNA minor groove electronegative potential genome map based on photo-chemical probing

The double-stranded DNA of the genome contains both sequence information directly relating to the protein and RNA coding as well as functional and structural information relating to protein recognition. Only recently is the importance of DNA shape in this recognition process being fully appreciated, and it also appears that minor groove electronegative potential may contribute significantly in guiding proteins to their cognate binding sites in the genome. Based on the photo-chemical probing results, we have derived an algorithm that predicts the minor groove electronegative potential in a DNA helix of any given sequence. We have validated this model on a series of protein-DNA binding sites known to involve minor groove electrostatic recognition as well as on stable nucleosome core complexes. The algorithm allows for the first time a full minor groove electrostatic description at the nucleotide resolution of any genome, and it is illustrated how such detailed studies of this sequence dependent, inherent property of the DNA may reflect on genome organization, gene expression and chromosomal condensation.