Neoepitopes prediction strategies: an integration of cancer genomics and immunoinformatics approaches.

A major near-term medical impact of the genomic technology revolution will be the elucidation of mechanisms of cancer pathogenesis, leading to improvements in the diagnosis of cancer and the selection of cancer treatment. Next-generation sequencing technologies have accelerated the characterization of a tumor, leading to the comprehensive discovery of all the major alterations in a given cancer genome, followed by the translation of this information using computational and immunoinformatics approaches to cancer diagnostics and therapeutic efforts. In the current article, we review various components of cancer immunoinformatics applied to a series of fields of cancer research, including computational tools for cancer mutation detection, cancer mutation and immunological databases, and computational vaccinology.

Immunodominant T-cell epitopes from the SARS-CoV-2 spike antigen reveal robust pre-existing T-cell immunity in unexposed individuals.

The COVID-19 pandemic has revealed a range of disease phenotypes in infected patients with asymptomatic, mild, or severe clinical outcomes, but the mechanisms that determine such variable outcomes remain unresolved. In this study, we identified immunodominant CD8 T-cell epitopes in the spike antigen using a novel TCR-binding algorithm. The predicted epitopes induced robust T-cell activation in unexposed donors demonstrating pre-existing CD4 and CD8 T-cell immunity to SARS-CoV-2 antigen. The T-cell reactivity to the predicted epitopes was higher than the Spike-S1 and S2 peptide pools in the unexposed donors. A key finding of our study is that pre-existing T-cell immunity to SARS-CoV-2 is contributed by TCRs that recognize common viral antigens such as Influenza and CMV, even though the viral epitopes lack sequence identity to the SARS-CoV-2 epitopes. This finding is in contrast to multiple published studies in which pre-existing T-cell immunity is suggested to arise from shared epitopes between SARS-CoV-2 and other common cold-causing coronaviruses. However, our findings suggest that SARS-CoV-2 reactive T-cells are likely to be present in many individuals because of prior exposure to flu and CMV viruses.

Multi-Omics Analysis to Characterize Cigarette Smoke Induced Molecular Alterations in Esophageal Cells.

Though smoking remains one of the established risk factors of esophageal squamous cell carcinoma, there is limited data on molecular alterations associated with cigarette smoke exposure in esophageal cells. To investigate molecular alterations associated with chronic exposure to cigarette smoke, non-neoplastic human esophageal epithelial cells were treated with cigarette smoke condensate (CSC) for up to 8 months. Chronic treatment with CSC increased cell proliferation and invasive ability of non-neoplastic esophageal cells. Whole exome sequence analysis of CSC treated cells revealed several mutations and copy number variations. This included loss of high mobility group nucleosomal binding domain 2 (HMGN2) and a missense variant in mediator complex subunit 1 (MED1). Both these genes play an important role in DNA repair. Global proteomic and phosphoproteomic profiling of CSC treated cells lead to the identification of 38 differentially expressed and 171 differentially phosphorylated proteins. Bioinformatics analysis of differentially expressed proteins and phosphoproteins revealed that most of these proteins are associated with DNA damage response pathway. Proteomics data revealed decreased expression of HMGN2 and hypophosphorylation of MED1. Exogenous expression of HMGN2 and MED1 lead to decreased proliferative and invasive ability of smoke exposed cells. Immunohistochemical labeling of HMGN2 in primary ESCC tumor tissue sections (from smokers) showed no detectable expression while strong to moderate staining of HMGN2 was observed in normal esophageal tissues. Our data suggests that cigarette smoke perturbs expression of proteins associated with DNA damage response pathways which might play a vital role in development of ESCC.

Mutational Landscape of Esophageal Squamous Cell Carcinoma in an Indian Cohort.

