Fr-TM-align: a new protein structural alignment method based on fragment alignments and the TM-score

©2008 Pandit and Skolnick; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This article is available from: http://www.biomedcentral.com/1471-2105/9/531doi:10.1186/1471-2105-9-531Background: Protein tertiary structure comparisons are employed in various fields of contemporary structural biology. Most structure comparison methods involve generation of an initial seed alignment, which is extended and/or refined to provide the best structural superposition between a pair of protein structures as assessed by a structure comparison metric. One such metric, the TM-score, was recently introduced to provide a combined structure quality measure of the coordinate root mean square deviation between a pair of structures and coverage. Using the TM-score, the TM-align structure alignment algorithm was developed that was often found to have better accuracy and coverage than the most commonly used structural alignment programs; however, there were a number of situations when this was not true. Results: To further improve structure alignment quality, the Fr-TM-align algorithm has been developed where aligned fragment pairs are used to generate the initial seed alignments that are then refined using dynamic programming to maximize the TM-score. For the assessment of the structural alignment quality from Fr-TM-align in comparison to other programs such as CE and TMalign, we examined various alignment quality assessment scores such as PSI and TM-score. The assessment showed that the structural alignment quality from Fr-TM-align is better in comparison to both CE and TM-align. On average, the structural alignments generated using Fr-TM-align have a higher TM-score (~9%) and coverage (~7%) in comparison to those generated by TM-align. Fr- TM-align uses an exhaustive procedure to generate initial seed alignments. Hence, the algorithm is computationally more expensive than TM-align. Conclusion: Fr-TM-align, a new algorithm that employs fragment alignment and assembly provides better structural alignments in comparison to TM-align. The source code and executables of Fr- TM-align are freely downloadable at: http://cssb.biology.gatech.edu/skolnick/files/FrTMalign/

The matrix model version of AGT conjecture and CIV-DV prepotential

Recently exact formulas were provided for partition function of conformal (multi-Penner) beta-ensemble in the Dijkgraaf-Vafa phase, which, if interpreted as Dotsenko-Fateev correlator of screenings and analytically continued in the number of screening insertions, represents generic Virasoro conformal blocks. Actually these formulas describe the lowest terms of the q_a-expansion, where q_a parameterize the shape of the Penner potential, and are exact in the filling numbers N_a. At the same time, the older theory of CIV-DV prepotential, straightforwardly extended to arbitrary beta and to non-polynomial potentials, provides an alternative expansion: in powers of N_a and exact in q_a. We check that the two expansions coincide in the overlapping region, i.e. for the lowest terms of expansions in both q_a and N_a. This coincidence is somewhat non-trivial, since the two methods use different integration contours: integrals in one case are of the B-function (Euler-Selberg) type, while in the other case they are Gaussian integrals.Comment: 27 pages, 1 figur

The genomic landscape of cutaneous SCC reveals drivers and a novel azathioprine associated mutational signature

Cutaneous squamous cell carcinoma (cSCC) has a high tumour mutational burden (50 mutations per megabase DNA pair). Here, we combine whole-exome analyses from 40 primary cSCC tumours, comprising 20 well-differentiated and 20 moderately/poorly differentiated tumours, with accompanying clinical data from a longitudinal study of immunosuppressed and immunocompetent patients and integrate this analysis with independent gene expression studies. We identify commonly mutated genes, copy number changes and altered pathways and processes. Comparisons with tumour differentiation status suggest events which may drive disease progression. Mutational signature analysis reveals the presence of a novel signature (signature 32), whose incidence correlates with chronic exposure to the immunosuppressive drug azathioprine. Characterisation of a panel of 15 cSCC tumour-derived cell lines reveals that they accurately reflect the mutational signatures and genomic alterations of primary tumours and provide a valuable resource for the validation of tumour drivers and therapeutic targets

Whole-genome sequencing of bladder cancers reveals somatic CDKN1A mutations and clinicopathological associations with mutation burden

