Genome-scale metabolic model of the rat liver predicts effects of diet restriction.

Mapping network analysis in cells and tissues can provide insights into metabolic adaptations to changes in external environment, pathological conditions, and nutrient deprivation. Here, we reconstructed a genome-scale metabolic network of the rat liver that will allow for exploration of systems-level physiology. The resulting in silico model (iRatLiver) contains 1,882 reactions, 1,448 metabolites, and 994 metabolic genes. We then used this model to characterize the response of the liver\u27s energy metabolism to a controlled perturbation in diet. Transcriptomics data were collected from the livers of Sprague Dawley rats at 4 or 14 days of being subjected to 15%, 30%, or 60% diet restriction. These data were integrated with the iRatLiver model to generate condition-specific metabolic models, allowing us to explore network differences under each condition. We observed different pathway usage between early and late time points. Network analysis identified several highly connected hub genes (Pklr, Hadha, Tkt, Pgm1, Tpi1, and Eno3) that showed differing trends between early and late time points. Taken together, our results suggest that the liver\u27s response varied with short- and long-term diet restriction. More broadly, we anticipate that the iRatLiver model can be exploited further to study metabolic changes in the liver under other conditions such as drug treatment, infection, and disease

Structural effects of the highly protective V127 polymorphism on human prion protein

Prion diseases, a group of incurable, lethal neurodegenerative disorders of mammals including humans, are caused by prions, assemblies of misfolded host prion protein (PrP). A single point mutation (G127V) in human PrP prevents prion disease, however the structural basis for its protective effect remains unknown. Here we show that the mutation alters and constrains the PrP backbone conformation preceding the PrP β-sheet, stabilising PrP dimer interactions by increasing intermolecular hydrogen bonding. It also markedly changes the solution dynamics of the β2-α2 loop, a region of PrP structure implicated in prion transmission and cross-species susceptibility. Both of these structural changes may affect access to protein conformers susceptible to prion formation and explain its profound effect on prion disease

Resolving the Ortholog Conjecture: Orthologs Tend to Be Weakly, but Significantly, More Similar in Function than Paralogs

The function of most proteins is not determined experimentally, but is extrapolated from homologs. According to the “ortholog conjecture”, or standard model of phylogenomics, protein function changes rapidly after duplication, leading to paralogs with different functions, while orthologs retain the ancestral function. We report here that a comparison of experimentally supported functional annotations among homologs from 13 genomes mostly supports this model. We show that to analyze GO annotation effectively, several confounding factors need to be controlled: authorship bias, variation of GO term frequency among species, variation of background similarity among species pairs, and propagated annotation bias. After controlling for these biases, we observe that orthologs have generally more similar functional annotations than paralogs. This is especially strong for sub-cellular localization. We observe only a weak decrease in functional similarity with increasing sequence divergence. These findings hold over a large diversity of species; notably orthologs from model organisms such as E. coli, yeast or mouse have conserved function with human proteins

An evaluation of urinary microRNA reveals a high sensitivity for bladder cancer

Background: Urinary biomarkers are needed to improve the care and reduce the cost of managing bladder cancer. Current biomarkers struggle to identify both high and low-grade cancers due to differing molecular pathways. Changes in microRNA (miR) expression are seen in urothelial carcinogenesis in a phenotype-specific manner. We hypothesised that urinary miRs reflecting low- and high-grade pathways could detect bladder cancers and overcome differences in genetic events seen within the disease. Methods: We investigated urinary samples (n ¼ 121) from patients with bladder cancer (n ¼ 68) and age-matched controls (n ¼ 53). Fifteen miRs were quantified using real-time PCR. Results: We found that miR is stable within urinary cells despite adverse handling and detected differential expression of 10 miRs from patients with cancer and controls (miRs 15a/15b/24-1/27b/100/135b/203/212/328/1224, ANOVA Po0.05). Individually, miR-1224-3p had the best individual performance with specificity, positive and negative predictive values and concordance of 83%, 83%, 75% and 77%, respectively. The combination of miRs-135b/15b/1224-3p detected bladder cancer with a high sensitivity (94.1%), sufficient specificity (51%) and was correct in 86% of patients (concordance). Conclusion: The use of this panel in patients with haematuria would have found 94% of urothelial cell carcinoma, while reducing cystoscopy rates by 26%. However, two invasive cancers (3%) would have been missed

Small molecule FGF receptor inhibitors block FGFR-dependent urothelial carcinoma growth in vitro and in vivo