Esophageal squamous cell carcinoma (ESCC) is the most common histological subtype of esophageal cancer in India. Cigarette smoking and chewing tobacco are known risk factors associated with ESCC. However, genomic alterations associated with ESCC in India are not well-characterized. In this study, we carried out exome sequencing to characterize the mutational landscape of ESCC tumors from subjects with a varied history of tobacco usage. Whole exome sequence analysis of ESCC from an Indian cohort revealed several genes that were mutated or had copy number changes. ESCC from tobacco chewers had a higher frequency of C:G > A:T transversions and 2-fold enrichment for mutation signature 4 compared to smokers and non-users of tobacco. Genes, such as TP53, CSMD3, SYNE1, PIK3CA, and NOTCH1 were found to be frequently mutated in Indian cohort. Mutually exclusive mutation patterns were observed in PIK3CA-NOTCH1, DNAH5-ZFHX4, MUC16-FAT1, and ZFHX4-NOTCH1 gene pairs. Recurrent amplifications were observed in 3q22-3q29, 11q13.3-q13.4, 7q22.1-q31.1, and 8q24 regions. Approximately 53% of tumors had genomic alterations in PIK3CA making this pathway a promising candidate for targeted therapy. In conclusion, we observe enrichment of mutation signature 4 in ESCC tumors from patients with a history of tobacco chewing. This is likely due to direct exposure of esophagus to tobacco carcinogens when it is chewed and swallowed. Genomic alterations were frequently observed in PIK3CA-AKT pathway members independent of the history of tobacco usage. PIK3CA pathway can be potentially targeted in ESCC which currently has no effective targeted therapeutic options.

MAP2K1 is a potential therapeutic target in erlotinib resistant head and neck squamous cell carcinoma.

Epidermal growth factor receptor (EGFR) targeted therapies have shown limited efficacy in head and neck squamous cell carcinoma (HNSCC) patients despite its overexpression. Identifying molecular mechanisms associated with acquired resistance to EGFR-TKIs such as erlotinib remains an unmet need and a therapeutic challenge. In this study, we employed an integrated multi-omics approach to delineate mechanisms associated with acquired resistance to erlotinib by carrying out whole exome sequencing, quantitative proteomic and phosphoproteomic profiling. We observed amplification of several genes including AXL kinase and transcription factor YAP1 resulting in protein overexpression. We also observed expression of constitutively active mutant MAP2K1 (p.K57E) in erlotinib resistant SCC-R cells. An integrated analysis of genomic, proteomic and phosphoproteomic data revealed alterations in MAPK pathway and its downstream targets in SCC-R cells. We demonstrate that erlotinib-resistant cells are sensitive to MAPK pathway inhibition. This study revealed multiple genetic, proteomic and phosphoproteomic alterations associated with erlotinib resistant SCC-R cells. Our data indicates that therapeutic targeting of MAPK pathway is an effective strategy for treating erlotinib-resistant HNSCC tumors.

Chronic Exposure to Chewing Tobacco Induces Metabolic Reprogramming and Cancer Stem Cell-Like Properties in Esophageal Epithelial Cells.

Tobacco in its smoke and smokeless form are major risk factors for esophageal squamous cell carcinoma (ESCC). However, molecular alterations associated with smokeless tobacco exposure are poorly understood. In the Indian subcontinent, tobacco is predominantly consumed in chewing form. An understanding of molecular alterations associated with chewing tobacco exposure is vital for identifying molecular markers and potential targets. We developed an in vitro cellular model by exposing non-transformed esophageal epithelial cells to chewing tobacco over an eight-month period. Chronic exposure to chewing tobacco led to increase in cell proliferation, invasive ability and anchorage independent growth, indicating cell transformation. Molecular alterations associated with chewing tobacco exposure were characterized by carrying out exome sequencing and quantitative proteomic profiling of parental cells and chewing tobacco exposed cells. Quantitative proteomic analysis revealed increased expression of cancer stem cell markers in tobacco treated cells. In addition, tobacco exposed cells showed the Oxidative Phosphorylation (OXPHOS) phenotype with decreased expression of enzymes associated with glycolytic pathway and increased expression of a large number of mitochondrial proteins involved in electron transport chain as well as enzymes of the tricarboxylic acid (TCA) cycle. Electron micrographs revealed increase in number and size of mitochondria. Based on these observations, we propose that chronic exposure of esophageal epithelial cells to tobacco leads to cancer stem cell-like phenotype. These cells show the characteristic OXPHOS phenotype, which can be potentially targeted as a therapeutic strategy.

