Challenges in neoantigen-directed therapeutics

A fundamental prerequisite for the efficacy of cancer immunotherapy is the presence of functional, antigen-specific T cells within the tumor. Neoantigen-directed therapy is a promising strategy that aims at targeting the host's immune response against tumor-specific antigens, thereby eradicating cancer cells. Initial forays have been made in clinical environments utilizing vaccines and adoptive cell therapy; however, many challenges lie ahead. We provide an in-depth overview of the current state of the field with an emphasis on in silico neoantigen discovery and the clinical aspects that need to be addressed to unlock the full potential of this therapy.

Entrainment to periodic initiation and transition rates in a computational model for gene translation

Periodic oscillations play an important role in many biomedical systems. Proper functioning of biological systems that respond to periodic signals requires the ability to synchronize with the periodic excitation. For example, the sleep/wake cycle is a manifestation of an internal timing system that synchronizes to the solar day. In the terminology of systems theory, the biological system must entrain or phase-lock to the periodic excitation. Entrainment is also important in synthetic biology. For example, connecting several artificial biological systems that entrain to a common clock may lead to a well-functioning modular system. The cell-cycle is a periodic program that regulates DNA synthesis and cell division. Recent biological studies suggest that cell-cycle related genes entrain to this periodic program at the gene translation level, leading to periodically-varying protein levels of these genes. The ribosome flow model (RFM) is a deterministic model obtained via a mean-field approximation of a stochastic model from statistical physics that has been used to model numerous processes including ribosome flow along the mRNA. Here we analyze the RFM under the assumption that the initiation and/or transition rates vary periodically with a common period . We show that the ribosome distribution profile in the RFM entrains to this periodic excitation. In particular, the protein synthesis pattern converges to a unique periodic solution with period . To the best of our knowledge, this is the first proof of entrainment in a mathematical model for translation that encapsulates aspects such as initiation and termination rates, ribosomal movement and interactions, and non-homogeneous elongation speeds along the mRNA. Our results support the conjecture that periodic oscillations in tRNA levels and other factors related to the translation process can induce periodic oscillations in protein levels, and may suggest a new approach for re-engineering genetic systems to obtain a desired, periodic, protein synthesis rate.

AtMBP-1, an alternative translation product of LOS2, affects abscisic acid responses and is modulated by the E3 ubiquitin ligase AtSAP5

The LOS2 gene in Arabidopsis encodes an enolase with 72% amino acid sequence identity with human ENO1. In mammalian cells, the α-enolase (ENO1) gene encodes both a 48 kDa glycolytic enzyme and a 37 kDa transcriptional suppressor protein that are targeted to different cellular compartments. The tumor suppressor c-myc binding protein (MBP-1), which is alternatively translated from the second start codon of ENO1 transcripts, is preferentially localized in nuclei while α-enolase is found in the cytoplasm. We report here that an Arabidopsis MBP-1-like protein (AtMBP-1) is alternatively translated from full-length LOS2 transcripts using a second start codon. Like mammalian MBP-1, this truncated form of LOS2 has little, if any, enolase activity, indicating that an intact N-terminal region of LOS2 is critical for catalysis. AtMBP-1 has a short half-life in vivo and is stabilized by the proteasome inhibitor MG132, indicating that it is degraded via the ubiquitin-dependent proteasome pathway. Arabidopsis plants that over-express AtMBP-1 are hypersensitive to abscisic acid (ABA) during seed germination and show defects in vegetative growth and lateral stem development. AtMBP-1 interacts directly with the E3 ubiquitin ligase AtSAP5 and co-expression of these proteins resulted in destabilization of AtMBP-1 in vivo and abolished the developmental defects associated with AtMBP-1 over-expression. Thus, AtMBP-1 is as a bona fide alternative translation product of LOS2. Accumulation of this factor is limited by ubiquitin-dependent destabilization, apparently mediated by AtSAP5.

Extensive transcript diversity and novel upstream open reading frame regulation in yeast

