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

Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans

Genome-wide association studies (GWAS) have identified numerous common prostate cancer (PrCa) susceptibility loci. We have fine-mapped 64 GWAS regions known at the conclusion of the iCOGS study using large-scale genotyping and imputation in 25 723 PrCa cases and 26 274 controls of European ancestry. We detected evidence for multiple independent signals at 16 regions, 12 of which contained additional newly identified significant associations. A single signal comprising a spectrum of correlated variation was observed at 39 regions; 35 of which are now described by a novel more significantly associated lead SNP, while the originally reported variant remained as the lead SNP only in 4 regions. We also confirmed two association signals in Europeans that had been previously reported only in East-Asian GWAS. Based on statistical evidence and linkage disequilibrium (LD) structure, we have curated and narrowed down the list of the most likely candidate causal variants for each region. Functional annotation using data from ENCODE filtered for PrCa cell lines and eQTL analysis demonstrated significant enrichment for overlap with bio-features within this set. By incorporating the novel risk variants identified here alongside the refined data for existing association signals, we estimate that these loci now explain ∼38.9% of the familial relative risk of PrCa, an 8.9% improvement over the previously reported GWAS tag SNPs. This suggests that a significant fraction of the heritability of PrCa may have been hidden during the discovery phase of GWAS, in particular due to the presence of multiple independent signals within the same regio

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

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

Rare and Complex Urology: Clinical Overview of ERN eUROGEN

Background: In 2017, the European Commission launched 24 European Reference Networks (ERNs). ERN eUROGEN is the network for urorectogenital diseases and complex conditions, and started with 29 full member healthcare providers (HCPs) in 11 countries. It then covered 19 different disease areas distributed over three workstreams (WSs). Objective: To provide an overview and identify challenges in data collection at European level of the ERN eUROGEN patient population treated by HCPs in the network. Design, setting, and participants: A retrospective cohort study was conducted of the 29 HCPs who were full members between 2013 and 2019. Outcome measurements and statistical analysis: Data were extracted from the original HCP applications and the ERN continuous monitoring system. Patient volumes, new patient numbers, and procedures were compared between different WSs, countries, and HCPs. Discrepancies between monitoring and application data were identified. Results and limitations: Between 2013 and 2019, 122 040 patients required long-term care within the 29 HCPs. The volume of patients treated and procedures undertaken per year increased over time. Large discrepancies were found between patient numbers contained in the application forms and those reported in the continuous monitoring system (0–1357% deviation). Conclusions: Patient numbers and procedures increased across ERN eUROGEN HCPs. Reliable data extraction appeared challenging, illustrated by the patient volume discrepancies between application forms and the continuous monitoring data. Improved disease definitions, re-evaluation of affiliated HCPs, and valid data extraction are needed for future improvements. Patient summary: We analysed the patient population with rare urorectogenital diseases or complex conditions within the ERN eUROGEN network between 2013 and 2019. Clinical activity was found to increase, but differences in patient numbers were evident between healthcare providers. In order to acquire valid patient numbers, both improved definitions of diagnostic codes and greater insight into the data-gathering process are required