Mutation bias within oncogene families is related to proliferation-specific codon usage.

It is well known that in cancer gene families some members are more frequently mutated in tumor samples than their family counterparts. A paradigmatic case of this phenomenon is KRAS from the RAS family. Different explanations have been proposed ranging from differential interaction with other proteins to preferential expression or localization. Interestingly, it has been described that despite the high amino acid identity between RAS family members, KRAS employs an intriguing differential codon usage. Here, we found that this phenomenon is not exclusive to the RAS family. Indeed, in the RAS family and other oncogene families with two or three members, the most prevalently mutated gene in tumor samples employs a differential codon usage that is characteristic of genes involved in proliferation. Prompted by these observations, we chose the RAS family to experimentally demonstrate that the translation efficiency of oncogenes that are preferentially mutated in tumor samples is increased in proliferative cells compared to quiescent cells. These results were further validated by assessing the translation efficiency of KRAS in cell lines that differ in their tRNA expression profile. These differences are related to the cell division rate of the studied cells and thus suggest an important role in context-specific oncogene expression regulation. Altogether, our study demonstrates that dynamic translation programs contribute to shaping the expression profiles of oncogenes. Therefore, we propose this codon bias as a regulatory layer to control cell context-specific expression and explain the differential prevalence of mutations in certain members of oncogene families.

Translational efficiency across healthy and tumor tissues is proliferation-related.

Different tissues express genes with particular codon usage and anticodon tRNA repertoires. However, the codon-anticodon co-adaptation in humans is not completely understood, nor is its effect on tissue-specific protein levels. Here, we first validated the accuracy of small RNA-seq for tRNA quantification across five human cell lines. We then analyzed the tRNA abundance of more than 8,000 tumor samples from TCGA, together with their paired mRNA-seq and proteomics data, to determine the Supply-to-Demand Adaptation. We thereby elucidate that the dynamic adaptation of the tRNA pool is largely related to the proliferative state across tissues. The distribution of such tRNA pools over the whole cellular translatome affects the subsequent translational efficiency, which functionally determines a condition-specific expression program both in healthy and tumor states. Furthermore, the aberrant translational efficiency of some codons in cancer, exemplified by ProCCA and GlyGGT, is associated with poor patient survival. The regulation of these tRNA profiles is partly explained by the tRNA gene copy numbers and their promoter DNA methylation.

An intermediate along the recovery stroke of myosin VI revealed by X-ray crystallography and molecular dynamics.

Myosins form a class of actin-based, ATPase motor proteins that mediate important cellular functions such as cargo transport and cell motility. Their functional cycle involves two large-scale swings of the lever arm: the force-generating powerstroke, which takes place on actin, and the recovery stroke during which the lever arm is reprimed into an armed configuration. Previous analyses of the prerecovery (postrigor) and postrecovery (prepowerstroke) states predicted that closure of switch II in the ATP binding site precedes the movement of the converter and the lever arm. Here, we report on a crystal structure of myosin VI, called pretransition state (PTS), which was solved at 2.2 Å resolution. Structural analysis and all-atom molecular dynamics simulations are consistent with PTS being an intermediate along the recovery stroke, where the Relay/SH1 elements adopt a postrecovery conformation, and switch II remains open. In this state, the converter appears to be largely uncoupled from the motor domain and explores an ensemble of partially reprimed configurations through extensive, reversible fluctuations. Moreover, we found that the free energy cost of hydrogen-bonding switch II to ATP is lowered by more than 10 kcal/mol compared with the prerecovery state. These results support the conclusion that closing of switch II does not initiate the recovery stroke transition in myosin VI. Rather, they suggest a mechanism in which lever arm repriming would be mostly driven by thermal fluctuations and eventually stabilized by the switch II interaction with the nucleotide in a ratchet-like fashion.

Systems level expression correlation of Ras GTPase regulators.

