Dynamic protein-DNA recognition: beyond what can be seen

Context-dependent, fuzzy protein interactions: Towards sequence-based insights

Predicting protein interactions in the cellular environment still remains a challenge in the AlphaFold era. Protein interactions, similarly to their structures, sample a continuum from ordered to disordered states, with specific partners in many bound configurations. A multiplicity of binding modes (MBM) enables transition between these states under different cellular conditions. This review focuses on how the cellular environment affects protein interactions, highlighting the molecular mechanisms, biophysical origin, and sequence-based principles of context-dependent, fuzzy interactions. It summarises experimental and computational approaches to address the challenge of interaction heterogeneity and its contribution to a wide range of biological functions. These insights will help in understanding complex cellular processes, involving conversions between protein assembly states, such as from liquid-like droplet state to the amyloid state.

Fuzzy recognition by the prokaryotic transcription factor HigA2 from Vibrio cholerae

Disordered protein sequences can exhibit different binding modes, ranging from well-ordered folding-upon-binding to highly dynamic fuzzy binding. The primary function of the intrinsically disordered region of the antitoxin HigA2 from Vibrio cholerae is to neutralize HigB2 toxin through ultra-high-affinity folding-upon-binding interaction. Here, we show that the same intrinsically disordered region can also mediate fuzzy interactions with its operator DNA and, through interplay with the folded helix-turn-helix domain, regulates transcription from the higBA2 operon. NMR, SAXS, ITC and in vivo experiments converge towards a consistent picture where a specific set of residues in the intrinsically disordered region mediate electrostatic and hydrophobic interactions while "hovering" over the DNA operator. Sensitivity of the intrinsically disordered region to scrambling the sequence, position-specific contacts and absence of redundant, multivalent interactions, point towards a more specific type of fuzzy binding. Our work demonstrates how a bacterial regulator achieves dual functionality by utilizing two distinct interaction modes within the same disordered sequence.

Mechanistic regulation of FOXO transcription factors in the nucleus

FOXO proteins represent evolutionarily conserved transcription factors (TFs) that play critical roles in responding to various physiological signals or pathological stimuli, either through transcription-dependent or -independent mechanisms. Dysfunction of these proteins have been implicated in numerous diseases, including cancer. Although the regulation of FOXO TFs shuttling between the cytoplasm and the nucleus has been extensively studied and reviewed, there's still a lack of a comprehensive review focusing on the intricate interactions between FOXO, DNA, and cofactors in the regulation of gene expression. In this review, we aim to summarize recent advances and provide a detailed understanding of the mechanism underlying FOXO proteins binding to target DNA. Additionally, we will discuss the challenges associated with pharmacological approaches in modulating FOXO function, and explore the dynamic association between TF, DNA, and RNA on chromatin. This review will contribute to a better understanding of mechanistic regulations of eukaryotic TFs within the nucleus.

A guideline for the distance measurement plans of site-directed spin labels for structural prediction of nucleic acids

CONTEXT AND RESULTS

Site-directed spin labeling (SDSL) combined with electron paramagnetic resonance spectroscopy methods has been successfully used to predict the structures of nucleic acids. These methods measure the distances between spin labels yielding distance equations that are solved using numerical algorithms to provide one or several structural predictions. In this work, the minimum number of SDSL distance measurements and distance measurement types required to predict a unique nucleic acid structure were investigated. Our results indicate that at least six distance measurements should be obtained given that the distance measurements do not connect one SDSL on one arm with more than three SDSLs on the other arm. Moreover, there may be a preference for 1-to-1 SLs distance measurements rather than 1-to-many SLs as the latter was linked to undefined structures discussed in this study.

METHODS

Pairs of double-helical arms of nucleic acid were simulated using the finite element software Pro/ENGINEER (PTC Inc., Boston, MA). In each simulation, a specific SDSL distance measurement plan was adopted and the resulting structure was tested for movability. Immovable structures indicate that this plan will potentially result in a unique structural prediction of the nucleic acid. All the possible plans for SDSL distance measurements were investigated either by direct measurement or by extrapolation.

Spatiotemporal and global profiling of DNA-protein interactions enables discovery of low-affinity transcription factors

Precise dissection of DNA-protein interactions is essential for elucidating the recognition basis, dynamics and gene regulation mechanism. However, global profiling of weak and dynamic DNA-protein interactions remains a long-standing challenge. Here, we establish the light-induced lysine (K) enabled crosslinking (LIKE-XL) strategy for spatiotemporal and global profiling of DNA-protein interactions. Harnessing unique abilities to capture weak and transient DNA-protein interactions, we demonstrate that LIKE-XL enables the discovery of low-affinity transcription-factor/DNA interactions via sequence-specific DNA baits, determining the binding sites for transcription factors that have been previously unknown. More importantly, we successfully decipher the dynamics of the transcription factor subproteome in response to drug treatment in a time-resolved manner, and find downstream target transcription factors from drug perturbations, providing insight into their dynamic transcriptional networks. The LIKE-XL strategy offers a complementary method to expand the DNA-protein profiling toolbox and map accurate DNA-protein interactomes that were previously inaccessible via non-covalent strategies, for better understanding of protein function in health and disease.

Small-angle x-ray and neutron scattering of MexR and its complex with DNA supports a conformational selection binding model

In this work, we used small-angle x-ray and neutron scattering to reveal the shape of the protein-DNA complex of the Pseudomonas aeruginosa transcriptional regulator MexR, a member of the multiple antibiotics resistance regulator (MarR) family, when bound to one of its native DNA binding sites. Several MarR-like proteins, including MexR, repress the expression of efflux pump proteins by binding to DNA on regulatory sites overlapping with promoter regions. When expressed, efflux proteins self-assemble to form multiprotein complexes and actively expel highly toxic compounds out of the host organism. The mutational pressure on efflux-regulating MarR family proteins is high since deficient DNA binding leads to constitutive expression of efflux pumps and thereby supports acquired multidrug resistance. Understanding the functional outcome of such mutations and their effects on DNA binding has been hampered by the scarcity of structural and dynamic characterization of both free and DNA-bound MarR proteins. Here, we show how combined neutron and x-ray small-angle scattering of both states in solution support a conformational selection model that enhances MexR asymmetry in binding to one of its promoter-overlapping DNA binding sites.

Correlating multi-functional role of cold shock domain proteins with intrinsically disordered regions