A cancer vaccine approach for personalized treatment of Lynch Syndrome.

Lynch syndrome (LS) is a cancer predisposition disorder wherein patients have a 70-80% lifetime risk of developing colorectal cancers (CRC). Finding germline mutations in predisposing genes allows for risk assessment of CRC development. Here we report a germline heterozygous frame-shift mutation in the mismatch repair MLH1 gene which was identified in members of two unrelated LS families. Since defects in DNA mismatch repair genes generate frame-shift mutations giving rise to highly immunogenic neoepitopes, we postulated that vaccination with these mutant peptide antigens could offer promising treatment options to LS patients. To this end we performed whole-exome and RNA seq analysis on the blood and tumour samples from an LS-CRC patient, and used our proprietary neoepitope prioritization pipeline OncoPeptVAC to select peptides, and confirm their immunogenicity in an ex vivo CD8+ T cell activation assay. Three neoepitopes derived from the tumour of this patient elicited a potent CD8+ T cell response. Furthermore, analysis of the tumour-associated immune infiltrate revealed CD8+ T cells expressing low levels of activation markers, suggesting mechanisms of immune suppression at play in this relapsed tumour. Taken together, our study paves the way towards development of a cancer vaccine to treat or delay the onset/relapse of LS-CRC.

Molecular alterations associated with chronic exposure to cigarette smoke and chewing tobacco in normal oral keratinocytes.

Tobacco usage is a known risk factor associated with development of oral cancer. It is mainly consumed in two different forms (smoking and chewing) that vary in their composition and methods of intake. Despite being the leading cause of oral cancer, molecular alterations induced by tobacco are poorly understood. We therefore sought to investigate the adverse effects of cigarette smoke/chewing tobacco exposure in oral keratinocytes (OKF6/TERT1). OKF6/TERT1 cells acquired oncogenic phenotype after treating with cigarette smoke/chewing tobacco for a period of 8 months. We employed whole exome sequencing (WES) and quantitative proteomics to investigate the molecular alterations in oral keratinocytes chronically exposed to smoke/ chewing tobacco. Exome sequencing revealed distinct mutational spectrum and copy number alterations in smoke/ chewing tobacco treated cells. We also observed differences in proteomic alterations. Proteins downstream of MAPK1 and EGFR were dysregulated in smoke and chewing tobacco exposed cells, respectively. This study can serve as a reference for fundamental damages on oral cells as a consequence of exposure to different forms of tobacco.

A Computational Approach Identifies Immunogenic Features of Prognosis in Human Cancers.

A large number of tumor intrinsic and extrinsic factors determine long-term survival in human cancers. In this study, we stratified 9120 tumors from 33 cancers with respect to their immune cell content and identified immunogenomic features associated with long-term survival. Our analysis demonstrates that tumors infiltrated by CD8+ T cells expressing higher levels of activation marker (PD1hi) along with TCR signaling genes and cytolytic T cell markers (IL2hi/TNF-αhi/IFN-γhi/GZMA-Bhi) extend survival, whereas survival benefit was absent for tumors infiltrated by anergic and hyperexhausted CD8+ T cells characterized by high expression of CTLA-4, TIM3, LAG3, and genes linked to PI3K signaling pathway. The computational approach of using robust and highly specific gene expression signatures to deconvolute the tumor microenvironment has important clinical applications, such as selecting patients who will benefit from checkpoint inhibitor treatment.

