Community Assessment of the Predictability of Cancer Protein and Phosphoprotein Levels from Genomics and Transcriptomics.

Cancer is driven by genomic alterations, but the processes causing this disease are largely performed by proteins. However, proteins are harder and more expensive to measure than genes and transcripts. To catalyze developments of methods to infer protein levels from other omics measurements, we leveraged crowdsourcing via the NCI-CPTAC DREAM proteogenomic challenge. We asked for methods to predict protein and phosphorylation levels from genomic and transcriptomic data in cancer patients. The best performance was achieved by an ensemble of models, including as predictors transcript level of the corresponding genes, interaction between genes, conservation across tumor types, and phosphosite proximity for phosphorylation prediction. Proteins from metabolic pathways and complexes were the best and worst predicted, respectively. The performance of even the best-performing model was modest, suggesting that many proteins are strongly regulated through translational control and degradation. Our results set a reference for the limitations of computational inference in proteogenomics. A record of this paper's transparent peer review process is included in the Supplemental Information.

Natural Compounds as Inhibitors of Tyrosyl-tRNA Synthetase.

Tyrosyl-tRNA synthetases (TyrRSs) as essential enzymes for all living organisms are good candidates for therapeutic target in the prevention and therapy of microbial infection. We examined the effect of various polyphenols, alkaloids, and terpenes-secondary metabolites produced by higher plants showing many beneficial properties for the human organism, on bacterial aminoacylation reaction. The most potent inhibitors of Escherichia coli TyrRS are epigallocatechin gallate, acacetin, kaempferide, and chrysin, whereas the enzymes from Staphylococcus aureus and Pseudomonas aeruginosa are inhibited mainly by acacetin and chrysin. Most of them act as competitive inhibitors. Structure-activity relationship showed that the most potent flavonoid inhibitors contain hydroxyl group at position 5 and 7 of A ring and OCH3 group at position 4' of B ring.

Computational characterisation of cancer molecular profiles derived using next generation sequencing.

Our current understanding of cancer genetics is grounded on the principle that cancer arises from a clone that has accumulated the requisite somatically acquired genetic aberrations, leading to the malignant transformation. It also results in aberrent of gene and protein expression. Next generation sequencing (NGS) or deep sequencing platforms are being used to create large catalogues of changes in copy numbers, mutations, structural variations, gene fusions, gene expression, and other types of information for cancer patients. However, inferring different types of biological changes from raw reads generated using the sequencing experiments is algorithmically and computationally challenging. In this article, we outline common steps for the quality control and processing of NGS data. We highlight the importance of accurate and application-specific alignment of these reads and the methodological steps and challenges in obtaining different types of information. We comment on the importance of integrating these data and building infrastructure to analyse it. We also provide exhaustive lists of available software to obtain information and point the readers to articles comparing software for deeper insight in specialised areas. We hope that the article will guide readers in choosing the right tools for analysing oncogenomic datasets.

Predicting early brain metastases based on clinicopathological factors and gene expression analysis in advanced HER2-positive breast cancer patients.

The overexpression or amplification of the human epidermal growth factor receptor 2 gene (HER2/neu) is associated with high risk of brain metastasis (BM). The identification of patients at highest immediate risk of BM could optimize screening and facilitate interventional trials. We performed gene expression analysis using complementary deoxyribonucleic acid-mediated annealing, selection, extension and ligation and real-time quantitative reverse transcription PCR (qRT-PCR) in primary tumor samples from two independent cohorts of advanced HER2 positive breast cancer patients. Additionally, we analyzed predictive relevance of clinicopathological factors in this series. Study group included discovery Cohort A (84 patients) and validation Cohort B (75 patients). The only independent variables associated with the development of early BM in both cohorts were the visceral location of first distant relapse [Cohort A: hazard ratio (HR) 7.4, 95 % CI 2.4-22.3; p < 0.001; Cohort B: HR 6.1, 95 % CI 1.5-25.6; p = 0.01] and the lack of trastuzumab administration in the metastatic setting (Cohort A: HR 5.0, 95 % CI 1.4-10.0; p = 0.009; Cohort B: HR 10.0, 95 % CI 2.0-100.0; p = 0.008). A profile including 13 genes was associated with early (≤36 months) symptomatic BM in the discovery cohort. This was refined by qRT-PCR to a 3-gene classifier (RAD51, HDGF, TPR) highly predictive of early BM (HR 5.3, 95 % CI 1.6-16.7; p = 0.005; multivariate analysis). However, predictive value of the classifier was not confirmed in the independent validation Cohort B. The presence of visceral metastases and the lack of trastuzumab administration in the metastatic setting apparently increase the likelihood of early BM in advanced HER2-positive breast cancer.