BACKGROUND: Proteins of the ubiquitously expressed core proteome are quantitatively correlated across multiple eukaryotic species. In addition, it was found that many protein paralogues exhibit expression anticorrelation, suggesting that the total level of protein with a given functionality must be kept constant. METHODS: We performed Spearman's rank correlation analyses of gene expression levels for the RAS GTPase subfamily and their regulatory GEF and GAP proteins across tissues and across individuals for each tissue. A large set of published data for normal tissues from a wide range of species, human cancer tissues and human cell lines was analysed. RESULTS: We show that although the multidomain regulatory proteins of Ras GTPases exhibit considerable tissue and individual gene expression variability, their total amounts are balanced in normal tissues. In a given tissue, the sum of activating (GEFs) and deactivating (GAPs) domains of Ras GTPases can vary considerably, but each person has balanced GEF and GAP levels. This balance is impaired in cell lines and in cancer tissues for some individuals. CONCLUSIONS: Our results are relevant for critical considerations of knock out experiments, where functionally related homologs may compensate for the down regulation of a protein.

Using protein-per-mRNA differences among human tissues in codon optimization.

BACKGROUND: Codon usage and nucleotide composition of coding sequences have profound effects on protein expression. However, while it is recognized that different tissues have distinct tRNA profiles and codon usages in their transcriptomes, the effect of tissue-specific codon optimality on protein synthesis remains elusive. RESULTS: We leverage existing state-of-the-art transcriptomics and proteomics datasets from the GTEx project and the Human Protein Atlas to compute the protein-to-mRNA ratios of 36 human tissues. Using this as a proxy of translational efficiency, we build a machine learning model that identifies codons enriched or depleted in specific tissues. We detect two clusters of tissues with an opposite pattern of codon preferences. We then use these identified patterns for the development of CUSTOM, a codon optimizer algorithm which suggests a synonymous codon design in order to optimize protein production in a tissue-specific manner. In human cell-line models, we provide evidence that codon optimization should take into account particularities of the translational machinery of the tissues in which the target proteins are expressed and that our approach can design genes with tissue-optimized expression profiles. CONCLUSIONS: We provide proof-of-concept evidence that codon preferences exist in tissue-specific protein synthesis and demonstrate its application to synthetic gene design. We show that CUSTOM can be of benefit in biological and biotechnological applications, such as in the design of tissue-targeted therapies and vaccines.

Genes enriched in A/T-ending codons are co-regulated and conserved across mammals.

Codon usage influences gene expression distinctly depending on the cell context. Yet, the importance of codon bias in the simultaneous turnover of specific groups of protein-coding genes remains to be investigated. Here, we find that genes enriched in A/T-ending codons are expressed more coordinately in general and across tissues and development than those enriched in G/C-ending codons. tRNA abundance measurements indicate that this coordination is linked to the expression changes of tRNA isoacceptors reading A/T-ending codons. Genes with similar codon composition are more likely to be part of the same protein complex, especially for genes with A/T-ending codons. The codon preferences of genes with A/T-ending codons are conserved among mammals and other vertebrates. We suggest that this orchestration contributes to tissue-specific and ontogenetic-specific expression, which can facilitate, for instance, timely protein complex formation.

Translational adaptation of human viruses to the tissues they infect.

Viruses need to hijack the translational machinery of the host cell for a productive infection to happen. However, given the dynamic landscape of tRNA pools among tissues, it is unclear whether different viruses infecting different tissues have adapted their codon usage toward their tropism. Here, we collect the coding sequences of 502 human-infecting viruses and determine that tropism explains changes in codon usage. Using the tRNA abundances across 23 human tissues from The Cancer Genome Atlas (TCGA), we build an in silico model of translational efficiency that validates the correspondence of the viral codon usage with the translational machinery of their tropism. For instance, we detect that severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is specifically adapted to the upper respiratory tract and alveoli. Furthermore, this correspondence is specifically defined in early viral proteins. The observed tissue-specific translational efficiency could be useful for the development of antiviral therapies and vaccines.

From research to rapid response: mass COVID-19 testing by volunteers at the Centre for Genomic Regulation.