Cold shock proteins (CSPs) are an ancient and conserved family of proteins. They are renowned for their role in response to low-temperature stress in bacteria and nucleic acid binding activities. In prokaryotes, cold and non-cold inducible CSPs are involved in various cellular and metabolic processes such as growth and development, osmotic oxidation, starvation, stress tolerance, and host cell invasion. In prokaryotes, cold shock condition reduces cell transcription and translation efficiency. Eukaryotic cold shock domain (CSD) proteins are evolved form of prokaryotic CSPs where CSD is flanked by N- and C-terminal domains. Eukaryotic CSPs are multi-functional proteins. CSPs also act as nucleic acid chaperons by preventing the formation of secondary structures in mRNA at low temperatures. In human, CSD proteins play a crucial role in the progression of breast cancer, colon cancer, lung cancer, and Alzheimer's disease. A well-defined three-dimensional structure of intrinsically disordered regions of CSPs family members is still undetermined. In this article, intrinsic disorder regions of CSPs have been explored systematically to understand the pleiotropic role of the cold shock family of proteins.

Intrinsically disordered proteins play diverse roles in cell signaling

Abstract

Signaling pathways allow cells to detect and respond to a wide variety of chemical (e.g. Ca²⁺ or chemokine proteins) and physical stimuli (e.g., sheer stress, light). Together, these pathways form an extensive communication network that regulates basic cell activities and coordinates the function of multiple cells or tissues. The process of cell signaling imposes many demands on the proteins that comprise these pathways, including the abilities to form active and inactive states, and to engage in multiple protein interactions. Furthermore, successful signaling often requires amplifying the signal, regulating or tuning the response to the signal, combining information sourced from multiple pathways, all while ensuring fidelity of the process. This sensitivity, adaptability, and tunability are possible, in part, due to the inclusion of intrinsically disordered regions in many proteins involved in cell signaling. The goal of this collection is to highlight the many roles of intrinsic disorder in cell signaling. Following an overview of resources that can be used to study intrinsically disordered proteins, this review highlights the critical role of intrinsically disordered proteins for signaling in widely diverse organisms (animals, plants, bacteria, fungi), in every category of cell signaling pathway (autocrine, juxtacrine, intracrine, paracrine, and endocrine) and at each stage (ligand, receptor, transducer, effector, terminator) in the cell signaling process. Thus, a cell signaling pathway cannot be fully described without understanding how intrinsically disordered protein regions contribute to its function. The ubiquitous presence of intrinsic disorder in different stages of diverse cell signaling pathways suggest that more mechanisms by which disorder modulates intra- and inter-cell signals remain to be discovered.

Graphical abstract

On the Effects of Disordered Tails, Supertertiary Structure and Quinary Interactions on the Folding and Function of Protein Domains

The vast majority of our current knowledge about the biochemical and biophysical properties of proteins derives from in vitro studies conducted on isolated globular domains. However, a very large fraction of the proteins expressed in the eukaryotic cell are structurally more complex. In particular, the discovery that up to 40% of the eukaryotic proteins are intrinsically disordered, or possess intrinsically disordered regions, and are highly dynamic entities lacking a well-defined three-dimensional structure, revolutionized the structure–function paradigm and our understanding of proteins. Moreover, proteins are mostly characterized by the presence of multiple domains, influencing each other by intramolecular interactions. Furthermore, proteins exert their function in a crowded intracellular milieu, transiently interacting with a myriad of other macromolecules. In this review we summarize the literature tackling these themes from both the theoretical and experimental perspectives, highlighting the effects on protein folding and function that are played by (i) flanking disordered tails; (ii) contiguous protein domains; (iii) interactions with the cellular environment, defined as quinary structures. We show that, in many cases, both the folding and function of protein domains is remarkably perturbed by the presence of these interactions, pinpointing the importance to increase the level of complexity of the experimental work and to extend the efforts to characterize protein domains in more complex contexts.

Tenuous transcriptional threshold of human sex determination. II. SRY exploits water-mediated clamp at the edge of ambiguity

Y-encoded transcription factor SRY initiates male differentiation in therian mammals. This factor contains a high-mobility-group (HMG) box, which mediates sequence-specific DNA binding with sharp DNA bending. A companion article in this issue described sex-reversal mutations at box position 72 (residue 127 in human SRY), invariant as Tyr among mammalian orthologs. Although not contacting DNA, the aromatic ring seals the domain's minor wing at a solvent-exposed junction with a basic tail. A seeming paradox was posed by the native-like biochemical properties of inherited Swyer variant Y72F: its near-native gene-regulatory activity is consistent with the father's male development, but at odds with the daughter's XY female somatic phenotype. Surprisingly, aromatic rings (Y72, F72 or W72) confer higher transcriptional activity than do basic or polar side chains generally observed at solvated DNA interfaces (Arg, Lys, His or Gln). Whereas biophysical studies (time-resolved fluorescence resonance energy transfer and heteronuclear NMR spectroscopy) uncovered only subtle perturbations, dissociation of the Y72F complex was markedly accelerated relative to wild-type. Studies of protein-DNA solvation by molecular-dynamics (MD) simulations of an homologous high-resolution crystal structure (SOX18) suggest that Y72 para-OH anchors a network of water molecules at the tail-DNA interface, perturbed in the variant in association with nonlocal conformational fluctuations. Loss of the Y72 anchor among SRY variants presumably "unclamps" its basic tail, leading to (a) rapid DNA dissociation despite native affinity and (b) attenuated transcriptional activity at the edge of sexual ambiguity. Conservation of Y72 suggests that this water-mediated clamp operates generally among SRY and metazoan SOX domains.

Fuzzy binding model of molecular interactions between tropoelastin and integrin alphaVbeta3

Tropoelastin is the highly flexible monomer subunit of elastin, required for the resilience of the extracellular matrix in elastic tissues. To elicit biological signaling, multiple sites on tropoelastin bind to cell surface integrins in a poorly understood multifactorial process. We constructed a full atomistic molecular model of the interactions between tropoelastin and integrin α_vβ₃ using ensemble-based computational methodologies. Conformational changes of integrin α_vβ₃ associated with outside-in signaling were more frequently facilitated in an ensemble in which tropoelastin bound the integrin's α1 helix rather than the upstream canonical binding site. Our findings support a model of fuzzy binding, whereby many tropoelastin conformations and defined sites cooperatively interact with multiple α_vβ₃ regions. This model explains prior experimental binding to distinct tropoelastin regions, domains 17 and 36, and points to the cooperative participation of domain 20. Our study highlights the utility of ensemble-based approaches in helping to understand the interactive mechanisms of functionally significant flexible proteins.

Phosphorus and sulfur SAD phasing of the nucleic acid-bound DNA-binding domain of interferon regulatory factor 4

Solution of the structure of the DNA-binding domain of interferon regulatory factor 4 bound to its interferon-stimulated response element by native intrinsic phosphorus and sulfur single-wavelength anomalous dispersion methods (native SAD) is described.

