Graph Theory and Networks in Biology

In this paper, we present a survey of the use of graph theoretical techniques in Biology. In particular, we discuss recent work on identifying and modelling the structure of bio-molecular networks, as well as the application of centrality measures to interaction networks and research on the hierarchical structure of such networks and network motifs. Work on the link between structural network properties and dynamics is also described, with emphasis on synchronization and disease propagation.Comment: 52 pages, 5 figures, Survey Pape

Gene expression patterns following unilateral traumatic brain injury reveals a local pro-inflammatory and remote anti-inflammatory response.

BackgroundTraumatic brain injury (TBI) results in irreversible damage at the site of impact and initiates cellular and molecular processes that lead to secondary neural injury in the surrounding tissue. We used microarray analysis to determine which genes, pathways and networks were significantly altered using a rat model of TBI. Adult rats received a unilateral controlled cortical impact (CCI) and were sacrificed 24 h post-injury. The ipsilateral hemi-brain tissue at the site of the injury, the corresponding contralateral hemi-brain tissue, and naïve (control) brain tissue were used for microarray analysis. Ingenuity Pathway Analysis (IPA) software was used to identify molecular pathways and networks that were associated with the altered gene expression in brain tissues following TBI.ResultsInspection of the top fifteen biological functions in IPA associated with TBI in the ipsilateral tissues revealed that all had an inflammatory component. IPA analysis also indicated that inflammatory genes were altered on the contralateral side, but many of the genes were inversely expressed compared to the ipsilateral side. The contralateral gene expression pattern suggests a remote anti-inflammatory molecular response. We created a network of the inversely expressed common (i.e., same gene changed on both sides of the brain) inflammatory response (IR) genes and those IR genes included in pathways and networks identified by IPA that changed on only one side. We ranked the genes by the number of direct connections each had in the network, creating a gene interaction hierarchy (GIH). Two well characterized signaling pathways, toll-like receptor/NF-kappaB signaling and JAK/STAT signaling, were prominent in our GIH.ConclusionsBioinformatic analysis of microarray data following TBI identified key molecular pathways and networks associated with neural injury following TBI. The GIH created here provides a starting point for investigating therapeutic targets in a ranked order that is somewhat different than what has been presented previously. In addition to being a vehicle for identifying potential targets for post-TBI therapeutic strategies, our findings can also provide a context for evaluating the potential of therapeutic agents currently in development

Automatic Network Fingerprinting through Single-Node Motifs