Classification of Protein Kinases Influenced by Conservation of Substrate Binding Residues.

With the advent of genome sequencing projects in the recent past, several kinases have come to light as regulating different signaling pathways. These kinases are generally classified into different subfamilies based on their sequence similarity with members of known subfamilies of kinases. A functional association is then defined to the kinase based on the subfamily to which it has been characterized. However, one of the key factors that give identity to a kinase in a subfamily is its ability to phosphorylate a given set of substrates. Substrate specificity of a kinase is largely determined by the residues at the substrate binding site. Though in general the sequence similarity based measure for classification more or less gives the preliminary idea on subfamily, understanding the molecular basis of kinase substrate recognition could further refine the classification scheme for kinases and render a better understanding of their functional role. In this analysis we emphasize on the possibility of using putative substrate binding information in the classification of a given kinase into a particular subfamily.

Sequence Analysis and Evolutionary Studies of Reelin Proteins.

The reelin gene is conserved across many vertebrate species, including humans. The protein product of this gene plays several important roles in early brain development and regulation of neural network plasticity of a matured brain structure. With an extended structure of 3461 amino acid sequences, consisting of eight reelin repeats, the human reelin sequence stands out as an exceptional model for evolutionary studies. In this study, sequence analysis of the human reelin and its homologues and reelin sequences from 104 other species is described in detail. Interesting sequence conservation patterns of individual repeats have been highlighted. Sequence phylogeny of the reelin sequences indicates a pattern similar to the evolution of the species, thereby serving as a highly conserved family for evolutionary purposes. Multiple sequence alignment of different reelin domain repeats, derived from homologues, suggests specific functions for individual repeats and high sequence conservation across reelin repeats from different organisms, albeit with few unusual domain architectures. A three-dimensional structural model of the full-length human reelin is now available that provides clues on residues at the dimer interface.

ECMIS: computational approach for the identification of hotspots at protein-protein interfaces.

BACKGROUND: Various methods have been developed to computationally predict hotspot residues at novel protein-protein interfaces. However, there are various challenges in obtaining accurate prediction. We have developed a novel method which uses different aspects of protein structure and sequence space at residue level to highlight interface residues crucial for the protein-protein complex formation. RESULTS: ECMIS (Energetic Conservation Mass Index and Spatial Clustering) algorithm was able to outperform existing hotspot identification methods. It was able to achieve around 80% accuracy with incredible increase in sensitivity and outperforms other existing methods. This method is even sensitive towards the hotspot residues contributing only small-scale hydrophobic interactions. CONCLUSION: Combination of diverse features of the protein viz. energy contribution, extent of conservation, location and surrounding environment, along with optimized weightage for each feature, was the key for the success of the algorithm. The academic version of the algorithm is available at http://caps.ncbs.res.in/download/ECMIS/ECMIS.zip.

Insights on pH-dependent conformational changes of mosquito odorant binding proteins by molecular dynamics simulations.

Chemical recognition plays an important role for the survival and reproduction of many insect species. Odorant binding proteins (OBPs) are the primary components of the insect olfactory mechanism and have been documented to play an important role in the host-seeking mechanism of mosquitoes. They are "transport proteins" believed to transport odorant molecules from the external environment to their respective membrane targets, the olfactory receptors. The mechanism by which this transport occurs in mosquitoes remains a conundrum in this field. Nevertheless, OBPs have proved to be amenable to conformational changes mediated by a pH change in other insect species. In this paper, the effect of pH on the conformational flexibility of mosquito OBPs is assessed computationally using molecular dynamics simulations of a mosquito OBP "CquiOBP1" bound to its pheromone 3OG (PDB ID: 3OGN). Conformational twist of a loop, driven by a set of well-characterized changes in intramolecular interactions of the loop, is demonstrated. The concomitant (i) closure of what is believed to be the entrance of the binding pocket, (ii) expansion of what could be an exit site, and (iii) migration of the ligand towards this putative exit site provide preliminary insights into the mechanism of ligand binding and release of these proteins in mosquitoes. The correlation of our results with previous experimental observations based on NMR studies help us provide a cardinal illustration on one of the probable dynamics and mechanism by which certain mosquito OBPs could deliver their ligand to their membrane-bound receptors at specific pH conditions.