To understand the diversity of transcripts in yeast (Saccharomyces cerevisiae) we analyzed the transcriptional landscapes for cells grown under 18 different environmental conditions. Each sample was analyzed using RNA-sequencing, and a total of 670,446,084 uniquely mapped reads and 377,263 poly-adenylated end tags were produced. Consistent with previous studies, we find that the majority of yeast genes are expressed under one or more different conditions. By directly comparing the 5′ and 3′ ends of the transcribed regions, we find extensive differences in transcript ends across many conditions, especially those of stationary phase, growth in grape juice, and salt stimulation, suggesting differential choice of transcription start and stop sites is pervasive in yeast. Relative to the exponential growth condition (i.e., YPAD), transcripts differing at the 5′ ends and 3′ ends are predicted to differ in their annotated start codon in 21 genes and their annotated stop codon in 63 genes. Many (431) upstream open reading frames (uORFs) are found in alternate 5′ ends and are significantly enriched in transcripts produced during the salt response. Mutational analysis of five genes with uORFs revealed that two sets of uORFs increase the expression of a reporter construct, indicating a role in activation which had not been reported previously, whereas two other uORFs decreased expression. In addition, RNA binding protein motifs are statistically enriched for alternate ends under many conditions. Overall, these results demonstrate enormous diversity of transcript ends, and that this heterogeneity is regulated under different environmental conditions. Moreover, transcript end diversity has important biological implications for the regulation of gene expression. In addition, our data also serve as a valuable resource for the scientific community.

On the steady-state distribution in the homogeneous ribosome flow model

A central biological process in all living organisms is gene translation. Developing a deeper understanding of this complex process may have ramifications to almost every biomedical discipline. Reuveni et al. recently proposed a new computational model of gene translation called the Ribosome Flow Model (RFM). In this paper, we consider a particular case of this model, called the Homogeneous Ribosome Flow Model (HRFM). From a biological viewpoint, this corresponds to the case where the transition rates of all the coding sequence codons are identical. This regime has been suggested recently based on experiments in mouse embryonic cells. We consider the steady-state distribution of the HRFM. We provide formulas that relate the different parameters of the model in steady state. We prove the following properties: 1) the ribosomal density profile is monotonically decreasing along the coding sequence; 2) the ribosomal density at each codon monotonically increases with the initiation rate; and 3) for a constant initiation rate, the translation rate monotonically decreases with the length of the coding sequence. In addition, we analyze the translation rate of the HRFM at the limit of very high and very low initiation rate, and provide explicit formulas for the translation rate in these two cases. We discuss the relationship between these theoretical results and biological findings on the translation process.

Stability analysis of the ribosome flow model

Gene translation is a central process in all living organisms. Developing a better understanding of this complex process may have ramifications to almost every biomedical discipline. Recently, Reuveni et al. proposed a new computational model of this process called the ribosome flow model (RFM). In this study, we show that the dynamical behavior of the RFM is relatively simple. There exists a unique equilibrium point e and every trajectory converges to e. Furthermore, convergence is monotone in the sense that the distance to e can never increase. This qualitative behavior is maintained for any feasible set of parameter values, suggesting that the RFM is highly robust. Our analysis is based on a contraction principle and the theory of monotone dynamical systems. These analysis tools may prove useful in studying other properties of the RFM as well as additional intracellular biological processes.

Purification and characterization of cytoplasmic NADP+-isocitrate dehydrogenase, and amplification of the NADP+-IDH gene from the wing-dimorphic sand field cricket, Gryllus firmus