Pivotal to the regulation of key cellular processes such as the transcription, replication and repair of DNA, DNA-binding proteins play vital roles in all aspects of genetic activity. The determination of high-quality structures of DNA-binding proteins, particularly those in complexes with DNA, provides crucial insights into the understanding of these processes. The presence in such complexes of phosphate-rich oligonucleotides offers the choice of a rapid method for the routine solution of DNA-binding proteins through the use of long-wavelength beamlines such as I23 at Diamond Light Source. This article reports the use of native intrinsic phosphorus and sulfur single-wavelength anomalous dispersion methods to solve the complex of the DNA-binding domain (DBD) of interferon regulatory factor 4 (IRF4) bound to its interferon-stimulated response element (ISRE). The structure unexpectedly shows three molecules of the IRF4 DBD bound to one ISRE. The sole reliance on native intrinsic anomalous scattering elements that belong to DNA–protein complexes renders the method of general applicability to a large number of such protein complexes that cannot be solved by molecular replacement or by other phasing methods.

High Conformational Flexibility of the E2F1/DP1/DNA Complex

The E2F1 transcription factor is a master regulator of cell-cycle progression whose uncontrolled activation contributes to tumor cells growth. E2F1 binds DNA as a heterodimer with DP partners, resulting in a multi-domain quaternary-structure complex composed of DNA binding domains, a coiled coil domain and a marked box domain separated by short linkers. Building on the 3D knowledge of the single domains of E2F and DPs, we characterized the structure and dynamics of the complete E2F1/DP1/DNA complex by a combination of small-angle X-ray scattering and molecular dynamics simulations. It shows an asymmetric contribution of the dynamics of the two proteins. Namely, the coiled-coil domain leans toward the DP1 side of the complex; the DP1 loop between α2 and α3 of the DBD partially populates a helical structure leaning far from the DNA and in the same direction of the coiled-coil domain; and the N-terminal disordered region of DP1, rich in basic residues, contributes to DNA binding stabilization. Intriguingly, tumor mutations in the flexible regions of the complex suggest that perturbation of protein dynamics could affect protein function in a context-dependent way. Our data suggest fundamental contributions of DP proteins in distinct aspects of E2F biology.

The GRAS gene family in watermelons: identification, characterization and expression analysis of different tissues and root-knot nematode infestations

The family of GRAS plant-specific transcription factor plays diverse roles in numerous biological processes. Despite the identification and characterization of GRAS genes family in dozens of plant species, until now, GRAS members in watermelon (Citrullus lanatus) have not been investigated comprehensively. In this study, using bioinformatic analysis, we identified 37 GRAS genes in the watermelon genome (ClGRAS). These genes are classified into 10 distinct subfamilies based on previous research, and unevenly distributed on 11 chromosomes. Furthermore, a complete analysis was conducted to characterize conserved motifs and gene structures, which revealed the members within same subfamily that have analogous conserved gene structure and motif composition. Additionally, the expression pattern of ClGRAS genes was characterized in fruit flesh and rind tissues during watermelon fruit development and under red light (RL) as well as root knot nematode infestation. Finally, for verification of the availability of public transcriptome data, we also evaluated the expression levels of randomly selected four ClGRAS genes under RL and nematode infection by using qRT-PCR. The qRT-PCR results indicated that several ClGRAS genes were differentially expressed, implying their vital role in RL induction of watermelon resistance against root-knot nematodes. The results obtained in this study could be useful in improving the quality of watermelon.

A Tale of Loops and Tails: The Role of Intrinsically Disordered Protein Regions in R-Loop Recognition and Phase Separation

R-loops are non-canonical, three-stranded nucleic acid structures composed of a DNA:RNA hybrid, a displaced single-stranded (ss)DNA, and a trailing ssRNA overhang. R-loops perform critical biological functions under both normal and disease conditions. To elucidate their cellular functions, we need to understand the mechanisms underlying R-loop formation, recognition, signaling, and resolution. Previous high-throughput screens identified multiple proteins that bind R-loops, with many of these proteins containing folded nucleic acid processing and binding domains that prevent (e.g., topoisomerases), resolve (e.g., helicases, nucleases), or recognize (e.g., KH, RRMs) R-loops. However, a significant number of these R-loop interacting Enzyme and Reader proteins also contain long stretches of intrinsically disordered regions (IDRs). The precise molecular and structural mechanisms by which the folded domains and IDRs synergize to recognize and process R-loops or modulate R-loop-mediated signaling have not been fully explored. While studying one such modular R-loop Reader, the Fragile X Protein (FMRP), we unexpectedly discovered that the C-terminal IDR (C-IDR) of FMRP is the predominant R-loop binding site, with the three N-terminal KH domains recognizing the trailing ssRNA overhang. Interestingly, the C-IDR of FMRP has recently been shown to undergo spontaneous Liquid-Liquid Phase Separation (LLPS) assembly by itself or in complex with another non-canonical nucleic acid structure, RNA G-quadruplex. Furthermore, we have recently shown that FMRP can suppress persistent R-loops that form during transcription, a process that is also enhanced by LLPS via the assembly of membraneless transcription factories. These exciting findings prompted us to explore the role of IDRs in R-loop processing and signaling proteins through a comprehensive bioinformatics and computational biology study. Here, we evaluated IDR prevalence, sequence composition and LLPS propensity for the known R-loop interactome. We observed that, like FMRP, the majority of the R-loop interactome, especially Readers, contains long IDRs that are highly enriched in low complexity sequences with biased amino acid composition, suggesting that these IDRs could directly interact with R-loops, rather than being “mere flexible linkers” connecting the “functional folded enzyme or binding domains”. Furthermore, our analysis shows that several proteins in the R-loop interactome are either predicted to or have been experimentally demonstrated to undergo LLPS or are known to be associated with phase separated membraneless organelles. Thus, our overall results present a thought-provoking hypothesis that IDRs in the R-loop interactome can provide a functional link between R-loop recognition via direct binding and downstream signaling through the assembly of LLPS-mediated membrane-less R-loop foci. The absence or dysregulation of the function of IDR-enriched R-loop interactors can potentially lead to severe genomic defects, such as the widespread R-loop-mediated DNA double strand breaks that we recently observed in Fragile X patient-derived cells.

A comparative study of protein-ssDNA interactions

Single-stranded DNA-binding proteins (SSBs) play crucial roles in DNA replication, recombination and repair, and serve as key players in the maintenance of genomic stability. While a number of SSBs bind single-stranded DNA (ssDNA) non-specifically, the others recognize and bind specific ssDNA sequences. The mechanisms underlying this binding discrepancy, however, are largely unknown. Here, we present a comparative study of protein-ssDNA interactions by annotating specific and non-specific SSBs and comparing structural features such as DNA-binding propensities and secondary structure types of residues in SSB-ssDNA interactions, protein-ssDNA hydrogen bonding and π-π interactions between specific and non-specific SSBs. Our results suggest that protein side chain-DNA base hydrogen bonds are the major contributors to protein-ssDNA binding specificity, while π-π interactions may mainly contribute to binding affinity. We also found the enrichment of aspartate in the specific SSBs, a key feature in specific protein-double-stranded DNA (dsDNA) interactions as reported in our previous study. In addition, no significant differences between specific and non-specific groups with respect of conformational changes upon ssDNA binding were found, suggesting that the flexibility of SSBs plays a lesser role than that of dsDNA-binding proteins in conferring binding specificity.

