RESUMO
Nonsense-mediated mRNA decay (NMD) is a eukaryotic RNA decay pathway with roles in cellular stress responses, differentiation, and viral defense. It functions in both quality control and post-transcriptional regulation of gene expression. NMD has also emerged as a modulator of cancer progression, although available evidence supports both a tumor suppressor and a pro-tumorigenic role, depending on the model. To further investigate the role of NMD in cancer, we knocked out the NMD factor SMG7 in the HT1080 human fibrosarcoma cell line, resulting in suppression of NMD function. We then compared the oncogenic properties of the parental cell line, the SMG7-knockout, and a rescue cell line in which we re-introduced both isoforms of SMG7. We also tested the effect of a drug inhibiting the NMD factor SMG1 to distinguish NMD-dependent effects from putative NMD-independent functions of SMG7. Using cell-based assays and a mouse xenograft tumor model, we showed that suppression of NMD function severely compromises the oncogenic phenotype. Molecular pathway analysis revealed that NMD suppression strongly reduces matrix metalloprotease 9 (MMP9) expression and that MMP9 re-expression partially rescues the oncogenic phenotype. Since MMP9 promotes cancer cell migration and invasion, metastasis and angiogenesis, its downregulation may contribute to the reduced tumorigenicity of NMD-suppressed cells. Collectively, our results highlight the potential value of NMD inhibition as a therapeutic approach.
RESUMO
Pharmacologic depletion of RNA-binding motif 39 (RBM39) using aryl sulfonamides represents a promising anti-cancer therapy but requires high levels of the adaptor protein DCAF15. Consequently, novel approaches to deplete RBM39 in an DCAF15-independent manner are required. Here, we uncover that RBM39 autoregulates via the inclusion of a poison exon into its own pre-mRNA and identify the cis-acting elements that govern this regulation. We also determine the NMR solution structures of RBM39's tandem RNA recognition motifs (RRM1 and RRM2) bound to their respective RNA targets, revealing how RRM1 recognises RNA stem loops whereas RRM2 binds specifically to single-stranded N(G/U)NUUUG. Our results support a model where RRM2 selects the 3'-splice site of a poison exon and the RRM3 and RS domain stabilise the U2 snRNP at the branchpoint. Our work provides molecular insights into RBM39-dependent 3'-splice site selection and constitutes a solid basis to design alternative anti-cancer therapies.
Assuntos
Neoplasias , Splicing de RNA , Splicing de RNA/genética , Motivos de Ligação ao RNA , Sítios de Splice de RNA , Homeostase , Fatores de Processamento de RNA/genética , Neoplasias/genéticaRESUMO
Early data suggested that CC-115, a clinical molecule, already known to inhibit the mammalian target of rapamycin kinase (TORK) and DNA-dependent protein kinase (DNA-PK) may have additional targets beyond TORK and DNA-PK. Therefore, we aimed to identify such target(s) and investigate a potential therapeutic applicability. Functional profiling of 141 cancer cell lines revealed inhibition of kinase suppressor of morphogenesis in genitalia 1 (SMG1), a key regulator of the RNA degradation mechanism nonsense-mediated mRNA decay (NMD), as an additional target of CC-115. CC-115 treatment showed a dose-dependent increase of SMG1-mediated NMD transcripts. A subset of cell lines, including multiple myeloma (MM) cell lines sensitive to the endoplasmic reticulum stress-inducing compound thapsigargin, were highly susceptible to SMG1 inhibition. CC-115 caused the induction of UPR transcripts and cell death by mitochondrial apoptosis, requiring the presence of BAX/BAK and caspase activity. Superior antitumor activity of CC-115 over TORK inhibitors in primary human MM cells and three xenograft mouse models appeared to be via inhibition of SMG1. Our data support further development of SMG1 inhibitors as possible therapeutics in MM.
Assuntos
Mieloma Múltiplo , Degradação do RNAm Mediada por Códon sem Sentido , Animais , Humanos , Camundongos , Linhagem Celular , DNA/metabolismo , Mamíferos/genética , Mamíferos/metabolismo , Mieloma Múltiplo/tratamento farmacológico , Mieloma Múltiplo/genética , Degradação do RNAm Mediada por Códon sem Sentido/genética , Proteínas Quinases/metabolismo , Proteínas Serina-Treonina Quinases/genética , Proteínas Serina-Treonina Quinases/metabolismoRESUMO
PURPOSE: Cereblon (CRBN), a substrate receptor of the E3 ubiquitin ligase complex CRL4CRBN, is the target of the small molecules lenalidomide and avadomide. Upon binding of the drugs, Aiolos and Ikaros are recruited to the E3 ligase, ubiquitylated, and subsequently degraded. In diffuse large B-cell lymphoma (DLBCL) cells, Aiolos and Ikaros are direct transcriptional repressors of interferon-stimulated genes (ISG) and degradation of these substrates results in increased ISG protein levels resulting in decreased proliferation and apoptosis. Herein, we aimed to uncover the mechanism(s) Aiolos and Ikaros use to repress ISG transcription and provide a mechanistic rationale for a combination strategy to enhance cell autonomous activities of CRBN modulators (CELMoD). EXPERIMENTAL DESIGN: We conducted paired RNA sequencing with histone modification and Aiolos/Ikaros chromatin immunoprecipitation sequencing to identify genes regulated by these transcription factors and to elucidate correlations to drug sensitivity. We confirmed Aiolos/Ikaros mediated transcriptional complex formation in DLBCL patient samples including those treated with avadomide. RESULTS: In DLBCL, the repression of ISG transcription is accomplished in part through recruitment of large transcriptional complexes such as the nucleosome remodeling and deacetylase, which modify the chromatin landscape of these promoters. A rational combination approach of avadomide with a specific histone deacetylase inhibitor leads to a significant increase in ISG transcription compared with either single agent, and synergistic antiproliferative activity in DLBCL cell lines. CONCLUSIONS: Our results provide a novel role for lineage factors Aiolos and Ikaros in DLBCL as well as further insight into the mechanism(s) of Aiolos and Ikaros-mediated transcriptional repression and unique therapeutic combination strategies.
Assuntos
Inibidores de Histona Desacetilases , Linfoma Difuso de Grandes Células B , Inibidores de Histona Desacetilases/farmacologia , Inibidores de Histona Desacetilases/uso terapêutico , Humanos , Fator de Transcrição Ikaros/genética , Fator de Transcrição Ikaros/metabolismo , Fatores Imunológicos/uso terapêutico , Lenalidomida/farmacologia , Linfoma Difuso de Grandes Células B/tratamento farmacológico , Linfoma Difuso de Grandes Células B/genética , Linfoma Difuso de Grandes Células B/metabolismo , Ubiquitina-Proteína Ligases/genéticaRESUMO
Mutations in the RNA-binding protein Fused in Sarcoma (FUS) cause early-onset amyotrophic lateral sclerosis (ALS). However, a detailed understanding of central RNA targets of FUS and their implications for disease remain elusive. Here, we use a unique blend of crosslinking and immunoprecipitation (CLIP) and NMR spectroscopy to identify and characterise physiological and pathological RNA targets of FUS. We find that U1 snRNA is the primary RNA target of FUS via its interaction with stem-loop 3 and provide atomic details of this RNA-mediated mode of interaction with the U1 snRNP. Furthermore, we show that ALS-associated FUS aberrantly contacts U1 snRNA at the Sm site with its zinc finger and traps snRNP biogenesis intermediates in human and murine motor neurons. Altogether, we present molecular insights into a FUS toxic gain-of-function involving direct and aberrant RNA-binding and strengthen the link between two motor neuron diseases, ALS and spinal muscular atrophy (SMA).
Assuntos
Esclerose Lateral Amiotrófica/metabolismo , RNA Nuclear Pequeno/metabolismo , Proteína FUS de Ligação a RNA/genética , Proteína FUS de Ligação a RNA/metabolismo , Ribonucleoproteína Nuclear Pequena U1/metabolismo , Esclerose Lateral Amiotrófica/genética , Animais , Linhagem Celular , Predisposição Genética para Doença/genética , Humanos , Camundongos , Camundongos Knockout , Modelos Moleculares , Neurônios Motores/metabolismo , Mutação , Domínios e Motivos de Interação entre Proteínas , RNA Nuclear Pequeno/química , Proteína FUS de Ligação a RNA/química , Ribonucleoproteína Nuclear Pequena U1/químicaRESUMO
Fused in sarcoma (FUS) is an RNA binding protein involved in regulating many aspects of RNA processing and linked to several neurodegenerative diseases. Transcriptomics studies indicate that FUS binds a large variety of RNA motifs, suggesting that FUS RNA binding might be quite complex. Here, we present solution structures of FUS zinc finger (ZnF) and RNA recognition motif (RRM) domains bound to RNA. These structures show a bipartite binding mode of FUS comprising of sequence-specific recognition of a NGGU motif via the ZnF and an unusual shape recognition of a stem-loop RNA via the RRM. In addition, sequence-independent interactions via the RGG repeats significantly increase binding affinity and promote destabilization of structured RNA conformation, enabling additional binding. We further show that disruption of the RRM and ZnF domains abolishes FUS function in splicing. Altogether, our results rationalize why deciphering the RNA binding mode of FUS has been so challenging.
Assuntos
Proteína FUS de Ligação a RNA/química , RNA/química , Sítios de Ligação , Células HeLa , Humanos , Modelos Moleculares , Motivos de Nucleotídeos , Ligação Proteica , Domínios e Motivos de Interação entre Proteínas , RNA/genética , RNA/metabolismo , Motivo de Reconhecimento de RNA , Splicing de RNA , Estabilidade de RNA , Proteína FUS de Ligação a RNA/genética , Proteína FUS de Ligação a RNA/metabolismo , Relação Estrutura-Atividade , Dedos de ZincoRESUMO
Interferon-γ (IFN-γ) plays an important role in innate and adaptive immunity against intracellular infections and is used clinically for the prevention and control of infections in chronic granulomatous disease (CGD) and inborn defects in the IFN-γ/interleukin (IL)-12 axis. Using transcriptome profiling (RNA-seq), we sought to identify differentially expressed genes, transcripts and exons in Epstein-Barr virus-transformed B lymphocytes (B-EBV) cells from CGD patients, IFN-γ receptor deficiency patients, and normal controls, treated in vitro with IFN-γ for 48 hours. Our results show that IFN-γ increased the expression of a diverse array of genes related to different cellular programs. In cells from normal controls and CGD patients, IFN-γ-induced expression of genes relevant to oxidative killing, nitric oxide synthase pathway, proteasome-mediated degradation, antigen presentation, chemoattraction, and cell adhesion. IFN-γ also upregulated genes involved in diverse stages of messenger RNA (mRNA) processing including pre-mRNA splicing, as well as others implicated in the folding, transport, and assembly of proteins. In particular, differential exon expression of WARS (encoding tryptophanyl-transfer RNA synthetase, which has an essential function in protein synthesis) induced by IFN-γ in normal and CGD cells suggests that this gene may have an important contribution to the benefits of IFN-γ treatment for CGD. Upregulation of mRNA and protein processing related genes in CGD and IFNRD cells could mediate some of the effects of IFN-γ treatment. These data support the concept that IFN-γ treatment may contribute to increased immune responses against pathogens through regulation of genes important for mRNA and protein processing.
Assuntos
Linfócitos B/metabolismo , Expressão Gênica/efeitos dos fármacos , Doença Granulomatosa Crônica/sangue , Doença Granulomatosa Crônica/genética , Interferon gama/farmacologia , Receptores de Interferon/deficiência , Linfócitos B/virologia , Linhagem Celular , Éxons/genética , Doença Granulomatosa Crônica/patologia , Herpesvirus Humano 4 , Humanos , Splicing de RNA/genética , RNA Mensageiro/genética , RNA-Seq , Transdução de Sinais/efeitos dos fármacos , Triptofano-tRNA Ligase/genética , Receptor de Interferon gamaRESUMO
In highly social species, individuals frequently face opportunities to cooperate. The molecular and neural mechanisms that integrate internal and external information prior to cooperative responses are not well understood. Using expression levels of egr-1, a genomic marker of neural activity, we quantified the neural response to an alloparental-care opportunity in a cooperatively breeding fish, a component of cooperative behaviour, across brain regions and time. In this species, alloparental care and submission are considered alternative strategies to appease dominants. We therefore investigated whether brood care and defence as well as submissive displays were associated with egr-1 expression. Finally, we predicted potential targets of the egr-1 transcription factor in the cichlid genome. This target prediction suggested that egr-1 regulates the expression of transcription factors involved in nervous system development, which could be implicated in social memory formation associated with cooperation. Egr-1 expression levels differed between test and control individuals and across time. Compared to a control, individuals experiencing the cooperation opportunity expressed less egr-1 in two brain regions, the cerebellum and the telencephalon. This down-regulation was independent of their behavioural reaction, i.e. whether they cooperated or not. However, within the subset of test individuals, egr-1 expression increased as a function of the amount of submissive behaviours, but not of cooperative behaviours, in the hypothalamus and potentially the telencephalon. These regions host structures that play a role in social decision-making; suggesting that egr-1 might be a suitable proxy for neural activation due to the social interaction component of the cooperation opportunity, rather than the actual alloparental care component.
Assuntos
Comportamento Animal/fisiologia , Encéfalo/metabolismo , Ciclídeos/metabolismo , Comportamento Social , Animais , Proteína 1 de Resposta de Crescimento Precoce/metabolismo , Proteínas de Peixes/metabolismo , Regulação da Expressão Gênica , ReproduçãoRESUMO
The term "nonsense-mediated mRNA decay" (NMD) originally described the degradation of mRNAs with premature translation-termination codons (PTCs), but its meaning has recently been extended to be a translation-dependent post-transcriptional regulator of gene expression affecting 3%-10% of all mRNAs. The degradation of NMD target mRNAs involves both exonucleolytic and endonucleolytic pathways in mammalian cells. While the latter is mediated by the endonuclease SMG6, the former pathway has been reported to require a complex of SMG5-SMG7 or SMG5-PNRC2 binding to UPF1. However, the existence, dominance, and mechanistic details of these exonucleolytic pathways are divisive. Therefore, we have investigated the possible exonucleolytic modes of mRNA decay in NMD by examining the roles of UPF1, SMG5, SMG7, and PNRC2 using a combination of functional assays and interaction mapping. Confirming previous work, we detected an interaction between SMG5 and SMG7 and also a functional need for this complex in NMD. In contrast, we found no evidence for the existence of a physical or functional interaction between SMG5 and PNRC2. Instead, we show that UPF1 interacts with PNRC2 and that it triggers 5'-3' exonucleolytic decay of reporter transcripts in tethering assays. PNRC2 interacts mainly with decapping factors and its knockdown does not affect the RNA levels of NMD reporters. We conclude that PNRC2 is probably an important mRNA decapping factor but that it does not appear to be required for NMD.
Assuntos
Proteínas de Transporte/metabolismo , Degradação do RNAm Mediada por Códon sem Sentido/genética , RNA Helicases/metabolismo , Receptores Citoplasmáticos e Nucleares/metabolismo , Transativadores/metabolismo , Proteínas de Transporte/genética , Linhagem Celular Tumoral , Códon sem Sentido/genética , Regulação da Expressão Gênica/genética , Células HeLa , Humanos , Ligação Proteica/genética , RNA Helicases/genética , Interferência de RNA , RNA Mensageiro/genética , RNA Interferente Pequeno/genética , Receptores Citoplasmáticos e Nucleares/genética , Transativadores/genética , Técnicas do Sistema de Duplo-HíbridoRESUMO
CRISPR/Cas9-based genome editing offers the possibility to knock out almost any gene of interest in an affordable and simple manner. The most common strategy is the introduction of a frameshift into the open reading frame (ORF) of the target gene which truncates the coding sequence (CDS) and targets the corresponding transcript for degradation by nonsense-mediated mRNA decay (NMD). However, we show that transcripts containing premature termination codons (PTCs) are not always degraded efficiently and can generate C-terminally truncated proteins which might have residual or dominant negative functions. Therefore, we recommend an alternative approach for knocking out genes, which combines CRISPR/Cas9 with gene traps (CRISPR-Trap) and is applicable to â¼50% of all spliced human protein-coding genes and a large subset of lncRNAs. CRISPR-Trap completely prevents the expression of the ORF and avoids expression of C-terminal truncated proteins. We demonstrate the feasibility of CRISPR-Trap by utilizing it to knock out several genes in different human cell lines. Finally, we also show that this approach can be used to efficiently generate gene replacements allowing for modulation of protein levels for otherwise lethal knockouts (KOs). Thus, CRISPR-Trap offers several advantages over conventional KO approaches and allows for generation of clean CRISPR/Cas9-based KOs.
Assuntos
Sistemas CRISPR-Cas , Repetições Palindrômicas Curtas Agrupadas e Regularmente Espaçadas , Técnicas de Inativação de Genes/métodos , Células HEK293 , Células HeLa , HumanosRESUMO
Cap binding protein 80 (Cbp80) is the larger subunit of the nuclear cap-binding complex (nCBC), which is known to play important roles in nuclear mRNA processing, export, stability and quality control events. Reducing Cbp80 mRNA levels in the female germline revealed that Cbp80 is also involved in defending the germline against transposable elements. Combining such knockdown experiments with large scale sequencing of small RNAs further showed that Cbp80 is involved in the initial biogenesis of piRNAs as well as in the secondary biogenesis pathway, the ping-pong amplification cycle. We further found that Cbp80 knockdown not only led to the upregulation of transposons, but also to delocalization of Piwi, Aub and Ago3, key factors in the piRNA biosynthesis pathway. Furthermore, compared to controls, levels of Piwi and Aub were also reduced upon knock down of Cbp80. On the other hand, with the same treatment we could not detect significant changes in levels or subcellular distribution (nuage localization) of piRNA precursor transcripts. This shows that Cbp80 plays an important role in the production and localization of the protein components of the piRNA pathway and it seems to be less important for the production and export of the piRNA precursor transcripts.
Assuntos
Proteínas Argonautas/genética , Proteínas de Drosophila/genética , Expressão Gênica , Complexo Proteico Nuclear de Ligação ao Cap/genética , Fatores de Iniciação de Peptídeos/genética , RNA Interferente Pequeno/genética , Animais , Animais Geneticamente Modificados , Proteínas Argonautas/metabolismo , Western Blotting , Núcleo Celular/genética , Núcleo Celular/metabolismo , Elementos de DNA Transponíveis/genética , Proteínas de Drosophila/metabolismo , Drosophila melanogaster/genética , Drosophila melanogaster/metabolismo , Feminino , Hibridização In Situ , Masculino , Microscopia Confocal , Ovário/crescimento & desenvolvimento , Ovário/metabolismo , Fatores de Iniciação de Peptídeos/metabolismo , Interferência de RNA , Reação em Cadeia da Polimerase Via Transcriptase ReversaRESUMO
Besides degrading aberrant mRNAs that harbor a premature translation termination codon (PTC), nonsense-mediated mRNA decay (NMD) also targets many seemingly "normal" mRNAs that encode for full-length proteins. To identify a bona fide set of such endogenous NMD targets in human cells, we applied a meta-analysis approach in which we combined transcriptome profiling of knockdowns and rescues of the three NMD factors UPF1, SMG6, and SMG7. We provide evidence that this combinatorial approach identifies NMD-targeted transcripts more reliably than previous attempts that focused on inactivation of single NMD factors. Our data revealed that SMG6 and SMG7 act on essentially the same transcripts, indicating extensive redundancy between the endo- and exonucleolytic decay routes. Besides mRNAs, we also identified as NMD targets many long noncoding RNAs as well as miRNA and snoRNA host genes. The NMD target feature with the most predictive value is an intron in the 3' UTR, followed by the presence of upstream open reading frames (uORFs) and long 3' UTRs. Furthermore, the 3' UTRs of NMD-targeted transcripts tend to have an increased GC content and to be phylogenetically less conserved when compared to 3' UTRs of NMD insensitive transcripts.
Assuntos
Proteínas de Transporte/metabolismo , Degradação do RNAm Mediada por Códon sem Sentido , RNA Mensageiro/metabolismo , Telomerase/metabolismo , Transativadores/metabolismo , Transcriptoma , Composição de Bases , Proteínas de Transporte/antagonistas & inibidores , Proteínas de Transporte/genética , Códon sem Sentido , Expressão Gênica , Células HeLa , Humanos , Íntrons , MicroRNAs/química , MicroRNAs/metabolismo , Ligação Proteica , RNA Helicases , RNA Longo não Codificante/química , RNA Longo não Codificante/metabolismo , RNA Mensageiro/química , RNA Interferente Pequeno/genética , RNA Interferente Pequeno/metabolismo , RNA Nucleolar Pequeno/química , RNA Nucleolar Pequeno/metabolismo , Telomerase/antagonistas & inibidores , Telomerase/genética , Transativadores/antagonistas & inibidores , Transativadores/genéticaRESUMO
mRNA (mRNA) transport focuses the expression of encoded proteins to specific regions within cells providing them with the means to assume specific functions and even identities. BicD and the mRNA binding protein Egl interact with the microtubule motor dynein to localize mRNAs in Drosophila. Because relatively few mRNA cargos were known, we isolated and identified Egl::GFP associated mRNAs. The top candidates were validated by qPCR, in situ hybridization and genetically by showing that their localization requires BicD. In young embryos these Egl target mRNAs are preferentially localized apically, between the plasma membrane and the blastoderm nuclei, but also in the pole plasm at the posterior pole. Egl targets expressed in the ovary were mostly enriched in the oocyte and some were apically localized in follicle cells. The identification of a large group of novel mRNAs associated with BicD/Egl points to several novel developmental and physiological functions of this dynein dependent localization machinery. The verified dataset also allowed us to develop a tool that predicts conserved A'-form-like stem loops that serve as localization elements in 3'UTRs.
Assuntos
Proteínas de Drosophila/metabolismo , RNA Mensageiro/genética , RNA Mensageiro/metabolismo , Transcriptoma , Animais , Sequência de Bases , Sítios de Ligação , Biologia Computacional , Drosophila melanogaster , Hibridização In Situ , Conformação de Ácido Nucleico , Transporte Proteico , Transporte de RNA , RNA Mensageiro/química , Proteínas de Ligação a RNA/metabolismoRESUMO
Fused in sarcoma (FUS) is a ubiquitously expressed RNA-binding protein proposed to function in various RNA metabolic pathways, including transcription regulation, pre-mRNA splicing, RNA transport and microRNA processing. Mutations in the FUS gene were identified in patients with amyotrophic lateral sclerosis (ALS), but the pathomechanisms by which these mutations cause ALS are not known. Here, we show that FUS interacts with the minor spliceosome constituent U11 snRNP, binds preferentially to minor introns and directly regulates their removal. Furthermore, a FUS knockout in neuroblastoma cells strongly disturbs the splicing of minor intron-containing mRNAs, among them mRNAs required for action potential transmission and for functional spinal motor units. Moreover, an ALS-associated FUS mutant that forms cytoplasmic aggregates inhibits splicing of minor introns by trapping U11 and U12 snRNAs in these aggregates. Collectively, our findings suggest a possible pathomechanism for ALS in which mutated FUS inhibits correct splicing of minor introns in mRNAs encoding proteins required for motor neuron survival.
Assuntos
Esclerose Lateral Amiotrófica/patologia , Íntrons , Proteínas Mutantes/genética , Proteínas Mutantes/metabolismo , Splicing de RNA , Proteína FUS de Ligação a RNA/genética , Proteína FUS de Ligação a RNA/metabolismo , Linhagem Celular , HumanosRESUMO
Accurate perception of taste information is crucial for animal survival. In adult Drosophila, gustatory receptor neurons (GRNs) perceive chemical stimuli of one specific gustatory modality associated with a stereotyped behavioural response, such as aversion or attraction. We show that GRNs of Drosophila larvae employ a surprisingly different mode of gustatory information coding. Using a novel method for calcium imaging in the larval gustatory system, we identify a multimodal GRN that responds to chemicals of different taste modalities with opposing valence, such as sweet sucrose and bitter denatonium, reliant on different sensory receptors. This multimodal neuron is essential for bitter compound avoidance, and its artificial activation is sufficient to mediate aversion. However, the neuron is also essential for the integration of taste blends. Our findings support a model for taste coding in larvae, in which distinct receptor proteins mediate different responses within the same, multimodal GRN.
Assuntos
Drosophila melanogaster , Larva/fisiologia , Papilas Gustativas/fisiologia , Paladar/fisiologia , Animais , Comportamento de Escolha , Proteínas de Drosophila/genética , Imunofluorescência , Estimulação Física , Compostos de Amônio Quaternário , Receptores de Superfície Celular/genética , Células Receptoras Sensoriais/fisiologia , SacaroseRESUMO
Relationships between evolutionary rates and gene properties on a genomic, functional, pathway, or system level are being explored to unravel the principles of the evolutionary process. In particular, functional network properties have been analyzed to recognize the constraints they may impose on the evolutionary fate of genes. Here we took as a case study the core metabolic network in human erythrocytes and we analyzed the relationship between the evolutionary rates of its genes and the metabolic flux distribution throughout it. We found that metabolic flux correlates with the ratio of nonsynonymous to synonymous substitution rates. Genes encoding enzymes that carry high fluxes have been more constrained in their evolution, while purifying selection is more relaxed in genes encoding enzymes carrying low metabolic fluxes. These results demonstrate the importance of considering the dynamical functioning of gene networks when assessing the action of selection on system-level properties.