Search

Abstract Background Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Acute myeloid leukemia (AML) and myeloid neoplasms develop through acquisition of somatic mutations that confer mutation-specific fitness advantages to hematopoietic stem and progenitor cells. However, our understanding of mutational effects remains limited to the resolution attainable within immunophenotypically and clinically accessible bulk cell populations. To decipher heterogeneous cellular fitness to preleukemic mutational perturbations, we performed single-cell RNA sequencing of eight different mouse models with driver mutations of myeloid malignancies, generating 269,048 single-cell profiles. Our analysis infers mutation-driven perturbations in cell abundance, cellular lineage fate, cellular metabolism, and gene expression at the continuous resolution, pinpointing cell populations with transcriptional alterations associated with differentiation bias. We further develop an 11-gene scoring system (Stem11) on the basis of preleukemic transcriptional signatures that predicts AML patient outcomes. Our results demonstrate that a single-cell-resolution deep characterization of preleukemic biology has the potential to enhance our understanding of AML heterogeneity and inform more effective risk stratification strategies., Competing Interests: Aspects of this work are included in United Kingdom patent application 2312684.0., (© 2023 The Author(s).)

The extraembryonic yolk sac (YS) ensures delivery of nutritional support and oxygen to the developing embryo but remains ill-defined in humans. We therefore assembled a comprehensive multiomic reference of the human YS from 3 to 8 postconception weeks by integrating single-cell protein and gene expression data. Beyond its recognized role as a site of hematopoiesis, we highlight roles in metabolism, coagulation, vascular development, and hematopoietic regulation. We reconstructed the emergence and decline of YS hematopoietic stem and progenitor cells from hemogenic endothelium and revealed a YS-specific accelerated route to macrophage production that seeds developing organs. The multiorgan functions of the YS are superseded as intraembryonic organs develop, effecting a multifaceted relay of vital functions as pregnancy proceeds.

