Analogue encoding of physicochemical properties of proteins in their cognate messenger RNAs

VaxOptiML: leveraging machine learning for accurate prediction of MHC-I and II epitopes for optimized cancer immunotherapy

Cancer immunotherapy hinges on accurate epitope prediction for advancing vaccine development. VaxOptiML (available at https://vaxoptiml.streamlit.app/ ) is an integrated pipeline designed to enhance epitope prediction and prioritization. This study aims to develop and deploy a robust tool for accurate prediction and prioritization of highly immunogenic and optimized MHC-I and MHC-II T-cell epitopes for cancer vaccine development and immunotherapy. Utilizing a curated dataset of experimentally validated epitopes and employing sophisticated machine learning techniques, VaxOptiML features three models: epitope prediction from target sequences, personalized HLA typing, and prioritization the predicted epitopes based on immunogenicity scores. Our rigorous data extraction, cleaning, and feature extraction processes, coupled with model building, yield exceptional accuracy, sensitivity, specificity, and F1 score, surpassing existing prediction methods. Comprehensive visual representations underscore VaxOptiML's robustness and efficacy in accelerating epitope discovery and vaccine design for cancer immunotherapy. Deployed via Streamlit for public use, VaxOptiML enhances accessibility and usability for researchers and clinicians, demonstrating significant potential in cancer immunotherapy.

Examining SRP pathway function in mRNA localization to the endoplasmic reticulum

Signal recognition particle (SRP) pathway function in protein translocation across the endoplasmic reticulum (ER) is well established; its role in RNA localization to the ER remains, however, unclear. In current models, mRNAs undergo translation- and SRP-dependent trafficking to the ER, with ER localization mediated via interactions between SRP-bound translating ribosomes and the ER-resident SRP receptor (SR), a heterodimeric complex comprising SRA, the SRP-binding subunit, and SRB, an integral membrane ER protein. To study SRP pathway function in RNA localization, SR knockout (KO) mammalian cell lines were generated and the consequences of SR KO on steady-state and dynamic mRNA localization examined. CRISPR/Cas9-mediated SRPRB KO resulted in profound destabilization of SRA. Pairing siRNA silencing of SRPRA in SRPRB KO cells yielded viable SR KO cells. Steady-state mRNA compositions and ER-localization patterns in parental and SR KO cells were determined by cell fractionation and deep sequencing. Notably, steady-state cytosol and ER mRNA compositions and partitioning patterns were largely unaltered by loss of SR expression. To examine SRP pathway function in RNA localization dynamics, the subcellular trafficking itineraries of newly exported mRNAs were determined by 4-thiouridine (4SU) pulse-labeling/4SU-seq/cell fractionation. Newly exported mRNAs were distinguished by high ER enrichment, with ER localization being SR-independent. Intriguingly, under conditions of translation initiation inhibition, the ER was the default localization site for all newly exported mRNAs. These data demonstrate that mRNA localization to the ER can be uncoupled from the SRP pathway function and reopen questions regarding the mechanism of RNA localization to the ER.

Preferred microenvironments of halogen bonds and hydrogen bonds revealed using statistics and QM/MM calculation studies

Hydrogen bonds (HBs) and halogen bonds (XBs) are two essential non-covalent interactions for molecular recognition and drug design. As proteins are heterogeneous in structure, the microenvironments of protein structures should have effects on the formation of HBs and XBs with ligands. However, there are no systematic studies reported on this effect to date. For quantitatively describing protein microenvironments, we defined the local hydrophobicities (LHs) and local dielectric constants (LDCs) in this study. With the defined parameters, we conducted an elaborate database survey on the basis of 22 011 ligand-protein structures to explore the microenvironmental preference of HBs (91 966 in total) and XBs (1436 in total). The statistics show that XBs prefer hydrophobic microenvironments compared to HBs. The polar residues like ASP are more likely to form HBs with ligands, while nonpolar residues such as PHE and MET prefer XBs. Both the LHs and LDCs (10.69 ± 4.36 for HBs; 8.86 ± 4.00 for XBs) demonstrate that XBs are prone to hydrophobic microenvironments compared with HBs with significant differences (p < 0.001), indicating that evaluating their strengths in the corresponding environments should be necessary. Quantum Mechanics-Molecular Mechanics (QM/MM) calculations reveal that in comparison with vacuum environments, the interaction energies of HBs and XBs are decreased to varying degrees given different microenvironments. In addition, the strengths of HBs are impaired more than those of XBs when the local dielectric constant's difference between the XB microenvironments and the HB microenvironments is large.

HDLBP binds ER-targeted mRNAs by multivalent interactions to promote protein synthesis of transmembrane and secreted proteins

The biological role of RNA-binding proteins in the secretory pathway is not well established. Here, we describe that human HDLBP/Vigilin directly interacts with more than 80% of ER-localized mRNAs. PAR-CLIP analysis reveals that these transcripts represent high affinity HDLBP substrates and are specifically bound in their coding sequences (CDS), in contrast to CDS/3’UTR-bound cytosolic mRNAs. HDLBP crosslinks strongly to long CU-rich motifs, which frequently reside in CDS of ER-localized mRNAs and result in high affinity multivalent interactions. In addition to HDLBP-ncRNA interactome, quantification of HDLBP-proximal proteome confirms association with components of the translational apparatus and the signal recognition particle. Absence of HDLBP results in decreased translation efficiency of HDLBP target mRNAs, impaired protein synthesis and secretion in model cell lines, as well as decreased tumor growth in a lung cancer mouse model. These results highlight a general function for HDLBP in the translation of ER-localized mRNAs and its relevance for tumor progression.

RNA binding protein HDLBP (or Vigilin) localizes in the endoplasmic reticulum (ER) membrane. Here the authors show that HDLBP contributes to translation of ER-targeted mRNAs.

The Origin of Genetic Code and Translation in the Framework of Current Concepts on the Origin of Life

The origin of genetic code and translation system is probably the central and most difficult problem in the investigations on the origin of life and one of the most complex problems in the evolutionary biology in general. There are multiple hypotheses on the emergence and development of existing genetic systems that propose the mechanisms for the origin and early evolution of genetic code, as well as for the emergence of replication and translation. Here, we discuss the most well-known of these hypotheses, although none of them provides a description of the early evolution of genetic systems without gaps and assumptions. The RNA world hypothesis is a currently prevailing scientific idea on the early evolution of biological and pre-biological structures, the main advantage of which is the assumption that RNAs as the first living systems were self-sufficient, i.e., capable of functioning as both catalysts and templates. However, this hypothesis has also significant limitations. In particular, no ribozymes with processive polymerase activity have been yet discovered or synthesized. Taking into account the mutual need of proteins and nucleic acids in each other in the current world, many authors propose the early evolution scenarios based on the co-evolution of these two classes of organic molecules. They postulate that the emergence of translation was necessary for the replication of nucleic acids, in contrast to the RNA world hypothesis, according to which the emergence of translation was preceded by the era of self-replicating RNAs. Although such scenarios are less parsimonious from the evolutionary point of view, since they require simultaneous emergence and evolution of two classes of organic molecules, as well as the emergence of synchronized replication and translation, their major advantage is that they explain the development of processive and much more accurate protein-dependent replication.

The Physics-Biology continuum challenges darwinism: Evolution is directed by the homeostasis-dependent bidirectional relation between genome and phenotype

The physics-biology continuum relies on the fact that life emerged from prebiotic molecules. Here, I argue that life emerged from the coupling between nucleic acid and protein synthesis during which proteins (or proto-phenotypes) maintained the physicochemical parameter equilibria (or proto-homeostasis) in the proximity of their encoding nucleic acids (or proto-genomes). This protected the proto-genome physicochemical integrity (i.e., atomic composition) from environmental physicochemical constraints, and therefore increased the probability of reproducing the proto-genome without variation. From there, genomes evolved depending on the biological activities they generated in response to environmental fluctuations. Thus, a genome maintaining homeostasis (i.e., internal physicochemical parameter equilibria), despite and in response to environmental fluctuations, maintains its physicochemical integrity and has therefore a higher probability to be reproduced without variation. Consequently, descendants have a higher probability to share the same phenotype than their parents. Otherwise, the genome is modified during replication as a consequence of the imbalance of the internal physicochemical parameters it generates, until new mutation-deriving biological activities maintain homeostasis in offspring. In summary, evolution depends on feedforward and feedback loops between genome and phenotype, as the internal physicochemical conditions that a genome generates ─ through its derived phenotype in response to environmental fluctuations ─ in turn either guarantee its stability or direct its variation. Evolution may not be explained by the Darwinism-derived, unidirectional principle (random mutations-phenotypes-natural selection) but rather by the bidirectional relationship between genome and phenotype, in which the phenotype in interaction with the environment directs the evolution of the genome it derives from.

Driving forces in the origins of life

What were the physico-chemical forces that drove the origins of life? We discuss four major prebiotic 'discoveries': persistent sampling of chemical reaction space; sequence-encodable foldable catalysts; assembly of functional pathways; and encapsulation and heritability. We describe how a 'proteins-first' world gives plausible mechanisms. We note the importance of hydrophobic and polar compositions of matter in these advances.

Localization elements and zip codes in the intracellular transport and localization of messenger RNAs in Saccharomyces cerevisiae

Intracellular trafficking and localization of mRNAs provide a mechanism of regulation of expression of genes with excellent spatial control. mRNA localization followed by localized translation appears to be a mechanism of targeted protein sorting to a specific cell-compartment, which is linked to the establishment of cell polarity, cell asymmetry, embryonic axis determination, and neuronal plasticity in metazoans. However, the complexity of the mechanism and the components of mRNA localization in higher organisms prompted the use of the unicellular organism Saccharomyces cerevisiae as a simplified model organism to study this vital process. Current knowledge indicates that a variety of mRNAs are asymmetrically and selectively localized to the tip of the bud of the daughter cells, to the vicinity of endoplasmic reticulum, mitochondria, and nucleus in this organism, which are connected to diverse cellular processes. Interestingly, specific cis-acting RNA localization elements (LEs) or RNA zip codes play a crucial role in the localization and trafficking of these localized mRNAs by providing critical binding sites for the specific RNA-binding proteins (RBPs). In this review, we present a comprehensive account of mRNA localization in S. cerevisiae, various types of localization elements influencing the mRNA localization, and the RBPs, which bind to these LEs to implement a number of vital physiological processes. Finally, we emphasize the significance of this process by highlighting their connection to several neuropathological disorders and cancers. This article is categorized under: RNA Export and Localization > RNA Localization.

Physicochemical Foundations of Life that Direct Evolution: Chance and Natural Selection are not Evolutionary Driving Forces

The current framework of evolutionary theory postulates that evolution relies on random mutations generating a diversity of phenotypes on which natural selection acts. This framework was established using a top-down approach as it originated from Darwinism, which is based on observations made of complex multicellular organisms and, then, modified to fit a DNA-centric view. In this article, it is argued that based on a bottom-up approach starting from the physicochemical properties of nucleic and amino acid polymers, we should reject the facts that (i) natural selection plays a dominant role in evolution and (ii) the probability of mutations is independent of the generated phenotype. It is shown that the adaptation of a phenotype to an environment does not correspond to organism fitness, but rather corresponds to maintaining the genome stability and integrity. In a stable environment, the phenotype maintains the stability of its originating genome and both (genome and phenotype) are reproduced identically. In an unstable environment (i.e., corresponding to variations in physicochemical parameters above a physiological range), the phenotype no longer maintains the stability of its originating genome, but instead influences its variations. Indeed, environment- and cellular-dependent physicochemical parameters define the probability of mutations in terms of frequency, nature, and location in a genome. Evolution is non-deterministic because it relies on probabilistic physicochemical rules, and evolution is driven by a bidirectional interplay between genome and phenotype in which the phenotype ensures the stability of its originating genome in a cellular and environmental physicochemical parameter-depending manner.

A secretion-enhancing cis regulatory targeting element (SECReTE) involved in mRNA localization and protein synthesis

The localization of mRNAs encoding secreted/membrane proteins (mSMPs) to the endoplasmic reticulum (ER) likely facilitates the co-translational translocation of secreted proteins. However, studies have shown that mSMP recruitment to the ER in eukaryotes can occur in a manner that is independent of the ribosome, translational control, and the signal recognition particle, although the mechanism remains largely unknown. Here, we identify a cis-acting RNA sequence motif that enhances mSMP localization to the ER and appears to increase mRNA stability, and both the synthesis and secretion of secretome proteins. Termed SECReTE, for secretion-enhancing cis regulatory targeting element, this motif is enriched in mRNAs encoding secretome proteins translated on the ER in eukaryotes and on the inner membrane of prokaryotes. SECReTE consists of ≥10 nucleotide triplet repeats enriched with pyrimidine (C/U) every third base (i.e. NNY, where N = any nucleotide, Y = pyrimidine) and can be present in the untranslated as well as the coding regions of the mRNA. Synonymous mutations that elevate the SECReTE count in a given mRNA (e.g. SUC2, HSP150, and CCW12) lead to an increase in protein secretion in yeast, while a reduction in count led to less secretion and physiological defects. Moreover, the addition of SECReTE to the 3'UTR of an mRNA for an exogenously expressed protein (e.g. GFP) led to its increased secretion from yeast cells. Thus, SECReTE constitutes a novel RNA motif that facilitates ER-localized mRNA translation and protein secretion.

RNA-protein interactions in an unstructured context

Despite their importance, our understanding of noncovalent RNA-protein interactions is incomplete. This especially concerns the binding between RNA and unstructured protein regions, a widespread class of such interactions. Here, we review the recent experimental and computational work on RNA-protein interactions in an unstructured context with a particular focus on how such interactions may be shaped by the intrinsic interaction affinities between individual nucleobases and protein side chains. Specifically, we articulate the claim that the universal genetic code reflects the binding specificity between nucleobases and protein side chains and that, in turn, the code may be seen as the Rosetta stone for understanding RNA-protein interactions in general.

Self-Consistent Framework Connecting Experimental Proxies of Protein Dynamics with Configurational Entropy

The recently developed NMR techniques enable estimation of protein configurational entropy change from the change in the average methyl order parameters. This experimental observable, however, does not directly measure the contribution of intramolecular couplings, protein main-chain motions, or angular dynamics. Here, we carry out a self-consistent computational analysis of the impact of these missing contributions on an extensive set of molecular dynamics simulations of different proteins undergoing binding. Specifically, we compare the configurational entropy change in protein complex formation as obtained by the maximum information spanning tree approximation (MIST), which treats the above entropy contributions directly, and the change in the average NMR methyl and NH order parameters. Our parallel implementation of MIST allows us to treat hard angular degrees of freedom as well as couplings up to full pairwise order explicitly, while still involving a high degree of sampling and tackling molecules of biologically relevant sizes. First, we demonstrate a remarkably strong linear relationship between the total configurational entropy change and the average change in both methyl and backbone-NH order parameters. Second, in contrast to canonical assumptions, we show that the main-chain and angular terms contribute significantly to the overall configurational entropy change and also scale linearly with it. Consequently, linear models starting from the average methyl order parameters are able to capture the contribution of main-chain and angular terms well. After applying the quantum-mechanical harmonic oscillator entropy formalism, we establish a similarly strong linear relationship for X-ray crystallographic B-factors. Finally, we demonstrate that the observed linear relationships remain robust against drastic undersampling and argue that they reflect an intrinsic property of compact proteins. Despite their remarkable strength, however, the above linear relationships yield estimates of configurational entropy change whose accuracy appears to be sufficient for qualitative applications only.

mRNA/protein sequence complementarity and its determinants: The impact of affinity scales

It has recently been demonstrated that the nucleobase-density profiles of mRNA coding sequences are related in a complementary manner to the nucleobase-affinity profiles of their cognate protein sequences. Based on this, it has been proposed that cognate mRNA/protein pairs may bind in a co-aligned manner, especially if unstructured. Here, we study the dependence of mRNA/protein sequence complementarity on the properties of the nucleobase/amino-acid affinity scales used. Specifically, we sample the space of randomly generated scales by employing a Monte Carlo strategy with a fitness function that depends directly on the level of complementarity. For model organisms representing all three domains of life, we show that even short searches reproducibly converge upon highly optimized scales, implying that the topology of the underlying fitness landscape is decidedly funnel-like. Furthermore, the optimized scales, generated without any consideration of the physicochemical attributes of nucleobases or amino acids, resemble closely the nucleobase/amino-acid binding affinity scales obtained from experimental structures of RNA-protein complexes. This provides support for the claim that mRNA/protein sequence complementarity may indeed be related to binding between the two. Finally, we characterize suboptimal scales and show that intermediate-to-high complementarity can be reached by substantially diverse scales, but with select amino acids contributing disproportionally. Our results expose the dependence of cognate mRNA/protein sequence complementarity on the properties of the underlying nucleobase/amino-acid affinity scales and provide quantitative constraints that any physical scales need to satisfy for the complementarity to hold.

Messenger RNAs and proteins, two essential types of biopolymers, have recently been shown to exhibit closely related, complementary physicochemical properties. Specifically, density profiles of certain groups in messenger RNA sequences directly match the affinity profiles for precisely those groups in protein sequences they encode. Based on this, it has been suggested that these molecules may interact with each other specifically and in a co-aligned fashion, especially when unstructured. Here, we explore different amino-acid scales used in the above analysis to assess which of their properties dictate the observed matching. Specifically, we define the constraints that need to be satisfied by physical scales for the complementarity to hold and show that the previously derived nucleobase/amino-acid affinity scales indeed satisfy these constraints. As a whole, our work provides a quantitative foundation for understanding the putative messenger RNA/protein complementarity with implications in different areas of RNA/protein biology including transcription, translation, splicing and viral assembly.

PHYSICO2: an UNIX based standalone procedure for computation of physicochemical, window-dependent and substitution based evolutionary properties of protein sequences along with automated block preparation tool, version 2

Automated genome sequencing procedure is enriching the sequence database very fast. To achieve a balance between the entry of sequences in the database and their analyses, efficient software is required. In this end PHYSICO2, compare to earlier PHYSICO and other public domain tools, is most efficient in that it i] extracts physicochemical, window-dependent and homologousposition-based-substitution (PWS) properties including positional and BLOCK-specific diversity and conservation, ii] provides users with optional-flexibility in setting relevant input-parameters, iii] helps users to prepare BLOCK-FASTA-file by the use of Automated Block Preparation Tool of the program, iv] performs fast, accurate and user-friendly analyses and v] redirects itemized outputs in excel format along with detailed methodology. The program package contains documentation describing application of methods. Overall the program acts as efficient PWS-analyzer and finds application in sequence-bioinformatics.