Cytoplasmic NADP(+)-isocitrate dehydrogenase (NADP(+)-IDH) has been purified and characterized, and its gene sequenced in many animal, plant, and yeast species. However, much less information is available on this enzyme-gene in insects. As a first step in investigating the biochemical and molecular mechanisms by which NADP(+)-IDH contributes to adaptations for flight vs. reproduction in insects, the enzyme was purified to homogeneity in the wing-dimorphic cricket, Gryllus firmus, characterized, and its corresponding gene sequenced. Using a combination of polyethylene glycol precipitation, Cibacron-Blue affinity chromatography, and hydrophobic interaction chromatography the enzyme was purified 291-fold (7% yield; specific activity = 15.8 µmol NADPH/min/mg protein). The purified enzyme exhibited a single band on SDS PAGE (46.3 kD), but consisted of two N-terminal amino acid sequences that differed in the first two amino acids. Purified enzyme exhibited standard Michaelis-Menten kinetics at pH 8.0 and 28° C (K(M(NADP+)) = 2.3 ± 0.4 µM; K(M(Na+-Isocitrate)) = 14.7 + 1.8 µM). Subunit molecular mass and K(M)S were similar to published values for NADP(+)-IDHs from a variety of vertebrate and two insect species. PCR amplification of an internal sequence using genomic DNA followed by 3' and 5' RACE yielded a nucleotide sequence of the mature protein and translated amino-acid sequences that exhibited high similarity (40-50% and 70-80%, respectively) to sequences from insect and vertebrate NADP(+)-IDHs. Two potential ATG start codons were identified. Both Nterminal amino-acid sequences matched the nucleotide sequence, consistent with both enzyme forms being transcribed from the same gene, although these variants could also be encoded by different genes. Bioinformatic analyses and differential centrifugation indicated that the majority, if not all, of the enzyme is cytoplasmic. The enzyme exhibited high specific activity in fat body, head and gut, and a single band on native PAGE.

Determinants of initiation codon selection during translation in mammalian cells

Factors affecting translation of mRNA contribute to the complexity of eukaryotic proteomes. In some cases, translation of a particular mRNA can generate multiple proteins. However, the factors that determine whether ribosomes initiate translation from the first AUG codon in the transcript, from a downstream codon, or from multiple sites are not completely understood. Various mRNA properties, including AUG codon-accessibility and 5′ leader length have been proposed as potential determinants that affect where ribosomes initiate translation. To explore this issue, we performed studies using synthetic mRNAs with two in-frame AUG codons−both in excellent context. Open reading frames initiating at AUG1 and AUG2 encode large and small isoforms of a reporter protein, respectively. Translation of such an mRNA in COS-7 cells was shown to be 5′ cap-dependent and to occur efficiently from both AUG codons. AUG codon-accessibility was modified by using two different elements: an antisense locked nucleic acid oligonucleotide and an exon-junction complex. When either element was used to mask AUG1, the ratio of the proteins synthesized changed, favoring the smaller (AUG2-initiated) protein. In addition, we observed that increased leader length by itself changed the ratio of the proteins and favored initiation at AUG1. These observations demonstrate that initiation codon selection is affected by various factors, including AUG codon-accessibility and 5′ leader length, and is not necessarily determined by the order of AUG codons (5′→3′). The modulation of AUG codon accessibility may provide a powerful means of translation regulation in eukaryotic cells.

Alternative translation start sites and hidden coding potential of eukaryotic mRNAs

It is widely suggested that a eukaryotic mRNA typically contains one translation start site and encodes a single functional protein product. However, according to current points of view on translation initiation mechanisms, eukaryotic ribosomes can recognize several alternative translation start sites and the number of experimentally verified examples of alternative translation is growing rapidly. Also, the frequent occurrence of alternative translation events and their functional significance are supported by the results of computational evaluations. The functional role of alternative translation and its contribution to eukaryotic proteome complexity are discussed.

uORFs, reinitiation and alternative translation start sites in human mRNAs

It is known that eukaryotic ribosomes are able to translate small ORFs and reinitiate translation at downstream start codons. However, this mechanism is widely considered to be inefficient and it is not commonly taken into account. We compiled a sample of human mRNAs containing small upstream ORFs overlapping with annotated protein coding sequences. Statistical analysis supported the hypothesis on reinitiation of translation at downstream AUG codons and functional significance of potential alternative ORFs. It may be assumed that some 5'UTR-located upstream ORFs can deliver ribosomes to alternative translation starts, and they should be taken into consideration in the prediction of human mRNA coding potential.