High-Throughput Affinity Measurements of Transcription Factor and DNA Mutations Reveal Affinity and Specificity Determinants

Transcription factors (TFs) bind regulatory DNA to control gene expression, and mutations to either TFs or DNA can alter binding affinities to rewire regulatory networks and drive phenotypic variation. While studies have profiled energetic effects of DNA mutations extensively, we lack similar information for TF variants. Here, we present STAMMP (simultaneous transcription factor affinity measurements via microfluidic protein arrays), a high-throughput microfluidic platform enabling quantitative characterization of hundreds of TF variants simultaneously. Measured affinities for ∼210 mutants of a model yeast TF (Pho4) interacting with 9 oligonucleotides (>1,800 K_ds) reveal that many combinations of mutations to poorly conserved TF residues and nucleotides flanking the core binding site alter but preserve physiological binding, providing a mechanism by which combinations of mutations in cis and trans could modulate TF binding to tune occupancies during evolution. Moreover, biochemical double-mutant cycles across the TF-DNA interface reveal molecular mechanisms driving recognition, linking sequence to function. A record of this paper's Transparent Peer Review process is included in the Supplemental Information.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

Capturing the Conformational Ensemble of the Mixed Folded Polyglutamine Protein Ataxin-3

Ataxin-3 is a deubiquitinase involved in protein quality control and other essential cellular functions. It preferentially interacts with polyubiquitin chains of four or more units attached to proteins delivered to the ubiquitin-proteasome system. Ataxin-3 is composed of an N-terminal Josephin domain and a flexible C terminus that contains two or three ubiquitin-interacting motifs (UIMs) and a polyglutamine tract, which, when expanded beyond a threshold, leads to protein aggregation and misfolding and causes spinocerebellar ataxia type 3. The high-resolution structure of the Josephin domain is available, but the structural and dynamical heterogeneity of ataxin-3 has so far hindered the structural description of the full-length protein. Here, we characterize non-expanded and expanded variants of ataxin-3 in terms of conformational ensembles adopted by the proteins in solution by jointly using experimental data from nuclear magnetic resonance and small-angle X-ray scattering with coarse-grained simulations. Our results pave the way to a molecular understanding of polyubiquitin recognition.

Aggregation and disaggregation features of the human proteome

Protein aggregates have negative implications in disease. While reductionist experiments have increased our understanding of aggregation processes, the systemic view in biological context is still limited. To extend this understanding, we used mass spectrometry-based proteomics to characterize aggregation and disaggregation in human cells after non-lethal heat shock. Aggregation-prone proteins were enriched in nuclear proteins, high proportion of intrinsically disordered regions, high molecular mass, high isoelectric point, and hydrophilic amino acids. During recovery, most aggregating proteins disaggregated with a rate proportional to the aggregation propensity: larger loss in solubility was counteracted by faster disaggregation. High amount of intrinsically disordered regions were associated with faster disaggregation. However, other characteristics enriched in aggregating proteins did not correlate with the disaggregation rates. In addition, we analyzed changes in protein thermal stability after heat shock. Soluble remnants of aggregated proteins were more thermally stable compared with control condition. Therefore, our results provide a rich resource of heat stress-related protein solubility data and can foster further studies related to protein aggregation diseases.

C1188D mutation abolishes specific recognition between MLL1-CXXC domain and CpG site by inducing conformational switch of flexible N-terminal

Mixed lineage leukemia protein (MLL1 protein) recognizes the CpG site via its CXXC domain and is frequently associated with leukemia. The specific recognition is abolished by C1188D mutation, which also prevents MLL-related leukemia. In this paper, multiple molecular dynamic (MD) simulations were performed to investigate the mechanism of recognition and influences of C1188D mutation. Started from fully dissociated DNA and MLL1-CXXC domain, remarkably, the center of mass (COM) of MLL1-CXXC domain quickly concentrates on the vicinity of the CpG site in all 53 short MD simulations. Extended simulations of the wild type showed that the native complex formed in 500 ns among 4 of 53 simulations. In contrast, the C1188D mutant COM distributed broadly around the DNA and the native complex was not observed in any of the extended simulations. Simulations on the apo MLL1-CXXC domain further suggest that the wild type protein remained predominantly in an open form that closely resembles its structure in the native complex whereas C1188D mutant formed predominantly compact structures in which the N- terminal bends to D1188. This conformational switch hinders the formation of encounter complex, thus abolishes the recognition. Our study also provides clues to the study mechanism of recognition, by the CXXC domain from proteins like DNA methyltransferase and ten-eleven translocation enzymes.

New insights into protein-DNA binding specificity from hydrogen bond based comparative study

Knowledge of protein-DNA binding specificity has important implications in understanding DNA metabolism, transcriptional regulation and developing therapeutic drugs. Previous studies demonstrated hydrogen bonds between amino acid side chains and DNA bases play major roles in specific protein-DNA interactions. In this paper, we investigated the roles of individual DNA strands and protein secondary structure types in specific protein-DNA recognition based on side chain-base hydrogen bonds. By comparing the contribution of each DNA strand to the overall binding specificity between DNA-binding proteins with different degrees of binding specificity, we found that highly specific DNA-binding proteins show balanced hydrogen bonding with each of the two DNA strands while multi-specific DNA binding proteins are generally biased towards one strand. Protein-base pair hydrogen bonds, in which both bases of a base pair are involved in forming hydrogen bonds with amino acid side chains, are more prevalent in the highly specific protein-DNA complexes than those in the multi-specific group. Amino acids involved in side chain-base hydrogen bonds favor strand and coil secondary structure types in highly specific DNA-binding proteins while multi-specific DNA-binding proteins prefer helices.

Genome-wide identification, structural analysis and expression profiles of GRAS gene family in orchardgrass

The GRAS gene family is a family of transcription factors that regulates plant growth and development. Despite being well-studied in many plant species, little is known about this gene family in orchardgrass (Dactylis glomerata L.), one of the top four economically important perennial forage grasses cultivated worldwide. We identified 46 GRAS genes in orchardgrass and analyzed their characteristics by phylogenetic, gene structural, motifs and expression patterns analysis. The phylogenetic analysis of eight species revealed that DgGRAS family had the evolutional conservation and closer homology relationship with the GRAS family of rice, barley and Brachypodium distachyon. Moreover, 46 DgGRAS proteins were divided into eight subfamilies based on the tree topology and rice or Arabidopsis classification, and LISCL subfamily was the largest one. Besides, we found that the motif 15 may be unique to the orchardgrass LISCL subfamily, and the motif 6 and motif 17 had indispensable functions in the orchardgrass LISCL subfamily. We further analyzed the expression profiles of DgGRAS genes at mature and seeding stage. And we found that DgGRAS17 played an important role in the growth and development no matter what stage it was at. DgGRAS5, DgGRAS28, DgGRAS31, DgGRAS42 and DgGRAS44 got involved in processes of the growth and development at seeding stage instead of mature stage. These results indicated that the major expression patterns and detailed functions of the DgGRAS genes varied with developmental stages. Taken together, this is the first systematic analysis of the GRAS gene family in the orchardgrass genome and the results provide insights into the potential functions of GRAS genes.