AVAILABILITY

PHYSICO2: is freely available at http://sourceforge.net/projects/physico2/ along with its documentation at https://sourceforge.net/projects/physico2/files/Documentation.pdf/download for all users.

Diversity and selectivity in mRNA translation on the endoplasmic reticulum

Pioneering electron microscopy studies defined two primary populations of ribosomes in eukaryotic cells: one freely dispersed through the cytoplasm and the other bound to the surface of the endoplasmic reticulum (ER). Subsequent investigations revealed a specialized function for each population, with secretory and integral membrane protein-encoding mRNAs translated on ER-bound ribosomes, and cytosolic protein synthesis was widely attributed to free ribosomes. Recent findings have challenged this view, and transcriptome-scale studies of mRNA distribution and translation have now demonstrated that ER-bound ribosomes also function in the translation of a large fraction of mRNAs that encode cytosolic proteins. These studies suggest a far more expansive role for the ER in transcriptome expression, where membrane and secretory protein synthesis represents one element of a multifaceted and dynamic contribution to post-transcriptional gene expression.

On the Contribution of Protein Spatial Organization to the Physicochemical Interconnection between Proteins and Their Cognate mRNAs

Early-stage evolutionary development of the universal genetic code remains a fundamental, open problem. One of the possible scenarios suggests that the code evolved in response to direct interactions between peptides and RNA oligonucleotides in the primordial environment. Recently, we have revealed a strong matching between base-binding preferences of modern protein sequences and the composition of their cognate mRNA coding sequences. These results point directly at the physicochemical foundation behind the code’s origin, but also support the possibility of direct complementary interactions between proteins and their cognate mRNAs, especially if the two are unstructured. Here, we analyze molecular-surface mapping of knowledge-based amino-acid/nucleobase interaction preferences for a set of complete, high-resolution protein structures and show that the connection between the two biopolymers could remain relevant even for structured, folded proteins. Specifically, protein surface loops are strongly enriched in residues with a high binding propensity for guanine and cytosine, while adenine- and uracil-preferring residues are uniformly distributed throughout protein structures. Moreover, compositional complementarity of cognate protein and mRNA sequences remains strong even after weighting protein sequence profiles by residue solvent exposure. Our results support the possibility that protein/mRNA sequence complementarity may also translate to cognate interactions between structured biopolymers.

Localization of mRNAs to the endoplasmic reticulum

Almost all cells use mRNA localization to establish spatial control of protein synthesis. One of the best-studied examples is the targeting and anchoring of mRNAs encoding secreted, organellar, and membrane-bound proteins to the surface of the endoplasmic reticulum (ER). In this review, we provide a historical perspective on the research that elucidated the canonical protein-mediated targeting of nascent-chain ribosome mRNA complexes to the surface of the ER. We then discuss subsequent studies which provided concrete evidence that a subpopulation of mRNAs utilize a translation-independent mechanism to localize to the surface of this organelle. This alternative mechanism operates alongside the signal recognition particle (SRP) mediated co-translational targeting pathway to promote proper mRNA localization to the ER. Recent work has uncovered trans-acting factors, such as the mRNA receptor p180, and cis-acting elements, such as transmembrane domain coding regions, that are responsible for this alternative mRNA localization process. Furthermore, some unanticipated observations have raised the possibility that this alternative pathway may be conserved from bacteria to mammalian cells.