Smooth muscle contamination analysis in clinical oncology gene expression research.

Gene expression profiling is one of the most explored methods for studying cancers and microarray data repositories have become a rich and important resource. The most common human cancers develop in organs that are walled by smooth muscles. The only method of sample extraction free of unintentional contamination with surrounding tissue is microdissection. Nevertheless, such an approach is implemented infrequently. In the light of the above, there is a possibility of smooth muscle contamination in a large portion of publicly available data. In this study, 2292 publicly available microarrays were analysed to develop a simple screening method for detecting smooth muscle contamination. Microarray Inspector software was used to perform the tests since it has the unique ability to use many selected genes and probesets in a single group as a tissue definition. Furthermore, the test was dataset-independent. Two strategies of tissue definition were explored and compared. The first one depended on Tissue Specific Genes Database (TiSGeD) and BioGPS web resources, which themselves were based on meta-analysis of thousands of microarrays. The second method was based on a differential gene expression analysis of a few hundred preselected arrays. The comparison of the two methods proved the latter to be superior. Among the tested samples of undefined contamination, nearly half were identified to possibly contain significant smooth muscle traces. The obtained results equip researches with a simple method of examining microarray data for smooth muscle contamination. The presented work serves as an example of how to create definitions when searching for other possible contaminations.

Double site saturation mutagenesis of the human cytochrome P450 2D6 results in regioselective steroid hydroxylation.

The human cytochrome P450 2D6 (CYP2D6) is one of the major human drug metabolizing enzymes and acts preferably on substrates containing a basic nitrogen atom. Testosterone - just as other steroids - is an atypical substrate and only poorly metabolized by CYP2D6. The present study intended to investigate the influence of the two active site residues 216 and 483 on the capability of CYP2D6 to hydroxylate steroids such as for example testosterone. All 400 possible combinatorial mutations at these two positions have been generated and expressed individually in Pichia pastoris. Employing whole-cell biotransformations coupled with HPLC-MS analysis the testosterone hydroxylase activity and regioselectivity of every single CYP2D6 variant was determined. Covering the whole sequence space, CYP2D6 variants with improved activity and so far unknown regio-preference in testosterone hydroxylation were identified. Most intriguingly and in contrast to previous literature reports about mutein F483I, the mutation F483G led to preferred hydroxylation at the 2ß-position, while the slow formation of 6ß-hydroxytestosterone, the main product of wild-type CYP2D6, was further reduced. Two point mutations have already been sufficient to convert CYP2D6 into a steroid hydroxylase with the highest ever reported testosterone hydroxylation rate for this enzyme, which is of the same order of magnitude as for the conversion of the standard substrate bufuralol by wild-type CYP2D6. Furthermore, this study is also an example for efficient human CYP engineering in P. pastoris for biocatalytic applications and to study so far unknown pharmacokinetic effects of individual and combined mutations in these key enzymes of the human drug metabolism.

Microarray Inspector: tissue cross contamination detection tool for microarray data.