Bladder cancers are a leading cause of death from malignancy. Molecular markers might predict disease progression and behaviour more accurately than the available prognostic factors. Here we use whole-genome sequencing to identify somatic mutations and chromosomal changes in 14 bladder cancers of different grades and stages. As well as detecting the known bladder cancer driver mutations, we report the identification of recurrent protein-inactivating mutations in CDKN1A and FAT1. The former are not mutually exclusive with TP53 mutations or MDM2 amplification, showing that CDKN1A dysfunction is not simply an alternative mechanism for p53 pathway inactivation. We find strong positive associations between higher tumour stage/grade and greater clonal diversity, the number of somatic mutations and the burden of copy number changes. In principle, the identification of sub-clones with greater diversity and/or mutation burden within early-stage or low-grade tumours could identify lesions with a high risk of invasive progression

Observation of a J^PC = 1-+ exotic resonance in diffractive dissociation of 190 GeV/c pi- into pi- pi- pi+

The COMPASS experiment at the CERN SPS has studied the diffractive dissociation of negative pions into the pi- pi- pi+ final state using a 190 GeV/c pion beam hitting a lead target. A partial wave analysis has been performed on a sample of 420000 events taken at values of the squared 4-momentum transfer t' between 0.1 and 1 GeV^2/c^2. The well-known resonances a1(1260), a2(1320), and pi2(1670) are clearly observed. In addition, the data show a significant natural parity exchange production of a resonance with spin-exotic quantum numbers J^PC = 1-+ at 1.66 GeV/c^2 decaying to rho pi. The resonant nature of this wave is evident from the mass-dependent phase differences to the J^PC = 2-+ and 1++ waves. From a mass-dependent fit a resonance mass of 1660 +- 10+0-64 MeV/c^2 and a width of 269+-21+42-64 MeV/c^2 is deduced.Comment: 7 page, 3 figures; version 2 gives some more details, data unchanged; version 3 updated authors, text shortened, data unchange

Normative productivity of the global vegetation

<p>Abstract</p> <p>Background</p> <p>The biosphere models of terrestrial productivity are essential for projecting climate change and assessing mitigation and adaptation options. Many of them have been developed in connection to the International Geosphere-Biosphere Program (IGBP) that backs the work of the Intergovernmental Panel on Climate Change (IPCC). In the end of 1990s, IGBP sponsored release of a data set summarizing the model outputs and setting certain norms for estimates of terrestrial productivity. Since a number of new models and new versions of old models were developed during the past decade, these normative data require updating.</p> <p>Results</p> <p>Here, we provide the series of updates that reflects evolution of biosphere models and demonstrates evolutional stability of the global and regional estimates of terrestrial productivity. Most of them fit well the long-living Miami model. At the same time we call attention to the emerging alternative: the global potential for net primary production of biomass may be as high as 70 PgC y^-1, the productivity of larch forest zone may be comparable to the productivity of taiga zone, and the productivity of rain-green forest zone may be comparable to the productivity of tropical rainforest zone.</p> <p>Conclusion</p> <p>The departure from Miami model's worldview mentioned above cannot be simply ignored. It requires thorough examination using modern observational tools and techniques for model-data fusion. Stability of normative knowledge is not its ultimate goal – the norms for estimates of terrestrial productivity must be evidence-based.</p

Amyloid-Mediated Sequestration of Essential Proteins Contributes to Mutant Huntingtin Toxicity in Yeast

