A text-mining system for extracting metabolic reactions from full-text articles

Background: Increasingly biological text mining research is focusing on the extraction of complex relationships relevant to the construction and curation of biological networks and pathways. However, one important category of pathway—metabolic pathways—has been largely neglected. Here we present a relatively simple method for extracting metabolic reaction information from free text that scores different permutations of assigned entities (enzymes and metabolites) within a given sentence based on the presence and location of stemmed keywords. This method extends an approach that has proved effective in the context of the extraction of protein–protein interactions. Results: When evaluated on a set of manually-curated metabolic pathways using standard performance criteria, our method performs surprisingly well. Precision and recall rates are comparable to those previously achieved for the well-known protein-protein interaction extraction task. Conclusions: We conclude that automated metabolic pathway construction is more tractable than has often been assumed, and that (as in the case of protein–protein interaction extraction) relatively simple text-mining approaches can prove surprisingly effective. It is hoped that these results will provide an impetus to further research and act as a useful benchmark for judging the performance of more sophisticated methods that are yet to be developed

Using structural motifs to identify proteins with DNA binding function

This work describes a method for predicting DNA binding function from structure using 3-dimensional templates. Proteins that bind DNA using small contiguous helix¿turn¿helix (HTH) motifs comprise a significant number of all DNA-binding proteins. A structural template library of seven HTH motifs has been created from non-homologous DNA-binding proteins in the Protein Data Bank. The templates were used to scan complete protein structures using an algorithm that calculated the root mean squared deviation (rmsd) for the optimal superposition of each template on each structure, based on Ca backbone coordinates. Distributions of rmsd values for known HTH-containing proteins (true hits) and non-HTH proteins (false hits) were calculated. A threshold value of 1.6 Å rmsd was selected that gave a true hit rate of 88.4% and a false positive rate of 0.7%. The false positive rate was further reduced to 0.5% by introducing an accessible surface area threshold value of 990 Å2 per HTH motif. The template library and the validated thresholds were used to make predictions for target proteins from a structural genomics project

Computational modelling of the binding of arachidonic acid to the human monooxygenase CYP2J2

An experimentally determined structure for human CYP2J2—a member of the cytochrome P450 family with significant and diverse roles across a number of tissues—does not yet exist. Our understanding of how CYP2J2 accommodates its cognate substrates and how it might be inhibited by other ligands thus relies on our ability to computationally predict such interactions using modelling techniques. In this study we present a computational investigation of the binding of arachidonic acid (AA) to CYP2J2 using homology modelling, induced fit docking (IFD) and molecular dynamics (MD) simulations. Our study reveals a catalytically competent binding mode for AA that is distinct from a recently published study that followed a different computational pipeline. Our proposed binding mode for AA is supported by crystal structures of complexes of related enzymes to inhibitors, and evolutionary conservation of a residue whose role appears essential for placing AA in the right site for catalysis

Untranslated Parts of Genes Interpreted: Making Heads or Tails of High-Throughput Transcriptomic Data via Computational Methods Computational methods to discover and quantify isoforms with alternative untranslated regions

In this review we highlight the importance of defining the untranslated parts of transcripts, and present a number of computational approaches for the discovery and quantification of alternative transcription start and poly‐adenylation events in high‐throughput transcriptomic data. The fate of eukaryotic transcripts is closely linked to their untranslated regions, which are determined by the position at which transcription starts and ends at a genomic locus. Although the extent of alternative transcription starts and alternative poly‐adenylation sites has been revealed by sequencing methods focused on the ends of transcripts, the application of these methods is not yet widely adopted by the community. We suggest that computational methods applied to standard high‐throughput technologies are a useful, albeit less accurate, alternative to the expertise‐demanding 5′ and 3′ sequencing and they are the only option for analysing legacy transcriptomic data. We review these methods here, focusing on technical challenges and arguing for the need to include better normalization of the data and more appropriate statistical models of the expected variation in the signal

KSHV SOX mediated host shutoff: the molecular mechanism underlying mRNA transcript processing

Onset of the lytic phase in the KSHV life cycle is accompanied by the rapid, global degradation of host (and viral) mRNA transcripts in a process termed host shutoff. Key to this destruction is the virally encoded alkaline exonuclease SOX. While SOX has been shown to possess an intrinsic RNase activity and a potential consensus sequence for endonucleolytic cleavage identified, the structures of the RNA substrates targeted remained unclear. Based on an analysis of three reported target transcripts, we were able to identify common structures and confirm that these are indeed degraded by SOX in vitro as well as predict the presence of such elements in the KSHV pre-microRNA transcript K12-2. From these studies, we were able to determine the crystal structure of SOX productively bound to a 31 nucleotide K12-2 fragment. This complex not only reveals the structural determinants required for RNA recognition and degradation but, together with biochemical and biophysical studies, reveals distinct roles for residues implicated in host shutoff. Our results further confirm that SOX and the host exoribonuclease Xrn1 act in concert to elicit the rapid degradation of mRNA substrates observed in vivo, and that the activities of the two ribonucleases are co-ordinated

Leflunomide treatment does not protect neural cells following Oxygen-Glucose Deprivation (OGD) in vitro