The functional diversity of structural disorder in plant proteins

Structural disorder in proteins is a widespread feature distributed in all domains of life, particularly abundant in eukaryotes, including plants. In these organisms, intrinsically disordered proteins (IDPs) perform a diversity of functions, participating as integrators of signaling networks, in transcriptional and post-transcriptional regulation, in metabolic control, in stress responses and in the formation of biomolecular condensates by liquid-liquid phase separation. Their roles impact the perception, propagation and control of various developmental and environmental cues, as well as the plant defense against abiotic and biotic adverse conditions. In this review, we focus on primary processes to exhibit a broad perspective of the relevance of IDPs in plant cell functions. The information here might help to incorporate this knowledge into a more dynamic view of plant cells, as well as open more questions and promote new ideas for a better understanding of plant life.

Structural Fuzziness of the RNA-Organizing Protein SERF Determines a Toxic Gain-of-interaction

The mechanisms by which protein complexes convert from functional to pathogenic are the subject of intensive research. Here, we report how functionally unfavorable protein interactions can be induced by structural fuzziness, i.e., by persisting conformational disorder in protein complexes. We show that extreme disorder in the bound state transforms the intrinsically disordered protein SERF1a from an RNA-organizing factor into a pathogenic enhancer of alpha-synuclein (aSyn) amyloid toxicity. We demonstrate that SERF1a promotes the incorporation of RNA into nucleoli and liquid-like artificial RNA-organelles by retaining an unusually high degree of conformational disorder in the RNA-bound state. However, this type of structural fuzziness also determines an undifferentiated interaction with aSyn. RNA and aSyn both bind to one identical, positively charged site of SERF1a by an analogous electrostatic binding mode, with similar binding affinities, and without any observable disorder-to-order transition. The absence of primary or secondary structure discriminants results in SERF1a being unable to select between nucleic acid and amyloidogenic protein, leading the pro-amyloid aSyn:SERF1a interaction to prevail in the cytosol under conditions of cellular stress. We suggest that fuzzy disorder in SERF1a complexes accounts for an adverse gain-of-interaction which favors toxic binding to aSyn at the expense of nontoxic RNA binding, thereby leading to a functionally distorted and pathogenic process. Thus, structural fuzziness constitutes a direct link between extreme conformational flexibility, amyloid aggregation, and the malfunctioning of RNA-associated cellular processes, three signatures of neurodegenerative proteinopathies.

Interactions of p53 with poly(ADP-ribose) and DNA induce distinct changes in protein structure as revealed by ATR-FTIR spectroscopy

Due to multiple domains and in part intrinsically disordered regions, structural analyses of p53 remain a challenging task, particularly in complex with DNA and other macromolecules. Here, we applied a novel attenuated total reflection Fourier transform infrared (ATR-FTIR) spectroscopic approach to investigate changes in secondary structure of full-length p53 induced by non-covalent interactions with DNA and poly(ADP-ribose) (PAR). To validate our approach, we confirmed a positive regulatory function of p53’s C-terminal domain (CTD) with regard to sequence-specific DNA binding and verified that the CTD mediates p53–PAR interaction. Further, we demonstrate that DNA and PAR interactions result in distinct structural changes of p53, indicating specific binding mechanisms via different domains. A time-dependent analysis of the interplay of DNA and PAR binding to p53 revealed that PAR represents p53’s preferred binding partner, which efficiently controls p53–DNA interaction. Moreover, we provide infrared spectroscopic data on PAR pointing to the absence of regular secondary structural elements. Finally, temperature-induced melting experiments via CD spectroscopy show that DNA binding stabilizes the structure of p53, while PAR binding can shift the irreversible formation of insoluble p53 aggregates to higher temperatures. In conclusion, this study provides detailed insights into the dynamic interplay of p53 binding to DNA and PAR at a formerly inaccessible molecular level.

The Intrinsically Disordered Loop in the USF1 bHLHZ Domain Modulates Its DNA-Binding Sequence Specificity in Hereditary Asthma

USF1, a basic region/helix-loop-helix/leucine zipper (bHLHZ) transcription factor, binds to the E-box in the PAI-1 (plasminogen activator inhibitor) promoter. Two alleles containing the E-box control PAI-1 transcription; these alleles are termed "4G" and "5G" based on the G tract flanking E-box. USF1-governed transcription of PAI-1 is elevated in heritable asthma sufferers: the 4G/4G genotype has the highest plasma levels of PAI-1. While USF1 uses its basic region to bind E-box, we found that it uses its 12 amino-acid loop to recognize the flanking sequence and discern the single-nucleotide difference between the alleles. We used the bacterial one-hybrid and electrophoretic mobility shift assays to assess protein-DNA recognition, and circular dichroism to examine protein secondary structure. We mutated Ser233 and Thr234 in the USF1 bHLHZ loop to Ala to generate S233A and T234A. Interestingly, USF1 bHLHZ, S233A, and T234A prefer the 5G sequence (USF1 bHLHZ K_d values 4.1 ± 0.3 nM and 7.0 ± 0.4 nM for 5G and 4G, respectively), whereas studies in stimulated human mast cells showed a preference for 4G. We replaced the 8 amino-acid loop of transcription factor Max bHLHZ with the 12 amino-acid USF1 loop: this mutant now distinguishes the 4G/5G polymorphism-while Max bHLHZ does not-confirming that USF1 differentiation of the 4G/5G is driven by the loop.

Intrinsically disordered domains: Sequence ➔ disorder ➔ function relationships

Disordered domains are long regions of intrinsic disorder that ideally have conserved sequences, conserved disorder, and conserved functions. These domains were first noticed in protein-protein interactions that are distinct from the interactions between two structured domains and the interactions between structured domains and linear motifs or molecular recognition features (MoRFs). So far, disordered domains have not been systematically characterized. Here, we present a bioinformatics investigation of the sequence-disorder-function relationships for a set of probable disordered domains (PDDs) identified from the Pfam database. All the Pfam seed proteins from those domains with at least one PDD sequence were collected. Most often, if a set contains one PDD sequence, then all members of the set are PDDs or nearly so. However, many seed sets have sequence collections that exhibit diverse proportions of predicted disorder and structure, thus giving the completely unexpected result that conserved sequences can vary substantially in predicted disorder and structure. In addition to the induction of structure by binding to protein partners, disordered domains are also induced to form structure by disulfide bond formation, by ion binding, and by complex formation with RNA or DNA. The two new findings, (a) that conserved sequences can vary substantially in their predicted disorder content and (b) that homologues from a single domain can evolve from structure to disorder (or vice versa), enrich our understanding of the sequence ➔ disorder ensemble ➔ function paradigm.