Association of putative members to family of mosquito odorant binding proteins: scoring scheme using fuzzy functional templates and cys residue positions.

Proteins may be related to each other very specifically as homologous subfamilies. Proteins can also be related to diverse proteins at the super family level. It has become highly important to characterize the existing sequence databases by their signatures to facilitate the function annotation of newly added sequences. The algorithm described here uses a scheme for the classification of odorant binding proteins on the basis of functional residues and Cys-pairing. The cysteine-based scoring scheme not only helps in unambiguously identifying families like odorant binding proteins (OBPs), but also aids in their classification at the subfamily level with reliable accuracy. The algorithm was also applied to yet another cysteine-rich family, where similar accuracy was observed that ensures the application of the protocol to other families.

Comparative genomics of odorant binding proteins in Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus.

About 1 million people in the world die each year from diseases spread by mosquitoes, and understanding the mechanism of host identification by the mosquitoes through olfaction is at stake. The role of odorant binding proteins (OBPs) in the primary molecular events of olfaction in mosquitoes is becoming an important focus of biological research in this area. Here, we present a comprehensive comparative genomics study of OBPs in the three disease-transmitting mosquito species Anopheles gambiae, Aedes aegypti, and Culex quinquefasciatus starting with the identification of 110 new OBPs in these three genomes. We have characterized their genomic distribution and orthologous and phylogenetic relationships. The diversity and expansion observed with respect to the Aedes and Culex genomes suggests that the OBP gene family acquired functional diversity concurrently with functional constraints posed on these two species. Sequences with unique features have been characterized such as the "two-domain OBPs" (previously known as Atypical OBPs) and "MinusC OBPs" in mosquito genomes. The extensive comparative genomics featured in this work hence provides useful primary insights into the role of OBPs in the molecular adaptations of mosquito olfactory system and could provide more clues for the identification of potential targets for insect repellants and attractants.

Analysis on conservation of disulphide bonds and their structural features in homologous protein domain families.

BACKGROUND: Disulphide bridges are well known to play key roles in stability, folding and functions of proteins. Introduction or deletion of disulphides by site-directed mutagenesis have produced varying effects on stability and folding depending upon the protein and location of disulphide in the 3-D structure. Given the lack of complete understanding it is worthwhile to learn from an analysis of extent of conservation of disulphides in homologous proteins. We have also addressed the question of what structural interactions replaces a disulphide in a homologue in another homologue. RESULTS: Using a dataset involving 34,752 pairwise comparisons of homologous protein domains corresponding to 300 protein domain families of known 3-D structures, we provide a comprehensive analysis of extent of conservation of disulphide bridges and their structural features. We report that only 54% of all the disulphide bonds compared between the homologous pairs are conserved, even if, a small fraction of the non-conserved disulphides do include cytoplasmic proteins. Also, only about one fourth of the distinct disulphides are conserved in all the members in protein families. We note that while conservation of disulphide is common in many families, disulphide bond mutations are quite prevalent. Interestingly, we note that there is no clear relationship between sequence identity between two homologous proteins and disulphide bond conservation. Our analysis on structural features at the sites where cysteines forming disulphide in one homologue are replaced by non-Cys residues show that the elimination of a disulphide in a homologue need not always result in stabilizing interactions between equivalent residues. CONCLUSION: We observe that in the homologous proteins, disulphide bonds are conserved only to a modest extent. Very interestingly, we note that extent of conservation of disulphide in homologous proteins is unrelated to the overall sequence identity between homologues. The non-conserved disulphides are often associated with variable structural features that were recruited to be associated with differentiation or specialisation of protein function.