Complex networks have been characterised by their specific connectivity patterns (network motifs), but their building blocks can also be identified and described by node-motifs---a combination of local network features. One technique to identify single node-motifs has been presented by Costa et al. (L. D. F. Costa, F. A. Rodrigues, C. C. Hilgetag, and M. Kaiser, Europhys. Lett., 87, 1, 2009). Here, we first suggest improvements to the method including how its parameters can be determined automatically. Such automatic routines make high-throughput studies of many networks feasible. Second, the new routines are validated in different network-series. Third, we provide an example of how the method can be used to analyse network time-series. In conclusion, we provide a robust method for systematically discovering and classifying characteristic nodes of a network. In contrast to classical motif analysis, our approach can identify individual components (here: nodes) that are specific to a network. Such special nodes, as hubs before, might be found to play critical roles in real-world networks.Comment: 16 pages (4 figures) plus supporting information 8 pages (5 figures

Predicting gene expression in the human malaria parasite Plasmodium falciparum using histone modification, nucleosome positioning, and 3D localization features.

Empirical evidence suggests that the malaria parasite Plasmodium falciparum employs a broad range of mechanisms to regulate gene transcription throughout the organism's complex life cycle. To better understand this regulatory machinery, we assembled a rich collection of genomic and epigenomic data sets, including information about transcription factor (TF) binding motifs, patterns of covalent histone modifications, nucleosome occupancy, GC content, and global 3D genome architecture. We used these data to train machine learning models to discriminate between high-expression and low-expression genes, focusing on three distinct stages of the red blood cell phase of the Plasmodium life cycle. Our results highlight the importance of histone modifications and 3D chromatin architecture in Plasmodium transcriptional regulation and suggest that AP2 transcription factors may play a limited regulatory role, perhaps operating in conjunction with epigenetic factors

Regulation of Microtubule Dynamics in Axon Regeneration: Insights from C. elegans.

The capacity of an axon to regenerate is regulated by its external environment and by cell-intrinsic factors. Studies in a variety of organisms suggest that alterations in axonal microtubule (MT) dynamics have potent effects on axon regeneration. We review recent findings on the regulation of MT dynamics during axon regeneration, focusing on the nematode Caenorhabditis elegans. In C. elegans the dual leucine zipper kinase (DLK) promotes axon regeneration, whereas the exchange factor for Arf6 (EFA-6) inhibits axon regeneration. Both DLK and EFA-6 respond to injury and control axon regeneration in part via MT dynamics. How the DLK and EFA-6 pathways are related is a topic of active investigation, as is the mechanism by which EFA-6 responds to axonal injury. We evaluate potential candidates, such as the MT affinity-regulating kinase PAR-1/MARK, in regulation of EFA-6 and axonal MT dynamics in regeneration

Quantitative profiling of selective Sox/POU pairing on hundreds of sequences in parallel by Coop-seq

© The Author(s) 2016. Published by Oxford University Press on behalf of Nucleic Acids Research. Cooperative binding of transcription factors is known to be important in the regulation of gene expression programs conferring cellular identities. However, current methods to measure cooperativity parameters have been laborious and therefore limited to studying only a few sequence variants at a time. We developed Coop-seq (cooperativity by sequencing) that is capable of efficiently and accurately determining the cooperativity parameters for hundreds of different DNA sequences in a single experiment. We apply Coop-seq to 12 dimer pairs from the Sox and POU families of transcription factors using 324 unique sequences with changed half-site orientation, altered spacing and discrete randomization within the binding elements. The study reveals specific dimerization profiles of different Sox factors with Oct4. By contrast, Oct4 and the three neural class III POU factors Brn2, Brn4 and Oct6 assemble with Sox2 in a surprisingly indistinguishable manner. Two novel half-site configurations can support functional Sox/Oct dimerization in addition to known composite motifs. Moreover, Coop-seq uncovers a nucleotide switch within the POU half-site when spacing is altered, which is mirrored in genomic loci bound by Sox2/Oct4 complexes.Link_to_subscribed_fulltex

Identification of a selective G1-phase benzimidazolone inhibitor by a senescence-targeted virtual screen using artificial neural networks

Cellular senescence is a barrier to tumorigenesis in normal cells and tumour cells undergo senescence responses to genotoxic stimuli, which is a potential target phenotype for cancer therapy. However, in this setting, mixed-mode responses are common with apoptosis the dominant effect. Hence, more selective senescence inducers are required. Here we report a machine learning-based in silico screen to identify potential senescence agonists. We built profiles of differentially affected biological process networks from expression data obtained under induced telomere dysfunction conditions in colorectal cancer cells and matched these to a panel of 17 protein targets with confirmatory screening data in PubChem. We trained a neural network using 3517 compounds identified as active or inactive against these targets. The resulting classification model was used to screen a virtual library of ~2M lead-like compounds. 147 virtual hits were acquired for validation in growth inhibition and senescence-associated β-galactosidase (SA-β-gal) assays. Among the found hits a benzimidazolone compound, CB-20903630, had low micromolar IC50 for growth inhibition of HCT116 cells and selectively induced SA-β-gal activity in the entire treated cell population without cytotoxicity or apoptosis induction. Growth suppression was mediated by G1 blockade involving increased p21 expression and suppressed cyclin B1, CDK1 and CDC25C. Additionally, the compound inhibited growth of multicellular spheroids and caused severe retardation of population kinetics in long term treatments. Preliminary structure-activity and structure clustering analyses are reported and expression analysis of CB-20903630 against other cell cycle suppressor compounds suggested a PI3K/AKT-inhibitor-like profile in normal cells, with different pathways affected in cancer cells