bHLH-PAS Proteins: Their Structure and Intrinsic Disorder

The basic helix–loop–helix/Per-ARNT-SIM (bHLH–PAS) proteins are a class of transcriptional regulators, commonly occurring in living organisms and highly conserved among vertebrates and invertebrates. These proteins exhibit a relatively well-conserved domain structure: the bHLH domain located at the N-terminus, followed by PAS-A and PAS-B domains. In contrast, their C-terminal fragments present significant variability in their primary structure and are unique for individual proteins. C-termini were shown to be responsible for the specific modulation of protein action. In this review, we present the current state of knowledge, based on NMR and X-ray analysis, concerning the structural properties of bHLH–PAS proteins. It is worth noting that all determined structures comprise only selected domains (bHLH and/or PAS). At the same time, substantial parts of proteins, comprising their long C-termini, have not been structurally characterized to date. Interestingly, these regions appear to be intrinsically disordered (IDRs) and are still a challenge to research. We aim to emphasize the significance of IDRs for the flexibility and function of bHLH–PAS proteins. Finally, we propose modern NMR methods for the structural characterization of the IDRs of bHLH–PAS proteins.

Can All-Atom Molecular Dynamics Simulations Quantitatively Describe Homeodomain-DNA Binding Equilibria?

We systematically investigate the applicability of a molecular dynamics-based setup for the calculations of standard binding free energies of biologically relevant protein-DNA complexes. The free energies are extracted from a potential of mean force calculated using umbrella sampling simulations. Two protein-DNA systems derived from a homeodomain transcription factor complex are studied in order to investigate the binding of both disordered and globular proteins. Free energies and trajectories obtained using two modern molecular mechanical force fields are compared to each other and to experimental data. The temperature dependence of the calculated standard binding free energies is investigated by performing all simulations over a range of temperatures. We show that the values of standard binding free energies obtained from these simulations are overestimated compared to experimental results. Significant differences are observed between the two protein-DNA systems and between the two force fields, which are explained by different propensities to form inter- and intramolecular contacts. The number of protein-DNA contacts increases with increasing temperature, in agreement with the experimentally known temperature dependence of enthalpies of binding. However, conclusions about the temperature dependence of the standard binding free energies cannot be made with confidence, as the differences among the values are on the order of statistical uncertainty.

Position-, disorder-, and salt-dependent diffusion in binding-coupled-folding of intrinsically disordered proteins

Successful extensions of protein-folding energy landscape theory to intrinsically disordered proteins' (IDPs') binding-coupled-folding transition can enormously simplify this biomolecular process into diffusion along a limited number of reaction coordinates, and the dynamics subsequently is described by Kramers' rate theory. As the critical pre-factor, the diffusion coefficient D has direct implications on the binding kinetics. Here, we employ a structure-based model (SBM) to calculate D in the binding-folding of an IDP prototype. We identify a strong position-dependent D during binding by applying a reaction coordinate that directly measures the fluctuations in a Cartesian configuration space. Using the malleability of the SBM, we modulate the degree of conformational disorder in an isolated IDP and determine complex effects of intrinsic disorder on D varying for different binding stages. Here, D tends to increase with disorder during initial binding but shows a non-monotonic relationship with disorder in terms of a decrease-followed-by-increase in D during the late binding stage. The salt concentration, which correlates with electrostatic interactions via Debye-Hückel theory in our SBM, also modulates D in a stepwise way. The speeding up of diffusion by electrostatic interactions is observed during the formation of the encounter complex at the beginning of binding, while the last diffusive binding dynamics is hindered by non-native salt bridges. Because D describes the diffusive speed locally, which implicitly reflects the roughness of the energy landscape, we are eventually able to portray the binding energy landscape, including that from IDPs' binding, then to binding with partial folding, and finally to rigid docking, as well as that under different environmental salt concentrations. Our theoretical results provide key mechanistic insights into IDPs' binding-folding, which is internally conformation- and externally salt-controlled with respect to diffusion.

Selection of DNA Cleavage Sites by Topoisomerase II Results from Enzyme-Induced Flexibility of DNA

Topoisomerase II cleaves DNA at preferred sequences with different efficiencies; however, the mechanism of cleavage site selection is not known. Here we used single-molecule fluorescence assays that monitor several critical steps of DNA-topoisomerase II interactions, including binding/dissociation, bending/straightening, and cleavage/religation, and reveal that the cleavage site is selected mainly during the bending step. Furthermore, despite the sensitivity of the bending rate to the DNA sequence, it is not an intrinsic property of the DNA itself. Rather, it is determined by protein-DNA interactions.

Towards a Stochastic Paradigm: From Fuzzy Ensembles to Cellular Functions

The deterministic sequence → structure → function relationship is not applicable to describe how proteins dynamically adapt to different cellular conditions. A stochastic model is required to capture functional promiscuity, redundant sequence motifs, dynamic interactions, or conformational heterogeneity, which facilitate the decision-making in regulatory processes, ranging from enzymes to membraneless cellular compartments. The fuzzy set theory offers a quantitative framework to address these problems. The fuzzy formalism allows the simultaneous involvement of proteins in multiple activities, the degree of which is given by the corresponding memberships. Adaptation is described via a fuzzy inference system, which relates heterogeneous conformational ensembles to different biological activities. Sequence redundancies (e.g., tandem motifs) can also be treated by fuzzy sets to characterize structural transitions affecting the heterogeneous interaction patterns (e.g., pathological fibrillization of stress granules). The proposed framework can provide quantitative protein models, under stochastic cellular conditions.

Plasticity in oligomerization, operator architecture, and DNA binding in the mode of action of a bacterial B12-based photoreceptor

Newly discovered bacterial photoreceptors called CarH sense light by using 5'-deoxyadenosylcobalamin (AdoCbl). They repress their own expression and that of genes for carotenoid synthesis by binding in the dark to operator DNA as AdoCbl-bound tetramers, whose light-induced disassembly relieves repression. High-resolution structures of Thermus thermophilus CarH_Tt have provided snapshots of the dark and light states and have revealed a unique DNA-binding mode whereby only three of four DNA-binding domains contact an operator comprising three tandem direct repeats. To gain further insights into CarH photoreceptors and employing biochemical, spectroscopic, mutational, and computational analyses, here we investigated CarH_Bm from Bacillus megaterium We found that apoCarH_Bm, unlike monomeric apoCarH_Tt, is an oligomeric molten globule that forms DNA-binding tetramers in the dark only upon AdoCbl binding, which requires a conserved W-X ₉-EH motif. Light relieved DNA binding by disrupting CarH_Bm tetramers to dimers, rather than to monomers as with CarH_Tt CarH_Bm operators resembled that of CarH_Tt, but were larger by one repeat and overlapped with the -35 or -10 promoter elements. This design persisted in a six-repeat, multipartite operator we discovered upstream of a gene encoding an Spx global redox-response regulator whose photoregulated expression links photooxidative and general redox responses in B. megaterium Interestingly, CarH_Bm recognized the smaller CarH_Tt operator, revealing an adaptability possibly related to the linker bridging the DNA- and AdoCbl-binding domains. Our findings highlight a remarkable plasticity in the mode of action of B₁₂-based CarH photoreceptors, important for their biological functions and development as optogenetic tools.