Mast cell accumulation is a hallmark of a number of diseases, including allergic asthma and systemic mastocytosis. Immunoglobulin E-mediated crosslinking of the FcεRI receptors causes mast cell activation and contributes to disease pathogenesis. The mast cell lineage is one of the least studied among the hematopoietic cell lineages, and controversies remain about whether FcεRI expression appears during the mast cell progenitor stage or during terminal mast cell maturation. Here, we used single-cell transcriptomics analysis to reveal a temporal association between the appearance of FcεRI and the mast cell gene signature in CD34+ hematopoietic progenitors in adult peripheral blood. In agreement with these data, the FcεRI+ hematopoietic progenitors formed morphologically, phenotypically, and functionally mature mast cells in long-term culture assays. Single-cell transcriptomics analysis further revealed the expression patterns of prospective cytokine receptors regulating development of mast cell progenitors. Culture assays showed that interleukin-3 (IL-3) and IL-5 promoted disparate effects on progenitor cell proliferation and survival, respectively, whereas IL-33 caused robust FcεRI downregulation. Taken together, we showed that FcεRI expression appears at the progenitor stage of mast cell differentiation in peripheral blood. We also showed that external stimuli regulate FcεRI expression of mast cell progenitors, providing a possible explanation for the variable FcεRI expression levels during mast cell development., (© 2022 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

Clonal hematopoiesis (CH), the clonal expansion of a blood stem cell and its progeny driven by somatic driver mutations, affects over a third of people, yet remains poorly understood. Here we analyze genetic data from 200,453 UK Biobank participants to map the landscape of inherited predisposition to CH, increasing the number of germline associations with CH in European-ancestry populations from 4 to 14. Genes at new loci implicate DNA damage repair (PARP1, ATM, CHEK2), hematopoietic stem cell migration/homing (CD164) and myeloid oncogenesis (SETBP1). Several associations were CH-subtype-specific including variants at TCL1A and CD164 that had opposite associations with DNMT3A- versus TET2-mutant CH, the two most common CH subtypes, proposing key roles for these two loci in CH development. Mendelian randomization analyses showed that smoking and longer leukocyte telomere length are causal risk factors for CH and that genetic predisposition to CH increases risks of myeloproliferative neoplasia, nonhematological malignancies, atrial fibrillation and blood epigenetic ageing., (© 2022. The Author(s).)

Clonal expansions driven by somatic mutations become pervasive across human tissues with age, including in the haematopoietic system, where the phenomenon is termed clonal haematopoiesis 1-4 . The understanding of how and when clonal haematopoiesis develops, the factors that govern its behaviour, how it interacts with ageing and how these variables relate to malignant progression remains limited 5,6 . Here we track 697 clonal haematopoiesis clones from 385 individuals 55 years of age or older over a median of 13 years. We find that 92.4% of clones expanded at a stable exponential rate over the study period, with different mutations driving substantially different growth rates, ranging from 5% (DNMT3A and TP53) to more than 50% per year (SRSF2 P95H ). Growth rates of clones with the same mutation differed by approximately ±5% per year, proportionately affecting slow drivers more substantially. By combining our time-series data with phylogenetic analysis of 1,731 whole-genome sequences of haematopoietic colonies from 7 individuals from an older age group, we reveal distinct patterns of lifelong clonal behaviour. DNMT3A-mutant clones preferentially expanded early in life and displayed slower growth in old age, in the context of an increasingly competitive oligoclonal landscape. By contrast, splicing gene mutations drove expansion only later in life, whereas TET2-mutant clones emerged across all ages. Finally, we show that mutations driving faster clonal growth carry a higher risk of malignant progression. Our findings characterize the lifelong natural history of clonal haematopoiesis and give fundamental insights into the interactions between somatic mutation, ageing and clonal selection., (© 2022. The Author(s).)

Advances in cancer genomics have revealed genomic classes of acute myeloid leukemia (AML) characterized by class-defining mutations, such as chimeric fusion genes or in genes such as NPM1, MLL, and CEBPA. These class-defining mutations frequently synergize with internal tandem duplications in FLT3 (FLT3-ITDs) to drive leukemogenesis. However, ∼20% of FLT3-ITD-positive AMLs bare no class-defining mutations, and mechanisms of leukemic transformation in these cases are unknown. To identify pathways that drive FLT3-ITD mutant AML in the absence of class-defining mutations, we performed an insertional mutagenesis (IM) screening in Flt3-ITD mice, using Sleeping Beauty transposons. All mice developed acute leukemia (predominantly AML) after a median of 73 days. Analysis of transposon insertions in 38 samples from Flt3-ITD/IM leukemic mice identified recurrent integrations at 22 loci, including Setbp1 (20/38), Ets1 (11/38), Ash1l (8/38), Notch1 (8/38), Erg (7/38), and Runx1 (5/38). Insertions at Setbp1 led exclusively to AML and activated a transcriptional program similar, but not identical, to those of NPM1-mutant and MLL-rearranged AMLs. Guide RNA targeting of Setbp1 was highly detrimental to Flt3ITD/+/Setbp1IM+, but not to Flt3ITD/+/Npm1cA/+, AMLs. Also, analysis of RNA-sequencing data from hundreds of human AMLs revealed that SETBP1 expression is significantly higher in FLT3-ITD AMLs lacking class-defining mutations. These findings propose that SETBP1 overexpression collaborates with FLT3-ITD to drive a subtype of human AML. To identify genetic vulnerabilities of these AMLs, we performed genome-wide CRISPR-Cas9 screening in Flt3ITD/+/Setbp1IM+ AMLs and identified potential therapeutic targets, including Kdm1a, Brd3, Ezh2, and Hmgcr. Our study gives new insights into epigenetic pathways that can drive AMLs lacking class-defining mutations and proposes therapeutic approaches against such cases., (© 2021 by The American Society of Hematology.)

Background: Glioma is the most common intrinsic brain tumor and also occurs in the spinal cord. Activating EGFR mutations are common in IDH1 wild-type gliomas. However, the cooperative partners of EGFR driving gliomagenesis remain poorly understood., Results: We explore EGFR-mutant glioma evolution in conditional mutant mice by whole-exome sequencing, transposon mutagenesis forward genetic screening, and transcriptomics. We show mutant EGFR is sufficient to initiate gliomagenesis in vivo, both in the brain and spinal cord. We identify significantly recurrent somatic alterations in these gliomas including mutant EGFR amplifications and Sub1, Trp53, and Tead2 loss-of-function mutations. Comprehensive functional characterization of 96 gliomas by genome-wide piggyBac insertional mutagenesis in vivo identifies 281 known and novel EGFR-cooperating driver genes, including Cdkn2a, Nf1, Spred1, and Nav3. Transcriptomics confirms transposon-mediated effects on expression of these genes. We validate the clinical relevance of new putative tumor suppressors by showing these are frequently altered in patients' gliomas, with prognostic implications. We discover shared and distinct driver mutations in brain and spinal gliomas and confirm in vivo differential tumor suppressive effects of Pten between these tumors. Functional validation with CRISPR-Cas9-induced mutations in novel genes Tead2, Spred1, and Nav3 demonstrates heightened EGFRvIII-glioma cell proliferation. Chemogenomic analysis of mutated glioma genes reveals potential drug targets, with several investigational drugs showing efficacy in vitro., Conclusion: Our work elucidates functional driver landscapes of EGFR-mutant gliomas, uncovering potential therapeutic strategies, and provides new tools for functional interrogation of gliomagenesis.

Although acquisition of leukemia-associated somatic mutations by 1 or more hematopoietic stem cells is inevitable with advancing age, its consequences are highly variable, ranging from clinically silent clonal hematopoiesis (CH) to leukemic progression. To investigate the influence of heritable factors on CH, we performed deep targeted sequencing of blood DNA from 52 monozygotic (MZ) and 27 dizygotic (DZ) twin pairs (aged 70-99 years). Using this highly sensitive approach, we identified CH (variant allele frequency ≥0.5%) in 62% of individuals. We did not observe higher concordance for CH within MZ twin pairs as compared with that within DZ twin pairs, or to that expected by chance. However, we did identify 2 MZ pairs in which both twins harbored identical rare somatic mutations, suggesting a shared cell of origin. Finally, in 3 MZ twin pairs harboring mutations in the same driver genes, serial blood samples taken 4 to 5 years apart showed substantial twin-to-twin variability in clonal trajectories. Our findings propose that the inherited genome does not exert a dominant influence on the behavior of adult CH and provide evidence that CH mutations may be acquired in utero., (© 2020 by The American Society of Hematology.)

Gene expression governs cell fate, and is regulated via a complex interplay of transcription factors and molecules that change chromatin structure. Advances in sequencing-based assays have enabled investigation of these processes genome-wide, leading to large datasets that combine information on the dynamics of gene expression, transcription factor binding and chromatin structure as cells differentiate. While numerous studies focus on the effects of these features on broader gene regulation, less work has been done on the mechanisms of gene-specific transcriptional control. In this study, we have focussed on the latter by integrating gene expression data for the in vitro differentiation of murine ES cells to macrophages and cardiomyocytes, with dynamic data on chromatin structure, epigenetics and transcription factor binding. Combining a novel strategy to identify communities of related control elements with a penalized regression approach, we developed individual models to identify the potential control elements predictive of the expression of each gene. Our models were compared to an existing method and evaluated using the existing literature and new experimental data from embryonic stem cell differentiation reporter assays. Our method is able to identify transcriptional control elements in a gene specific manner that reflect known regulatory relationships and to generate useful hypotheses for further testing., Competing Interests: The authors have declared that no competing interests exist.

Toxoplasma gondii is associated with physiological effects in the host. Dysregulation of catecholamines in the central nervous system has previously been observed in chronically infected animals. In the study described here, the noradrenergic system was found to be suppressed with decreased levels of norepinephrine (NE) in brains of infected animals and in infected human and rat neural cells in vitro The mechanism responsible for the NE suppression was found to be downregulation of dopamine β-hydroxylase (DBH) gene expression, encoding the enzyme that synthesizes norepinephrine from dopamine, with downregulation observed in vitro and in infected brain tissue, particularly in the dorsal locus coeruleus/pons region. The downregulation was sex specific, with males expressing reduced DBH mRNA levels whereas females were unchanged. Rather, DBH expression correlated with estrogen receptor in the female rat brains for this estrogen-regulated gene. DBH silencing was not a general response of neurons to infection, as human cytomegalovirus did not downregulate DBH expression. The noradrenergic-linked behaviors of sociability and arousal were altered in chronically infected animals, with a high correlation between DBH expression and infection intensity. A decrease in DBH expression in noradrenergic neurons can elevate dopamine levels, which provides a possible explanation for mixed observations of changes in this neurotransmitter with infection. Decreased NE is consistent with the loss of coordination and motor impairments associated with toxoplasmosis. Further, the altered norepinephrine synthesis observed here may, in part, explain behavioral effects of infection and associations with mental illness., (Copyright © 2019 American Society for Microbiology.)

We recently identified the splicing kinase gene SRPK1 as a genetic vulnerability of acute myeloid leukemia (AML). Here, we show that genetic or pharmacological inhibition of SRPK1 leads to cell cycle arrest, leukemic cell differentiation and prolonged survival of mice transplanted with MLL-rearranged AML. RNA-seq analysis demonstrates that SRPK1 inhibition leads to altered isoform levels of many genes including several with established roles in leukemogenesis such as MYB, BRD4 and MED24. We focus on BRD4 as its main isoforms have distinct molecular properties and find that SRPK1 inhibition produces a significant switch from the short to the long isoform at the mRNA and protein levels. This was associated with BRD4 eviction from genomic loci involved in leukemogenesis including BCL2 and MYC. We go on to show that this switch mediates at least part of the anti-leukemic effects of SRPK1 inhibition. Our findings reveal that SRPK1 represents a plausible new therapeutic target against AML.

Mutations in the gene encoding isocitrate dehydrogenase 2 (IDH2) occur in several types of cancer, including acute myeloid leukemia (AML). In model systems, mutant IDH2 causes hematopoietic differentiation arrest. Enasidenib, a selective small-molecule inhibitor of mutant IDH2, produces a clinical response in 40% of treated patients with relapsed/refractory AML by promoting leukemic cell differentiation. Here, we studied the clonal basis of response and acquired resistance to enasidenib treatment. Using sequential patient samples, we determined the clonal structure of hematopoietic cell populations at different stages of differentiation. Before therapy, IDH2-mutant clones showed variable differentiation arrest. Enasidenib treatment promoted hematopoietic differentiation from either terminal or ancestral mutant clones; less frequently, treatment promoted differentiation of nonmutant cells. Analysis of paired diagnosis/relapse samples did not identify second-site mutations in IDH2 at relapse. Instead, relapse arose by clonal evolution or selection of terminal or ancestral clones, thus highlighting multiple bypass pathways that could potentially be targeted to restore differentiation arrest. These results show how mapping of clonal structure in cell populations at different stages of differentiation can reveal the response and evolution of clones during treatment response and relapse.

In metazoans, tissues experiencing proteotoxic stress induce "transcellular chaperone signaling" (TCS) that activates molecular chaperones, such as hsp-90, in distal tissues. How this form of inter-tissue communication is mediated to upregulate systemic chaperone expression and whether it can be utilized to protect against protein misfolding diseases remain open questions. Using C. elegans, we identified key components of a systemic stress signaling pathway that links the innate immune response with proteostasis maintenance. We show that mild perturbation of proteostasis in the neurons or the intestine activates TCS via the GATA zinc-finger transcription factor PQM-1. PQM-1 coordinates neuron-activated TCS via the innate immunity-associated transmembrane protein CLEC-41, whereas intestine-activated TCS depends on the aspartic protease ASP-12. Both TCS pathways can induce hsp-90 in muscle cells and facilitate amelioration of Aβ 3-42 -associated toxicity. This may have powerful implications for the treatment of diseases related to proteostasis dysfunction., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

Underlying the development of malaria parasites within erythrocytes and the resulting pathogenicity is a hardwired program that secures proper timing of gene transcription and production of functionally relevant proteins. How stage-specific gene expression is orchestrated in vivo remains unclear. Here, using the assay for transposase accessible chromatin sequencing (ATAC-seq), we identified ∼4,000 regulatory regions in P. falciparum intraerythrocytic stages. The vast majority of these sites are located within 2 kb upstream of transcribed genes and their chromatin accessibility pattern correlates positively with abundance of the respective mRNA transcript. Importantly, these regions are sufficient to drive stage-specific reporter gene expression and DNA motifs enriched in stage-specific sets of regulatory regions interact with members of the P. falciparum AP2 transcription factor family. Collectively, this study provides initial insights into the in vivo gene regulatory network of P. falciparum intraerythrocytic stages and should serve as a valuable resource for future studies., (Copyright © 2018 The Author(s). Published by Elsevier Inc. All rights reserved.)

NPM1 mutations define the commonest subgroup of acute myeloid leukemia (AML) and frequently co-occur with FLT3 internal tandem duplications (ITD) or, less commonly, NRAS or KRAS mutations. Co-occurrence of mutant NPM1 with FLT3-ITD carries a significantly worse prognosis than NPM1-RAS combinations. To understand the molecular basis of these observations, we compare the effects of the 2 combinations on hematopoiesis and leukemogenesis in knock-in mice. Early effects of these mutations on hematopoiesis show that compound Npm1 cA/+ ;Nras G12D/+ or Npm1 cA ;Flt3 ITD share a number of features: Hox gene overexpression, enhanced self-renewal, expansion of hematopoietic progenitors, and myeloid differentiation bias. However, Npm1 cA ;Flt3 ITD mutants displayed significantly higher peripheral leukocyte counts, early depletion of common lymphoid progenitors, and a monocytic bias in comparison with the granulocytic bias in Npm1 cA/+ ;Nras G12D/+ mutants. Underlying this was a striking molecular synergy manifested as a dramatically altered gene expression profile in Npm1 cA ;Flt3 ITD , but not Npm1 cA/+ ;Nras G12D/+ , progenitors compared with wild-type. Both double-mutant models developed high-penetrance AML, although latency was significantly longer with Npm1 cA/+ ;Nras G12D/+ During AML evolution, both models acquired additional copies of the mutant Flt3 or Nras alleles, but only Npm1 cA/+ ;Nras G12D/+ mice showed acquisition of other human AML mutations, including IDH1 R132Q. We also find, using primary Cas9-expressing AMLs, that Hoxa genes and selected interactors or downstream targets are required for survival of both types of double-mutant AML. Our results show that molecular complementarity underlies the higher frequency and significantly worse prognosis associated with NPM1 c/ FLT3-ITD vs NPM1/NRAS-G12D- mutant AML and functionally confirm the role of HOXA genes in NPM1c-driven AML., (© 2017 by The American Society of Hematology.)

A number of real-world systems have common underlying patterns among them and deducing these patterns is important for us in order to understand the environment around us. These patterns in some instances are apparent upon observation while in many others especially those found in nature are well hidden. Moreover, the inherent stochasticity in these systems introduces sufficient noise that we need models capable to handling it in order to decipher the underlying pattern. Hidden Markov model (HMM) is a probabilistic model that is frequently used for studying the hidden patterns in an observed sequence or sets of observed sequences. Since its conception in the late 1960s it has been extensively applied in biology to capture patterns in various disciplines ranging from small DNA and protein molecules, their structure and architecture that forms the basis of life to multicellular levels such as movement analysis in humans. This chapter aims at a gentle introduction to the theory of HMM, the statistical problems usually associated with HMMs and their uses in biology.

Background: Transcriptional regulation of gene expression is essential for cellular differentiation and function, and defects in the process are associated with cancer. The ENCODE project has mapped potential regulatory sites across the complete genome in many cell types, and these regions have been shown to harbour many of the somatic mutations that occur in cancer cells, suggesting that their effects may drive cancer initiation and development. The ENCODE data suggests a very large number of regulatory sites, and methods are needed to identify those that are most relevant and to connect them to the genes that they control., Methods: Predictive models of gene expression were developed by integrating the ENCODE data for regulation, including transcription factor binding and DNase1 hypersensitivity, with RNA-seq data for gene expression. A penalized regression method was used to identify the most predictive potential regulatory sites for each transcript. Known cancer somatic mutations from the COSMIC database were mapped to potential regulatory sites, and we examined differences in the mapping frequencies associated with sites chosen in regulatory models and other (rejected) sites. The effects of potential confounders, for example replication timing, were considered., Results: Cancer somatic mutations preferentially occupy those regulatory regions chosen in our models as most predictive of gene expression., Conclusion: Our methods have identified a significantly reduced set of regulatory sites that are enriched in cancer somatic mutations and are more predictive of gene expression. This has significance for the mechanistic interpretation of cancer mutations, and the understanding of genetic regulation.

Metazoan development involves the successive activation and silencing of specific gene expression programs and is driven by tissue-specific transcription factors programming the chromatin landscape. To understand how this process executes an entire developmental pathway, we generated global gene expression, chromatin accessibility, histone modification, and transcription factor binding data from purified embryonic stem cell-derived cells representing six sequential stages of hematopoietic specification and differentiation. Our data reveal the nature of regulatory elements driving differential gene expression and inform how transcription factor binding impacts on promoter activity. We present a dynamic core regulatory network model for hematopoietic specification and demonstrate its utility for the design of reprogramming experiments. Functional studies motivated by our genome-wide data uncovered a stage-specific role for TEAD/YAP factors in mammalian hematopoietic specification. Our study presents a powerful resource for studying hematopoiesis and demonstrates how such data advance our understanding of mammalian development., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Dps (DNA-binding protein from starved cells) are dodecameric assemblies belonging to the ferritin family that can bind DNA, carry out ferroxidation, and store iron in their shells. The ferritin-like trimeric pore harbors the channel for the entry and exit of iron. By representing the structure of Dps as a network we have identified a charge-driven interface formed by a histidine aspartate cluster at the pore interface unique to Mycobacterium smegmatis Dps protein, MsDps2. Site-directed mutagenesis was employed to generate mutants to disrupt the charged interactions. Kinetics of iron uptake/release of the wild type and mutants were compared. Crystal structures were solved at a resolution of 1.8-2.2 Å for the various mutants to compare structural alterations vis à vis the wild type protein. The substitutions at the pore interface resulted in alterations in the side chain conformations leading to an overall weakening of the interface network, especially in cases of substitutions that alter the charge at the pore interface. Contrary to earlier findings where conserved aspartate residues were found crucial for iron release, we propose here that in the case of MsDps2, it is the interplay of negative-positive potentials at the pore that enables proper functioning of the protein. In similar studies in ferritins, negative and positive patches near the iron exit pore were found to be important in iron uptake/release kinetics. The unique ionic cluster in MsDps2 makes it a suitable candidate to act as nano-delivery vehicle, as these gated pores can be manipulated to exhibit conformations allowing for slow or fast rates of iron release.

Cell fate decisions during haematopoiesis are governed by lineage-specific transcription factors, such as RUNX1, SCL/TAL1, FLI1 and C/EBP family members. To gain insight into how these transcription factors regulate the activation of haematopoietic genes during embryonic development, we measured the genome-wide dynamics of transcription factor assembly on their target genes during the RUNX1-dependent transition from haemogenic endothelium (HE) to haematopoietic progenitors. Using a Runx1-/- embryonic stem cell differentiation model expressing an inducible Runx1 gene, we show that in the absence of RUNX1, haematopoietic genes bind SCL/TAL1, FLI1 and C/EBPβ and that this early priming is required for correct temporal expression of the myeloid master regulator PU.1 and its downstream targets. After induction, RUNX1 binds to numerous de novo sites, initiating a local increase in histone acetylation and rapid global alterations in the binding patterns of SCL/TAL1 and FLI1. The acquisition of haematopoietic fate controlled by Runx1 therefore does not represent the establishment of a new regulatory layer on top of a pre-existing HE program but instead entails global reorganization of lineage-specific transcription factor assemblies.

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or "sequence conservation" as the basis for their understanding. Recently "interaction energy" based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the "interaction conservation" viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

The three-dimensional structure of a protein is formed and maintained by the noncovalent interactions among the amino-acid residues of the polypeptide chain. These interactions can be represented collectively in the form of a network. So far, such networks have been investigated by considering the connections based on distances between the amino-acid residues. Here we present a method of constructing the structure network based on interaction energies among the amino-acid residues in the protein. We have investigated the properties of such protein energy-based networks (PENs) and have shown correlations to protein structural features such as the clusters of residues involved in stability, formation of secondary and super-secondary structural units. Further we demonstrate that the analysis of PENs in terms of parameters such as hubs and shortest paths can provide a variety of biologically important information, such as the residues crucial for stabilizing the folded units and the paths of communication between distal residues in the protein. Finally, the energy regimes for different levels of stabilization in the protein structure have clearly emerged from the PEN analysis., (Copyright © 2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The irreversible dodecamerization of native Dps trimers from Mycobacterium smegmatis, in vitro, is known to be directly associated with the bimodal function of this protein. Hence it is important to explore this pathway at the molecular level. Two types of trimers, Trimer A (tA) and Trimer B (tB), can be derived from the dodecamer due to the inherent 3-fold symmetry of the spherical crystal structure. These derived trimers were expressed as protein structure graphs (PSGs) using the computed interaction strength among the residues. Interface clusters which were identified from PSGs allowed us to convincingly predict E146 and F47 for further mutation studies. Various single and double mutants were constructed and characterized. We were finally able to generate a single mutant F47E impaired in dodecamerization and a double mutant E146AF47E as native monomer in solution. These two observed results suggest that the two trimers are important for dodecamerization and that the residues selected are important for the structural stability of the protein in vitro.

Protein-DNA interactions are crucial for many cellular processes. Now with the increased availability of structures of protein-DNA complexes, gaining deeper insights into the nature of protein-DNA interactions has become possible. Earlier, investigations have characterized the interface properties by considering pairwise interactions. However, the information communicated along the interfaces is rarely a pairwise phenomenon, and we feel that a global picture can be obtained by considering a protein-DNA complex as a network of noncovalently interacting systems. Furthermore, most of the earlier investigations have been carried out from the protein point of view (protein-centric), and the present network approach aims to combine both the protein-centric and the DNA-centric points of view. Part of the study involves the development of methodology to investigate protein-DNA graphs/networks with the development of key parameters. A network representation provides a holistic view of the interacting surface and has been reported here for the first time. The second part of the study involves the analyses of these graphs in terms of clusters of interacting residues and the identification of highly connected residues (hubs) along the protein-DNA interface. A predominance of deoxyribose-amino acid clusters in beta-sheet proteins, distinction of the interface clusters in helix-turn-helix, and the zipper-type proteins would not have been possible by conventional pairwise interaction analysis. Additionally, we propose a potential classification scheme for a set of protein-DNA complexes on the basis of the protein-DNA interface clusters. This provides a general idea of how the proteins interact with the different components of DNA in different complexes. Thus, we believe that the present graph-based method provides a deeper insight into the analysis of the protein-DNA recognition mechanisms by throwing more light on the nature and the specificity of these interactions.