BACKGROUND: Activating mutations of FGFR3 are frequently identified in superficial urothelial carcinoma (UC) and increased expression of FGFR1 and FGFR3 are common in both superficial and invasive UC. METHODS: The effects of inhibition of receptor activity by three small molecule inhibitors (PD173074, TKI-258 and SU5402) were investigated in a panel of bladder tumour cell lines with known FGFR expression levels and FGFR3 mutation status. RESULTS: All inhibitors prevented activation of FGFR3, and inhibited downstream MAPK pathway signalling. Response was related to FGFR3 and/or FGFR1 expression levels. Cell lines with the highest levels of FGFR expression showed the greatest response and little or no effect was measured in normal human urothelial cells or in UC cell lines with activating RAS gene mutations. In sensitive cell lines, the drugs induced cell cycle arrest and/or apoptosis. IC(50) values for PD173074 and TKI-258 were in the nanomolar concentration range compared with micromolar concentrations for SU5402. PD173074 showed the greatest effects in vitro and in vivo significantly delayed the growth of subcutaneous bladder tumour xenografts. CONCLUSION: These results indicate that inhibition of FGFR1 and wild-type or mutant FGFR3 may represent a useful therapeutic approach in patients with both non-muscle invasive and muscle invasive UC

Quantitative sequence-function relationships in proteins based on gene ontology

<p>Abstract</p> <p>Background</p> <p>The relationship between divergence of amino-acid sequence and divergence of function among homologous proteins is complex. The assumption that homologs share function – the basis of transfer of annotations in databases – must therefore be regarded with caution. Here, we present a quantitative study of sequence and function divergence, based on the Gene Ontology classification of function. We determined the relationship between sequence divergence and function divergence in 6828 protein families from the PFAM database. Within families there is a broad range of sequence similarity from very closely related proteins – for instance, orthologs in different mammals – to very distantly-related proteins at the limit of reliable recognition of homology.</p> <p>Results</p> <p>We correlated the divergence in sequences determined from pairwise alignments, and the divergence in function determined by path lengths in the Gene Ontology graph, taking into account the fact that many proteins have multiple functions. Our results show that, among homologous proteins, the proportion of divergent functions decreases dramatically above a threshold of sequence similarity at about 50% residue identity. For proteins with more than 50% residue identity, transfer of annotation between homologs will lead to an erroneous attribution with a totally dissimilar function in fewer than 6% of cases. This means that for very similar proteins (about 50 % identical residues) the chance of completely incorrect annotation is low; however, because of the phenomenon of recruitment, it is still non-zero.</p> <p>Conclusion</p> <p>Our results describe general features of the evolution of protein function, and serve as a guide to the reliability of annotation transfer, based on the closeness of the relationship between a new protein and its nearest annotated relative.</p

Searching the protein structure database for ligand-binding site similarities using CPASS v.2

<p>Abstract</p> <p>Background</p> <p>A recent analysis of protein sequences deposited in the NCBI RefSeq database indicates that ~8.5 million protein sequences are encoded in prokaryotic and eukaryotic genomes, where ~30% are explicitly annotated as "hypothetical" or "uncharacterized" protein. Our Comparison of Protein Active-Site Structures (CPASS v.2) database and software compares the sequence and structural characteristics of experimentally determined ligand binding sites to infer a functional relationship in the absence of global sequence or structure similarity. CPASS is an important component of our Functional Annotation Screening Technology by NMR (FAST-NMR) protocol and has been successfully applied to aid the annotation of a number of proteins of unknown function.</p> <p>Findings</p> <p>We report a major upgrade to our CPASS software and database that significantly improves its broad utility. CPASS v.2 is designed with a layered architecture to increase flexibility and portability that also enables job distribution over the Open Science Grid (OSG) to increase speed. Similarly, the CPASS interface was enhanced to provide more user flexibility in submitting a CPASS query. CPASS v.2 now allows for both automatic and manual definition of ligand-binding sites and permits pair-wise, one versus all, one versus list, or list versus list comparisons. Solvent accessible surface area, ligand root-mean square difference, and Cβ distances have been incorporated into the CPASS similarity function to improve the quality of the results. The CPASS database has also been updated.</p> <p>Conclusions</p> <p>CPASS v.2 is more than an order of magnitude faster than the original implementation, and allows for multiple simultaneous job submissions. Similarly, the CPASS database of ligand-defined binding sites has increased in size by ~ 38%, dramatically increasing the likelihood of a positive search result. The modification to the CPASS similarity function is effective in reducing CPASS similarity scores for false positives by ~30%, while leaving true positives unaffected. Importantly, receiver operating characteristics (ROC) curves demonstrate the high correlation between CPASS similarity scores and an accurate functional assignment. As indicated by distribution curves, scores ≥ 30% infer a functional similarity. Software URL: <url>http://cpass.unl.edu</url>.</p