How a single 5-methylation of cytosine regulates the recognition of C/EBPβ transcription factor: a molecular dynamic simulation study

CpG methylation can regulate gene expression by altering the specific binding of protein and DNA. In order to understand how a single 5mC regulates protein-DNA interactions, we have compared the structures and dynamics of CEBP/βprotein-DNA complexes before and after methylation, and the results indicate that even a single 5mC can regulate protein-DNA recognition by steric-hindrance effect of methyl group and changing the hydrogen bond interactions. The interactions between the methyl group, mCpG motif, and the conserved residue arginine make the protein read out the variation of local environment, which further enhances the specific recognition and affects the base pair stacking. The stacking interactions can propagate along the backbone of DNA and lead to long-range allosteric effects, including obvious conformational variations for DNA base pairs.

A low-complexity region in the YTH domain protein Mmi1 enhances RNA binding

Mmi1 is an essential RNA-binding protein in the fission yeast Schizosaccharomyces pombe that eliminates meiotic transcripts during normal vegetative growth. Mmi1 contains a YTH domain that binds specific RNA sequences, targeting mRNAs for degradation. The YTH domain of Mmi1 uses a noncanonical RNA-binding surface that includes contacts outside the conserved fold. Here, we report that an N-terminal extension that is proximal to the YTH domain enhances RNA binding. Using X-ray crystallography, NMR, and biophysical methods, we show that this low-complexity region becomes more ordered upon RNA binding. This enhances the affinity of the interaction of the Mmi1 YTH domain with specific RNAs by reducing the dissociation rate of the Mmi1-RNA complex. We propose that the low-complexity region influences RNA binding indirectly by reducing dynamic motions of the RNA-binding groove and stabilizing a conformation of the YTH domain that binds to RNA with high affinity. Taken together, our work reveals how a low-complexity region proximal to a conserved folded domain can adopt an ordered structure to aid nucleic acid binding.

Fuzziness in Protein Interactions-A Historical Perspective