BACKGROUND: Polyglutamine expansion is responsible for several neurodegenerative disorders, among which Huntington disease is the most well-known. Studies in the yeast model demonstrated that both aggregation and toxicity of a huntingtin (htt) protein with an expanded polyglutamine region strictly depend on the presence of the prion form of Rnq1 protein ([PIN+]), which has a glutamine/asparagine-rich domain. PRINCIPAL FINDINGS: Here, we showed that aggregation and toxicity of mutant htt depended on [PIN+] only quantitatively: the presence of [PIN+] elevated the toxicity and the levels of htt detergent-insoluble polymers. In cells lacking [PIN+], toxicity of mutant htt was due to the polymerization and inactivation of the essential glutamine/asparagine-rich Sup35 protein and related inactivation of another essential protein, Sup45, most probably via its sequestration into Sup35 aggregates. However, inhibition of growth of [PIN+] cells depended on Sup35/Sup45 depletion only partially, suggesting that there are other sources of mutant htt toxicity in yeast. CONCLUSIONS: The obtained data suggest that induced polymerization of essential glutamine/asparagine-rich proteins and related sequestration of other proteins which interact with these polymers represent an essential source of htt toxicity

Functional classification of proteins based on projection of amino acid sequences: application for prediction of protein kinase substrates

<p>Abstract</p> <p>Background</p> <p>The knowledge about proteins with specific interaction capacity to the protein partners is very important for the modeling of cell signaling networks. However, the experimentally-derived data are sufficiently not complete for the reconstruction of signaling pathways. This problem can be solved by the network enrichment with predicted protein interactions. The previously published <it>in silico </it>method PAAS was applied for prediction of interactions between protein kinases and their substrates.</p> <p>Results</p> <p>We used the method for recognition of the protein classes defined by the interaction with the same protein partners. 1021 protein kinase substrates classified by 45 kinases were extracted from the Phospho.ELM database and used as a training set. The reasonable accuracy of prediction calculated by leave-one-out cross validation procedure was observed in the majority of kinase-specificity classes. The random multiple splitting of the studied set onto the test and training set had also led to satisfactory results. The kinase substrate specificity for 186 proteins extracted from TRANSPATH^®database was predicted by PAAS method. Several kinase-substrate interactions described in this database were correctly predicted. Using the previously developed ExPlain™ system for the reconstruction of signal transduction pathways, we showed that addition of the newly predicted interactions enabled us to find the possible path between signal trigger, TNF-alpha, and its target genes in the cell.</p> <p>Conclusions</p> <p>It was shown that the predictions of protein kinase substrates by PAAS were suitable for the enrichment of signaling pathway networks and identification of the novel signaling pathways. The on-line version of PAAS for prediction of protein kinase substrates is freely available at <url>http://www.ibmc.msk.ru/PAAS/</url>.</p

Immunotherapy for head and neck squamous cell carcinoma

Over the past years immuno-oncology has evolved and become a novel promising strategy for cancer therapy. Immune checkpoint inhibitors such as pembrolizumab or nivolumab, which target the interaction between programmed death receptor 1/programmed death ligand 1 (PD-1/PDL-1) and PDL-2, have been recently approved for the treatment of various malignancies and are currently being investigated in clinical phase III trials for head and neck squamous cell carcinoma (HNSCC). Data available from these trials indicate substantial activity accompanied by a favorable safety and toxicity profile in this patient population. This review article focuses on the molecular background, gives an overview of current clinical data of checkpoint inhibitors in HNSCC, and points out future challenges such as the need for appropriate biomarkers for these novel compounds.(VLID)362080

Genomic architecture and evolution of clear cell renal cell carcinomas defined by multiregion sequencing

Clear cell renal carcinomas (ccRCCs) can display intratumor heterogeneity (ITH). We applied multiregion exome sequencing (M-seq) to resolve the genetic architecture and evolutionary histories of ten ccRCCs. Ultra-deep sequencing identified ITH in all cases. We found that 73–75% of identified ccRCC driver aberrations were subclonal, confounding estimates of driver mutation prevalence. ITH increased with the number of biopsies analyzed, without evidence of saturation in most tumors. Chromosome 3p loss and VHL aberrations were the only ubiquitous events. The proportion of C>T transitions at CpG sites increased during tumor progression. M-seq permits the temporal resolution of ccRCC evolution and refines mutational signatures occurring during tumor development