Neonatal hypoxia-ischemia (HI) affects 2–3 per 1000 live births in developed countries and up to 26 per 1000 live births in developing countries. It is estimated that of the 750,000 infants experiencing a hypoxic-ischemic event during birth per year, more than 400,000 will be severely affected. As treatment options are limited, rapidly identifying new therapeutic avenues is critical, and repurposing drugs al- ready in clinical use offers a fast-track route to clinic. One emerging avenue for therapeutic intervention in neonatal HI is to target mitochondrial dysfunction, which occurs early in the development of brain injury. Mitochondrial dynamics are particularly affected, with mitochondrial fragmentation occurring at the expense of the pro-fusion protein Optic Atrophy (OPA)1. OPA1, together with mitofusins (MFN)1/2, are required for membrane fusion, and therefore, protecting their function may also safeguard mitochon- drial dynamics. Leflunomide, an FDA-approved immunosuppressant, was recently identified as an acti- vator of MFN2 with partial effects on OPA1 expression. We, therefore, treated C17.2 cells with Lefluno- mide before or after oxygen-glucose deprivation, an in vitro mimic of HI, to determine its efficacy as a neuroprotection and inhibitor of mitochondrial dysfunction. Leflunomide increased baseline OPA1 but not MFN2 expression in C17.2 cells. However, Leflunomide was unable to promote cell survival follow- ing OGD. Equally, there was no obvious effect on mitochondrial morphology or bioenergetics. These data align with studies suggesting that the tissue and mitochondrial protein profile of the target cell/tissue are critical for taking advantage of the therapeutic actions of Leflunomide

Dysregulation of alternative poly-adenylation as a potential player in Autism Spectrum Disorder

We present here the hypothesis that alternative poly-adenylation (APA) is dysregulated in the brains of individuals affected by Autism Spectrum Disorder (ASD), due to disruptions in the calcium signaling networks. APA, the process of selecting different poly-adenylation sites on the same gene, yielding transcripts with different-length 3′ untranslated regions (UTRs), has been documented in different tissues, stages of development and pathologic conditions. Differential use of poly-adenylation sites has been shown to regulate the function, stability, localization and translation efficiency of target RNAs. However, the role of APA remains rather unexplored in neurodevelopmental conditions. In the human brain, where transcripts have the longest 3′ UTRs and are thus likely to be under more complex post-transcriptional regulation, erratic APA could be particularly detrimental. In the context of ASD, a condition that affects individuals in markedly different ways and whose symptoms exhibit a spectrum of severity, APA dysregulation could be amplified or dampened depending on the individual and the extent of the effect on specific genes would likely vary with genetic and environmental factors. If this hypothesis is correct, dysregulated APA events might be responsible for certain aspects of the phenotypes associated with ASD. Evidence supporting our hypothesis is derived from standard RNA-seq transcriptomic data but we suggest that future experiments should focus on techniques that probe the actual poly-adenylation site (3′ sequencing). To address issues arising from the use of post-mortem tissue and low numbers of heterogeneous samples affected by confounding factors (such as the age, gender and health of the individuals), carefully controlled in vitro systems will be required to model the effect of calcium signaling dysregulation in the ASD brain

Visualisation of variable binding pockets on protein surfaces by probabilistic analysis of related structure sets

Background: Protein structures provide a valuable resource for rational drug design. For a protein with no known ligand, computational tools can predict surface pockets that are of suitable size and shape to accommodate a complementary small-molecule drug. However, pocket prediction against single static structures may miss features of pockets that arise from proteins’ dynamic behaviour. In particular, ligand-binding conformations can be observed as transiently populated states of the apo protein, so it is possible to gain insight into ligand-bound forms by considering conformational variation in apo proteins. This variation can be explored by considering sets of related structures: computationally generated conformers, solution NMR ensembles, multiple crystal structures, homologues or homology models. It is non-trivial to compare pockets, either from different programs or across sets of structures. For a single structure, difficulties arise in defining particular pocket’s boundaries. For a set of conformationally distinct structures the challenge is how to make reasonable comparisons between them given that a perfect structural alignment is not possible. Results: We have developed a computational method, Provar, that provides a consistent representation of predicted binding pockets across sets of related protein structures. The outputs are probabilities that each atom or residue of the protein borders a predicted pocket. These probabilities can be readily visualised on a protein using existing molecular graphics software. We show how Provar simplifies comparison of the outputs of different pocket prediction algorithms, of pockets across multiple simulated conformations and between homologous structures. We demonstrate the benefits of use of multiple structures for protein-ligand and protein-protein interface analysis on a set of complexes and consider three case studies in detail: i) analysis of a kinase superfamily highlights the conserved occurrence of surface pockets at the active and regulatory sites ii) a simulated ensemble of unliganded Bcl2 structures reveals extensions of a known ligand-binding pocket not apparent in the apo crystal structure; iii) visualisations of interleukin-2 and its homologues highlight conserved pockets at the known receptor interfaces and regions whose conformation is known to change on inhibitor binding. Conclusions: Through post-processing of the output of a variety of pocket prediction software, Provar provides a flexible approach to the analysis and visualization of the persistence or variability of pockets in sets of related protein structures

KSHV SOX mediated host shutoff: the molecular mechanism underlying mRNA transcript processing.

Onset of the lytic phase in the KSHV life cycle is accompanied by the rapid, global degradation of host (and viral) mRNA transcripts in a process termed host shutoff. Key to this destruction is the virally encoded alkaline exonuclease SOX. While SOX has been shown to possess an intrinsic RNase activity and a potential consensus sequence for endonucleolytic cleavage identified, the structures of the RNA substrates targeted remained unclear. Based on an analysis of three reported target transcripts, we were able to identify common structures and confirm that these are indeed degraded by SOX in vitro as well as predict the presence of such elements in the KSHV pre-microRNA transcript K12-2. From these studies, we were able to determine the crystal structure of SOX productively bound to a 31 nucleotide K12-2 fragment. This complex not only reveals the structural determinants required for RNA recognition and degradation but, together with biochemical and biophysical studies, reveals distinct roles for residues implicated in host shutoff. Our results further confirm that SOX and the host exoribonuclease Xrn1 act in concert to elicit the rapid degradation of mRNA substrates observed in vivo, and that the activities of the two ribonucleases are co-ordinated