The proposal that coupled folding to binding is not an obligatory mechanism for intrinsically disordered (ID) proteins was put forward 10 years ago. The notion of fuzziness implies that conformational heterogeneity can be maintained upon interactions of ID proteins, which has a functional impact either on regulated assembly or activity of the corresponding complexes. Here I review how the concept has evolved in the past decade, via increasing experimental data providing insights into the mechanisms, pathways and regulatory modes. The effects of structural diversity and transient contacts on protein assemblies have been collected and systematically analyzed (Fuzzy Complexes Database, http://protdyn-database.org). Fuzziness has also been exploited as a framework to decipher molecular organization of higher-order protein structures. Quantification of conformational heterogeneity opens exciting future perspectives for drug discovery from small molecule-ID protein interactions to supramolecular assemblies.

Roles of conformational disorder and downhill folding in modulating protein-DNA recognition

Transcription factors are thought to efficiently search for their target DNA site via a combination of conventional 3D diffusion and 1D diffusion along the DNA molecule mediated by non-specific electrostatic interactions. This process requires the DNA-binding protein to quickly exchange between a search competent and a target recognition mode, but little is known as to how these two binding modes are encoded in the conformational properties of the protein. Here, we investigate this issue on the engrailed homeodomain (EngHD), a DNA-binding domain that folds ultrafast and exhibits a complex conformational behavior consistent with the downhill folding scenario. We explore the interplay between folding and DNA recognition using a coarse-grained computational model that allows us to manipulate the folding properties of the protein and monitor its non-specific and specific binding to DNA. We find that conformational disorder increases the search efficiency of EngHD by promoting a fast gliding search mode in addition to sliding. When gliding, EngHD remains loosely bound to DNA moving linearly along its length. A partially disordered EngHD also binds more dynamically to the target site, reducing the half-life of the specific complex via a spring-loaded mechanism. These findings apply to all conditions leading to partial disorder. However, we also find that at physiologically relevant temperatures EngHD is well folded and can only obtain the conformational flexibility required to accelerate 1D diffusion when it folds/unfolds within the downhill scenario (crossing a marginal free energy barrier). In addition, the conformational flexibility of native downhill EngHD enables its fast reconfiguration to lock into the specific binding site upon arrival, thereby affording finer control of the on- and off-rates of the specific complex. Our results provide key mechanistic insights into how DNA-binding domains optimize specific DNA recognition through the control of their conformational dynamics and folding mechanism.

Determination of helix orientations in a flexible DNA by multi-frequency EPR spectroscopy

Distance measurements are performed between a pair of spin labels attached to nucleic acids using Pulsed Electron-Electron Double Resonance (PELDOR, also called DEER) spectroscopy which is a complementary tool to other structure determination methods in structural biology. The rigid spin label Ç, when incorporated pairwise into two helical parts of a nucleic acid molecule, allows the determination of both the mutual orientation and the distance between those labels, since Ç moves rigidly with the helix to which it is attached. We have developed a two-step protocol to investigate the conformational flexibility of flexible nucleic acid molecules by multi-frequency PELDOR. In the first step, a library with a broad collection of conformers, which are in agreement with topological constraints, NMR restraints and distances derived from PELDOR, was created. In the second step, a weighted structural ensemble of these conformers was chosen, such that it fits the multi-frequency PELDOR time traces of all doubly Ç-labelled samples simultaneously. This ensemble reflects the global structure and the conformational flexibility of the two-way DNA junction. We demonstrate this approach on a flexible bent DNA molecule, consisting of two short helical parts with a five adenine bulge at the center. The kink and twist motions between both helical parts were quantitatively determined and showed high flexibility, in agreement with a Förster Resonance Energy Transfer (FRET) study on a similar bent DNA motif. The approach presented here should be useful to describe the relative orientation of helical motifs and the conformational flexibility of nucleic acid structures, both alone and in complexes with proteins and other molecules.

Oligomer formation and G-quadruplex binding by purified murine Rif1 protein, a key organizer of higher-order chromatin architecture

Rap1-interacting protein 1 (Rif1) regulates telomere length in budding yeast. We previously reported that, in metazoans and fission yeast, Rif1 also plays pivotal roles in controlling genome-wide DNA replication timing. We proposed that Rif1 may assemble chromatin compartments that contain specific replication-timing domains by promoting chromatin loop formation. Rif1 also is involved in DNA lesion repair, restart after replication fork collapse, anti-apoptosis activities, replicative senescence, and transcriptional regulation. Although multiple physiological functions of Rif1 have been characterized, biochemical and structural information on mammalian Rif1 is limited, mainly because of difficulties in purifying the full-length protein. Here, we expressed and purified the 2418-amino-acid-long, full-length murine Rif1 as well as its partially truncated variants in human 293T cells. Hydrodynamic analyses indicated that Rif1 forms elongated or extended homo-oligomers in solution, consistent with the presence of a HEAT-type helical repeat segment known to adopt an elongated shape. We also observed that the purified murine Rif1 bound G-quadruplex (G4) DNA with high specificity and affinity, as was previously shown for Rif1 from fission yeast. Both the N-terminal (HEAT-repeat) and C-terminal segments were involved in oligomer formation and specifically bound G4 DNA, and the central intrinsically disordered polypeptide segment increased the affinity for G4. Of note, pulldown assays revealed that Rif1 simultaneously binds multiple G4 molecules. Our findings support a model in which Rif1 modulates chromatin loop structures through binding to multiple G4 assemblies and by holding chromatin fibers together.

Conformational heterogeneity in tails of DNA-binding proteins is augmented by proline containing repeats

A cationic terminal extension or tail is a common feature of many DNA-binding proteins. We show that a particular type of tail rich in proline, alanine and lysine belongs to the class of 'flexible disorder' and consists of characteristic pentapeptide repeats. Our designed peptides, (AAKKA)_1-4 and (PAKKA)_1-4, represent the tails of several bacterial DNA-binding proteins. Enhanced conformational sampling of these representative peptides using accelerated molecular dynamic simulations supported by circular dichroism spectroscopy and nuclear magnetic resonance studies demonstrates the role of frequent and interspersed prolines in augmenting conformational heterogeneity of the peptide backbone. Analysis of circular variance of backbone dihedral angles indicates alternating regions of relative rigidity and flexibility along the peptide sequence due to prolines. Preferred placement of lysines in the regions of higher backbone flexibility might improve DNA-binding by conformational selection. Our results could be relevant for rational de novo design of disordered peptides.

Evidence for a Strong Correlation Between Transcription Factor Protein Disorder and Organismic Complexity

Studies of diverse phylogenetic lineages reveal that protein disorder increases in concert with organismic complexity but that differences nevertheless exist among lineages. To gain insight into this phenomenology, we analyzed all of the transcription factor (TF) families for which sequences are known for 17 species spanning bacteria, yeast, algae, land plants, and animals and for which the number of different cell types has been reported in the primary literature. Although the fraction of disordered residues in TF sequences is often moderately or poorly correlated with organismic complexity as gauged by cell-type number (r² < 0.5), an unbiased and phylogenetically broad analysis shows that organismic complexity is positively and strongly correlated with the total number of TFs, the number of their spliced variants and their total disordered residues content (r² > 0.8). Furthermore, the correlation between the fraction of disordered residues and cell-type number becomes stronger when confined to the TF families participating in cell cycle, cell size, cell division, cell differentiation, or cell proliferation, and other important developmental processes. The data also indicate that evolutionarily simpler organisms allow for the detection of subtle differences in the conserved IDRs of TFs as well as changes in variable IDRs, which can influence the DNA recognition and multifunctionality of TFs through direct or indirect mechanisms. Although strong correlations cannot be taken as evidence for cause-and-effect relationships, we interpret our data to indicate that increasing TF disorder likely was an important factor contributing to the evolution of organismic complexity and not merely a concurrent unrelated effect of increasing organismic complexity.

Identification and expression of GRAS family genes in maize (Zea mays L.)

GRAS transcriptional factors have diverse functions in plant growth and development, and are named after the first three transcription factors, namely, GAI (GIBBERELLIC ACID INSENSITIVE), RGA (REPRESSOR OF GAI) and SCR (SCARECROW) identified in this family. Knowledge of the GRAS gene family in maize remains was largely unknown, and their characterization is necessary to understand their importance in the maize life cycle. This study identified 86 GRAS genes in maize, and further characterized with phylogenetics, gene structural analysis, genomic loci, and expression patterns. The 86 GRAS genes were divided into 8 groups (SCL3, HAM, LS, SCR, DELLA, SHR, PAT1 and LISCL) by phylogenetic analysis. Most of the maize GRAS genes contain one exon (80.23%) and closely related members in the phylogenetic tree had similar structure and motif composition. Different motifs especially in the N-terminus might be the sources of their functional divergence. Segmental- and tandem-duplication occurred in this family leading to expansion of maize GRAS genes and the expression patterns of the duplicated genes in the heat map according to the published microarray data were very similar. Quantitative RT-PCR (qRT-PCR) results demonstrated that the expression level of genes in different tissues were different, suggesting their differential roles in plant growth and development. The data set expands our knowledge to understanding the function of GRAS genes in maize, an important crop plant in the world.

A Comprehensive Survey of the Roles of Highly Disordered Proteins in Type 2 Diabetes

Type 2 diabetes mellitus (T2DM) is a chronic and progressive disease that is strongly associated with hyperglycemia (high blood sugar) related to either insulin resistance or insufficient insulin production. Among the various molecular events and players implicated in the manifestation and development of diabetes mellitus, proteins play several important roles. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database has information on 34 human proteins experimentally shown to be related to the T2DM pathogenesis. It is known that many proteins associated with different human maladies are intrinsically disordered as a whole, or contain intrinsically disordered regions. The presented study shows that T2DM is not an exception to this rule, and many proteins known to be associated with pathogenesis of this malady are intrinsically disordered. The multiparametric bioinformatics analysis utilizing several computational tools for the intrinsic disorder characterization revealed that IRS1, IRS2, IRS4, MAFA, PDX1, ADIPO, PIK3R2, PIK3R5, SoCS1, and SoCS3 are expected to be highly disordered, whereas VDCC, SoCS2, SoCS4, JNK9, PRKCZ, PRKCE, insulin, GCK, JNK8, JNK10, PYK, INSR, TNF-α, MAPK3, and Kir6.2 are classified as moderately disordered proteins, and GLUT2, GLUT4, mTOR, SUR1, MAPK1, IKKA, PRKCD, PIK3CB, and PIK3CA are predicted as mostly ordered. More focused computational analyses and intensive literature mining were conducted for a set of highly disordered proteins related to T2DM. The resulting work represents a comprehensive survey describing the major biological functions of these proteins and functional roles of their intrinsically disordered regions, which are frequently engaged in protein–protein interactions, and contain sites of various posttranslational modifications (PTMs). It is also shown that intrinsic disorder-associated PTMs may play important roles in controlling the functions of these proteins. Consideration of the T2DM proteins from the perspective of intrinsic disorder provides useful information that can potentially lead to future experimental studies that may uncover latent and novel pathways associated with the disease.