Microarray technology changed the landscape of contemporary life sciences by providing vast amounts of expression data. Researchers are building up repositories of experiment results with various conditions and samples which serve the scientific community as a precious resource. Ensuring that the sample is of high quality is of utmost importance to this effort. The task is complicated by the fact that in many cases datasets lack information concerning pre-experimental quality assessment. Transcription profiling of tissue samples may be invalidated by an error caused by heterogeneity of the material. The risk of tissue cross contamination is especially high in oncological studies, where it is often difficult to extract the sample. Therefore, there is a need of developing a method detecting tissue contamination in a post-experimental phase. We propose Microarray Inspector: customizable, user-friendly software that enables easy detection of samples containing mixed tissue types. The advantage of the tool is that it uses raw expression data files and analyses each array independently. In addition, the system allows the user to adjust the criteria of the analysis to conform to individual needs and research requirements. The final output of the program contains comfortable to read reports about tissue contamination assessment with detailed information about the test parameters and results. Microarray Inspector provides a list of contaminant biomarkers needed in the analysis of adipose tissue contamination. Using real data (datasets from public repositories) and our tool, we confirmed high specificity of the software in detecting contamination. The results indicated the presence of adipose tissue admixture in a range from approximately 4% to 13% in several tested surgical samples.

Mapping the substrate binding site of phenylacetone monooxygenase from Thermobifida fusca by mutational analysis.

Baeyer-Villiger monooxygenases catalyze oxidations that are of interest for biocatalytic applications. Among these enzymes, phenylacetone monooxygenase (PAMO) from Thermobifida fusca is the only protein showing remarkable stability. While related enzymes often present a broad substrate scope, PAMO accepts only a limited number of substrates. Due to the absence of a substrate in the elucidated crystal structure of PAMO, the substrate binding site of this protein has not yet been defined. In this study, a structural model of cyclopentanone monooxygenase, which acts on a broad range of compounds, has been prepared and compared with the structure of PAMO. This revealed 15 amino acid positions in the active site of PAMO that may account for its relatively narrow substrate specificity. We designed and analyzed 30 single and multiple mutants in order to verify the role of these positions. Extensive substrate screening revealed several mutants that displayed increased activity and altered regio- or enantioselectivity in Baeyer-Villiger reactions and sulfoxidations. Further substrate profiling resulted in the identification of mutants with improved catalytic properties toward synthetically attractive compounds. Moreover, the thermostability of the mutants was not compromised in comparison to that of the wild-type enzyme. Our data demonstrate that the positions identified within the active site of PAMO, namely, V54, I67, Q152, and A435, contribute to the substrate specificity of this enzyme. These findings will aid in more dedicated and effective redesign of PAMO and related monooxygenases toward an expanded substrate scope.

Detailed mechanism of squalene epoxidase inhibition by terbinafine.

Squalene epoxidase (SE) is a key flavin adenine dinucleotide (FAD)-dependent enzyme of ergosterol and cholesterol biosynthetic pathways and an attractive potential target for drugs used to inhibit the growth of pathogenic fungi or to lower cholesterol level. Although many studies on allylamine drugs activity have been published during the last 30 years, up until now no detailed mechanism of the squalene epoxidase inhibition has been presented. Our study brings such a model at atomic resolution in the case of yeast Saccharomyces cerevisiae . Presented data resulting from modeling studies are in excellent agreement with experimental findings. A fully atomic three-dimensional (3D) model of squalene epoxidase (EC 1.14.99.7) from S. cerevisiae was built with the help of 3D-Jury approach and further screened based on data known from mutation experiments leading to terbinafine resistance. Docking studies followed by molecular dynamics simulations and quantum interaction energy calculations [MP2/6-31G(d)] resulted in the identification of the terbinafine-squalene epoxidase mode of interaction. In the energetically most likely orientation of terbinafine its interaction energy with the protein is ca. 120 kJ/mol. In the favorable position the terbinafine lipophilic moiety is located vertically inside the squalene epoxidase binding pocket with the tert-butyl group oriented toward its center. Such a position results in the SE conformational changes and prevents the natural substrate from being able to bind to the enzyme's active site. That would explain the noncompetitive manner of SE inhibition. We found that the strongest interaction between terbinafine and SE stems from hydrogen bonding between hydrogen-bond donors, hydroxyl group of Tyr90 and amine nitrogen atom of terbinafine. Moreover, strong attractive interactions were recorded for amino acids whose mutations resulted in terbinafine resistance. Our results, elucidating at a molecular level the mode of terbinafine inhibitory activity, can be utilized in designing more potent or selective antifungal drugs or even medicines lowering cholesterol in humans.