The COVID-19 pandemic has posed and is continuously posing enormous societal and health challenges worldwide. The research community has mobilized to develop novel projects to find a cure or a vaccine, as well as to contribute to mass testing, which has been a critical measure to contain the infection in several countries. Through this article, we share our experiences and learnings as a group of volunteers at the Centre for Genomic Regulation (CRG) in Barcelona, Spain. As members of the ORFEU project, an initiative by the Government of Catalonia to achieve mass testing of people at risk and contain the epidemic in Spain, we share our motivations, challenges and the key lessons learnt, which we feel will help better prepare the global society to address similar situations in the future.

Interaction Dynamics Determine Signaling and Output Pathway Responses.

The understanding of interaction dynamics in signaling pathways can shed light on pathway architecture and provide insights into targets for intervention. Here, we explored the relevance of kinetic rate constants of a key upstream osmosensor in the yeast high-osmolarity glycerol-mitogen-activated protein kinase (HOG-MAPK) pathway to signaling output responses. We created mutant pairs of the Sln1-Ypd1 complex interface that caused major compensating changes in the association (kon) and dissociation (koff) rate constants (kinetic perturbations) but only moderate changes in the overall complex affinity (Kd). Yeast cells carrying a Sln1-Ypd1 mutant pair with moderate increases in kon and koff displayed a lower threshold of HOG pathway activation than wild-type cells. Mutants with higher kon and koff rates gave rise to higher basal signaling and gene expression but impaired osmoadaptation. Thus, the kon and koff rates of the components in the Sln1 osmosensor determine proper signaling dynamics and osmoadaptation.

The yin-yang of kinase activation and unfolding explains the peculiarity of Val600 in the activation segment of BRAF.

Many driver mutations in cancer are specific in that they occur at significantly higher rates than - presumably - functionally alternative mutations. For example, V600E in the BRAF hydrophobic activation segment (AS) pocket accounts for >95% of all kinase mutations. While many hypotheses tried to explain such significant mutation patterns, conclusive explanations are lacking. Here, we use experimental and in silico structure-energy statistical analyses, to elucidate why the V600E mutation, but no other mutation at this, or any other positions in BRAF's hydrophobic pocket, is predominant. We find that BRAF mutation frequencies depend on the equilibrium between the destabilization of the hydrophobic pocket, the overall folding energy, the activation of the kinase and the number of bases required to change the corresponding amino acid. Using a random forest classifier, we quantitatively dissected the parameters contributing to BRAF AS cancer frequencies. These findings can be applied to genome-wide association studies and prediction models.

Tuneable endogenous mammalian target complementation via multiplexed plasmid-based recombineering.

Understanding the quantitative functional consequences of human disease mutations requires silencing of endogenous genes and expression of mutants at close to physiological levels. Changing protein levels above or below these levels is also important for system perturbation and modelling. Fast design optimization demands flexible interchangeable cassettes for endogenous gene silencing and tuneable expression. Here, we introduce 'TEMTAC', a multigene recombineering and delivery system for simultaneous siRNA-based knockdown and regulated mutant (or other variant) expression with different dynamic ranges. We show its applicability by confirming known phenotypic effects for selected mutations for BRAF, HRAS, and SHP2.

How actin initiates the motor activity of Myosin.

Fundamental to cellular processes are directional movements driven by molecular motors. A common theme for these and other molecular machines driven by ATP is that controlled release of hydrolysis products is essential for using the chemical energy efficiently. Mechanochemical transduction by myosin motors on actin is coupled to unknown structural changes that result in the sequential release of inorganic phosphate (Pi) and MgADP. We present here a myosin structure possessing an actin-binding interface and a tunnel (back door) that creates an escape route for Pi with a minimal rotation of the myosin lever arm that drives movements. We propose that this state represents the beginning of the powerstroke on actin and that Pi translocation from the nucleotide pocket triggered by actin binding initiates myosin force generation. This elucidates how actin initiates force generation and movement and may represent a strategy common to many molecular machines.