RESUMEN
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.
RESUMEN
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.
Asunto(s)
Neoplasias , Empalme del ARN , Empalme del ARN/genética , Motivos de Unión al ARN , Sitios de Empalme de ARN , Homeostasis , Factores de Empalme de ARN/genética , Neoplasias/genéticaRESUMEN
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.
Asunto(s)
Mieloma Múltiple , Degradación de ARNm Mediada por Codón sin Sentido , Animales , Humanos , Ratones , Línea Celular , ADN/metabolismo , Mamíferos/genética , Mamíferos/metabolismo , Mieloma Múltiple/tratamiento farmacológico , Mieloma Múltiple/genética , Degradación de ARNm Mediada por Codón sin Sentido/genética , Proteínas Quinasas/metabolismo , Proteínas Serina-Treonina Quinasas/genética , Proteínas Serina-Treonina Quinasas/metabolismoRESUMEN
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.
Asunto(s)
Inhibidores de Histona Desacetilasas , Linfoma de Células B Grandes Difuso , Inhibidores de Histona Desacetilasas/farmacología , Inhibidores de Histona Desacetilasas/uso terapéutico , Humanos , Factor de Transcripción Ikaros/genética , Factor de Transcripción Ikaros/metabolismo , Factores Inmunológicos/uso terapéutico , Lenalidomida/farmacología , Linfoma de Células B Grandes Difuso/tratamiento farmacológico , Linfoma de Células B Grandes Difuso/genética , Linfoma de Células B Grandes Difuso/metabolismo , Ubiquitina-Proteína Ligasas/genéticaRESUMEN
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).
Asunto(s)
Esclerosis Amiotrófica Lateral/metabolismo , ARN Nuclear Pequeño/metabolismo , Proteína FUS de Unión a ARN/genética , Proteína FUS de Unión a ARN/metabolismo , Ribonucleoproteína Nuclear Pequeña U1/metabolismo , Esclerosis Amiotrófica Lateral/genética , Animales , Línea Celular , Predisposición Genética a la Enfermedad/genética , Humanos , Ratones , Ratones Noqueados , Modelos Moleculares , Neuronas Motoras/metabolismo , Mutación , Dominios y Motivos de Interacción de Proteínas , ARN Nuclear Pequeño/química , Proteína FUS de Unión a ARN/química , Ribonucleoproteína Nuclear Pequeña U1/químicaRESUMEN
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.
Asunto(s)
Proteína FUS de Unión a ARN/química , ARN/química , Sitios de Unión , Células HeLa , Humanos , Modelos Moleculares , Motivos de Nucleótidos , Unión Proteica , Dominios y Motivos de Interacción de Proteínas , ARN/genética , ARN/metabolismo , Motivo de Reconocimiento de ARN , Empalme del ARN , Estabilidad del ARN , Proteína FUS de Unión a ARN/genética , Proteína FUS de Unión a ARN/metabolismo , Relación Estructura-Actividad , Dedos de ZincRESUMEN
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.
Asunto(s)
Linfocitos B/metabolismo , Expresión Génica/efectos de los fármacos , Enfermedad Granulomatosa Crónica/sangre , Enfermedad Granulomatosa Crónica/genética , Interferón gamma/farmacología , Receptores de Interferón/deficiencia , Linfocitos B/virología , Línea Celular , Exones/genética , Enfermedad Granulomatosa Crónica/patología , Herpesvirus Humano 4 , Humanos , Empalme del ARN/genética , ARN Mensajero/genética , RNA-Seq , Transducción de Señal/efectos de los fármacos , Triptófano-ARNt Ligasa/genética , Receptor de Interferón gammaRESUMEN
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.
Asunto(s)
Conducta Animal/fisiología , Encéfalo/metabolismo , Cíclidos/metabolismo , Conducta Social , Animales , Proteína 1 de la Respuesta de Crecimiento Precoz/metabolismo , Proteínas de Peces/metabolismo , Regulación de la Expresión Génica , ReproducciónRESUMEN
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.
Asunto(s)
Proteínas Portadoras/metabolismo , Degradación de ARNm Mediada por Codón sin Sentido/genética , ARN Helicasas/metabolismo , Receptores Citoplasmáticos y Nucleares/metabolismo , Transactivadores/metabolismo , Proteínas Portadoras/genética , Línea Celular Tumoral , Codón sin Sentido/genética , Regulación de la Expresión Génica/genética , Células HeLa , Humanos , Unión Proteica/genética , ARN Helicasas/genética , Interferencia de ARN , ARN Mensajero/genética , ARN Interferente Pequeño/genética , Receptores Citoplasmáticos y Nucleares/genética , Transactivadores/genética , Técnicas del Sistema de Dos HíbridosRESUMEN
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.
Asunto(s)
Sistemas CRISPR-Cas , Repeticiones Palindrómicas Cortas Agrupadas y Regularmente Espaciadas , Técnicas de Inactivación de Genes/métodos , Células HEK293 , Células HeLa , HumanosRESUMEN
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.
Asunto(s)
Proteínas Argonautas/genética , Proteínas de Drosophila/genética , Expresión Génica , Complejo Proteico Nuclear de Unión a la Caperuza/genética , Factores de Iniciación de Péptidos/genética , ARN Interferente Pequeño/genética , Animales , Animales Modificados Genéticamente , Proteínas Argonautas/metabolismo , Western Blotting , Núcleo Celular/genética , Núcleo Celular/metabolismo , Elementos Transponibles de ADN/genética , Proteínas de Drosophila/metabolismo , Drosophila melanogaster/genética , Drosophila melanogaster/metabolismo , Femenino , Hibridación in Situ , Masculino , Microscopía Confocal , Ovario/crecimiento & desarrollo , Ovario/metabolismo , Factores de Iniciación de Péptidos/metabolismo , Interferencia de ARN , Reacción en Cadena de la Polimerasa de Transcriptasa InversaRESUMEN
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.
Asunto(s)
Proteínas Portadoras/metabolismo , Degradación de ARNm Mediada por Codón sin Sentido , ARN Mensajero/metabolismo , Telomerasa/metabolismo , Transactivadores/metabolismo , Transcriptoma , Composición de Base , Proteínas Portadoras/antagonistas & inhibidores , Proteínas Portadoras/genética , Codón sin Sentido , Expresión Génica , Células HeLa , Humanos , Intrones , MicroARNs/química , MicroARNs/metabolismo , Unión Proteica , ARN Helicasas , ARN Largo no Codificante/química , ARN Largo no Codificante/metabolismo , ARN Mensajero/química , ARN Interferente Pequeño/genética , ARN Interferente Pequeño/metabolismo , ARN Nucleolar Pequeño/química , ARN Nucleolar Pequeño/metabolismo , Telomerasa/antagonistas & inhibidores , Telomerasa/genética , Transactivadores/antagonistas & inhibidores , Transactivadores/genéticaRESUMEN
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.
Asunto(s)
Proteínas de Drosophila/metabolismo , ARN Mensajero/genética , ARN Mensajero/metabolismo , Transcriptoma , Animales , Secuencia de Bases , Sitios de Unión , Biología Computacional , Drosophila melanogaster , Hibridación in Situ , Conformación de Ácido Nucleico , Transporte de Proteínas , Transporte de ARN , ARN Mensajero/química , Proteínas de Unión al ARN/metabolismoRESUMEN
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.
Asunto(s)
Esclerosis Amiotrófica Lateral/patología , Intrones , Proteínas Mutantes/genética , Proteínas Mutantes/metabolismo , Empalme del ARN , Proteína FUS de Unión a ARN/genética , Proteína FUS de Unión a ARN/metabolismo , Línea Celular , HumanosRESUMEN
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.
Asunto(s)
Drosophila melanogaster , Larva/fisiología , Papilas Gustativas/fisiología , Gusto/fisiología , Animales , Conducta de Elección , Proteínas de Drosophila/genética , Técnica del Anticuerpo Fluorescente , Estimulación Física , Compuestos de Amonio Cuaternario , Receptores de Superficie Celular/genética , Células Receptoras Sensoriales/fisiología , SacarosaRESUMEN
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.