Computational modeling of eukaryotic mRNA turnover

Cytoplasmic mRNA decay and quality control machineries in eukaryotes

mRNA degradation pathways have key regulatory roles in gene expression. The intrinsic stability of mRNAs in the cytoplasm of eukaryotic cells varies widely in a gene- and isoform-dependent manner and can be regulated by cellular cues, such as kinase signalling, to control mRNA levels and spatiotemporal dynamics of gene expression. Moreover, specialized quality control pathways exist to rid cells of non-functional mRNAs produced by errors in mRNA processing or mRNA damage that negatively impact translation. Recent advances in structural, single-molecule and genome-wide methods have provided new insights into the central machineries that carry out mRNA turnover, the mechanisms by which mRNAs are targeted for degradation and the general principles that govern mRNA stability at a global level. This improved understanding of mRNA degradation in the cytoplasm of eukaryotic cells is finding practical applications in the design of therapeutic mRNAs.

RNA degradation triggered by decapping is largely independent of initial deadenylation

RNA stability, important for eukaryotic gene expression, is thought to depend on deadenylation rates, with shortened poly(A) tails triggering decapping and 5' to 3' degradation. In contrast to this view, recent large-scale studies indicate that the most unstable mRNAs have, on average, long poly(A) tails. To clarify the role of deadenylation in mRNA decay, we first modeled mRNA poly(A) tail kinetics and mRNA stability in yeast. Independent of deadenylation rates, differences in mRNA decapping rates alone were sufficient to explain current large-scale results. To test the hypothesis that deadenylation and decapping are uncoupled, we used rapid depletion of decapping and deadenylation enzymes and measured changes in mRNA levels, poly(A) length and stability, both transcriptome-wide and with individual reporters. These experiments revealed that perturbations in poly(A) tail length did not correlate with variations in mRNA stability. Thus, while deadenylation may be critical for specific regulatory mechanisms, our results suggest that for most yeast mRNAs, it is not critical for mRNA decapping and degradation.

The role of polymers in enabling RNAi-based technology for sustainable pest management

The growing global food demand, coupled with the limitations of traditional pest control methods, has driven the search for innovative and sustainable solutions in agricultural pest management. In this review, we highlight polymeric nanocarriers for their potential to deliver double-stranded RNA (dsRNA) and control pests through the gene-silencing mechanism of RNA interference (RNAi). Polymer-dsRNA systems have shown promise in protecting dsRNA, facilitating cellular uptake, and ensuring precise release. Despite these advances, challenges such as scalability, cost-efficiency, regulatory approval, and public acceptance persist, necessitating further research to overcome these obstacles and fully unlock the potential of RNAi in sustainable agriculture.

Analysis of a detailed multi-stage model of stochastic gene expression using queueing theory and model reduction

We introduce a biologically detailed, stochastic model of gene expression describing the multiple rate-limiting steps of transcription, nuclear pre-mRNA processing, nuclear mRNA export, cytoplasmic mRNA degradation and translation of mRNA into protein. The processes in sub-cellular compartments are described by an arbitrary number of processing stages, thus accounting for a significantly finer molecular description of gene expression than conventional models such as the telegraph, two-stage and three-stage models of gene expression. We use two distinct tools, queueing theory and model reduction using the slow-scale linear-noise approximation, to derive exact or approximate analytic expressions for the moments or distributions of nuclear mRNA, cytoplasmic mRNA and protein fluctuations, as well as lower bounds for their Fano factors in steady-state conditions. We use these to study the phase diagram of the stochastic model; in particular we derive parametric conditions determining three types of transitions in the properties of mRNA fluctuations: from sub-Poissonian to super-Poissonian noise, from high noise in the nucleus to high noise in the cytoplasm, and from a monotonic increase to a monotonic decrease of the Fano factor with the number of processing stages. In contrast, protein fluctuations are always super-Poissonian and show weak dependence on the number of mRNA processing stages. Our results delineate the region of parameter space where conventional models give qualitatively incorrect results and provide insight into how the number of processing stages, e.g. the number of rate-limiting steps in initiation, splicing and mRNA degradation, shape stochastic gene expression by modulation of molecular memory.

Decoupled degradation and translation enables noise modulation by poly(A) tails

Poly(A) tails are crucial for mRNA translation and degradation, but the exact relationship between tail length and mRNA kinetics remains unclear. Here, we employ a small library of identical mRNAs that differ only in their poly(A)-tail length to examine their behavior in human embryonic kidney cells. We find that tail length strongly correlates with mRNA degradation rates but is decoupled from translation. Interestingly, an optimal tail length of ∼100 nt displays the highest translation rate, which is identical to the average endogenous tail length measured by nanopore sequencing. Furthermore, poly(A)-tail length variability-a feature of endogenous mRNAs-impacts translation efficiency but not mRNA degradation rates. Stochastic modeling combined with single-cell tracking reveals that poly(A) tails provide cells with an independent handle to tune gene expression fluctuations by decoupling mRNA degradation and translation. Together, this work contributes to the basic understanding of gene expression regulation and has potential applications in nucleic acid therapeutics.

Current insight into the role of mRNA decay pathways in fungal pathogenesis

Pathogenic fungal species can cause superficial and mucosal infections, to potentially fatal systemic or invasive infections in humans. These infections are more common in immunocompromised or critically ill patients and have a significant morbidity and fatality rate. Fungal pathogens utilize several strategies to adapt the host environment resulting in efficient and comprehensive alterations in their cellular metabolism. Fungal virulence is regulated by several factors and post-transcriptional regulation mechanisms involving mRNA molecules are one of them. Post-transcriptional controls have emerged as critical regulatory mechanisms involved in the pathogenesis of fungal species. The untranslated upstream and downstream regions of the mRNA, as well as RNA-binding proteins, regulate morphogenesis and virulence by controlling mRNA degradation and stability. The limited number of available therapeutic drugs, the emergence of multidrug resistance, and high death rates associated with systemic fungal illnesses pose a serious risk to human health. Therefore, new antifungal treatments that specifically target mRNA pathway components can decrease fungal pathogenicity and when combined increase the effectiveness of currently available antifungal drugs. This review summarizes the mRNA degradation pathways and their role in fungal pathogenesis.

Solving stochastic gene-expression models using queueing theory: A tutorial review

Stochastic models of gene expression are typically formulated using the chemical master equation, which can be solved exactly or approximately using a repertoire of analytical methods. Here, we provide a tutorial review of an alternative approach based on queueing theory that has rarely been used in the literature of gene expression. We discuss the interpretation of six types of infinite-server queues from the angle of stochastic single-cell biology and provide analytical expressions for the stationary and nonstationary distributions and/or moments of mRNA/protein numbers and bounds on the Fano factor. This approach may enable the solution of complex models that have hitherto evaded analytical solution.

Distilling dynamical knowledge from stochastic reaction networks

Stochastic reaction networks are widely used in the modeling of stochastic systems across diverse domains such as biology, chemistry, physics, and ecology. However, the comprehension of the dynamic behaviors inherent in stochastic reaction networks is a formidable undertaking, primarily due to the exponential growth in the number of possible states or trajectories as the state space dimension increases. In this study, we introduce a knowledge distillation method based on reinforcement learning principles, aimed at compressing the dynamical knowledge encoded in stochastic reaction networks into a singular neural network construct. The trained neural network possesses the capability to accurately predict the state conditional joint probability distribution that corresponds to the given query contexts, when prompted with rate parameters, initial conditions, and time values. This obviates the need to track the dynamical process, enabling the direct estimation of normalized state and trajectory probabilities, without necessitating the integration over the complete state space. By applying our method to representative examples, we have observed a high degree of accuracy in both multimodal and high-dimensional systems. Additionally, the trained neural network can serve as a foundational model for developing efficient algorithms for parameter inference and trajectory ensemble generation. These results collectively underscore the efficacy of our approach as a universal means of distilling knowledge from stochastic reaction networks. Importantly, our methodology also spotlights the potential utility in harnessing a singular, pretrained, large-scale model to encapsulate the solution space underpinning a wide spectrum of stochastic dynamical systems.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Gene expression inherently gives rise to stochastic variation ("noise") in the production of gene products. Minimizing noise is crucial for ensuring reliable cellular functions. However, noise cannot be suppressed below a certain intrinsic limit. For constitutively expressed genes, this limit is typically assumed to be Poissonian noise, wherein the variance in mRNA numbers is equal to their mean. Here, we demonstrate that several cell division genes in fission yeast exhibit mRNA variances significantly below this limit. The reduced variance can be explained by a gene expression model incorporating multiple transcription and mRNA degradation steps. Notably, in this sub-Poissonian regime, distinct from Poissonian or super-Poissonian regimes, cytoplasmic noise is effectively suppressed through a higher mRNA export rate. Our findings redefine the lower limit of eukaryotic gene expression noise and uncover molecular requirements for achieving ultralow noise, which is expected to be important for vital cellular functions.

Contributions of Ccr4 and Gcn2 to the Translational Response of C. neoformans to Host-Relevant Stressors and Integrated Stress Response Induction

In response to the host environment, the human pathogen Cryptococcus neoformans must rapidly reprogram its translatome from one which promotes growth to one which is responsive to host stress. In this study, we investigate the two events which comprise translatome reprogramming: the removal of abundant, pro-growth mRNAs from the translating pool, and the regulated entry of stress-responsive mRNAs into the translating pool. Removal of pro-growth mRNAs from the translating pool is controlled primarily by two regulatory mechanisms, repression of translation initiation via Gcn2, and decay mediated by Ccr4. We determined that translatome reprogramming in response to oxidative stress requires both Gcn2 and Ccr4, whereas the response to temperature requires only Ccr4. Additionally, we assessed ribosome collision in response to host-relevant stress and found that collided ribosomes accumulated during temperature stress but not during oxidative stress. The phosphorylation of eIF2α that occurred as a result of translational stress led us to investigate the induction of the integrated stress response (ISR). We found that eIF2α phosphorylation varied in response to the type and magnitude of stress, yet all tested conditions induced translation of the ISR transcription factor Gcn4. However, Gcn4 translation did not necessarily result in canonical Gcn4-dependent transcription. Finally, we define the ISR regulon in response to oxidative stress. In conclusion, this study begins to reveal the translational regulation in response to host-relevant stressors in an environmental fungus which is capable of adapting to the environment inside the human host. IMPORTANCE Cryptococcus neoformans is a human pathogen capable of causing devastating infections. It must rapidly adapt to changing environments as it leaves its niche in the soil and enters the human lung. Previous work has demonstrated a need to reprogram gene expression at the level of translation to promote stress adaptation. In this work, we investigate the contributions and interplay of the major mechanisms that regulate entry of new mRNAs into the pool (translation initiation) and the clearance of unneeded mRNAs from the pool (mRNA decay). One result of this reprogramming is the induction of the integrated stress response (ISR) regulon. Surprisingly, all stresses tested led to the production of the ISR transcription factor Gcn4, but not necessarily to transcription of ISR target genes. Furthermore, stresses result in differential levels of ribosome collisions, but these are not necessarily predictive of initiation repression as has been suggested in the model yeast.

The minimal intrinsic stochasticity of constitutively expressed eukaryotic genes is sub-Poissonian

Stochastic variation in gene products ("noise") is an inescapable by-product of gene expression. Noise must be minimized to allow for the reliable execution of cellular functions. However, noise cannot be suppressed beyond an intrinsic lower limit. For constitutively expressed genes, this limit is believed to be Poissonian, meaning that the variance in mRNA numbers cannot be lower than their mean. Here, we show that several cell division genes in fission yeast have mRNA variances significantly below this limit, which cannot be explained by the classical gene expression model for low-noise genes. Our analysis reveals that multiple steps in both transcription and mRNA degradation are essential to explain this sub-Poissonian variance. The sub-Poissonian regime differs qualitatively from previously characterized noise regimes, a hallmark being that cytoplasmic noise is reduced when the mRNA export rate increases. Our study re-defines the lower limit of eukaryotic gene expression noise and identifies molecular requirements for ultra-low noise which are expected to support essential cell functions.

Rapid nuclear deadenylation of mammalian messenger RNA

Poly(A) tails protect RNAs from degradation and their deadenylation rates determine RNA stability. Although poly(A) tails are generated in the nucleus, deadenylation of tails has mostly been investigated within the cytoplasm. Here, we combined long-read sequencing with metabolic labeling, splicing inhibition and cell fractionation experiments to quantify, separately, the genesis and trimming of nuclear and cytoplasmic tails in vitro and in vivo. We present evidence for genome-wide, nuclear synthesis of tails longer than 200 nt, which are rapidly shortened after transcription. Our data suggests that rapid deadenylation is a nuclear process, and that different classes of transcripts and even transcript isoforms have distinct nuclear tail lengths. For example, many long-noncoding RNAs retain long poly(A) tails. Modeling deadenylation dynamics predicts nuclear deadenylation about 10 times faster than cytoplasmic deadenylation. In summary, our data suggests that nuclear deadenylation might be a key mechanism for regulating mRNA stability, abundance, and subcellular localization.

Approximating solutions of the Chemical Master equation using neural networks

The Chemical Master Equation (CME) provides an accurate description of stochastic biochemical reaction networks in well-mixed conditions, but it cannot be solved analytically for most systems of practical interest. Although Monte Carlo methods provide a principled means to probe system dynamics, the large number of simulations typically required can render the estimation of molecule number distributions and other quantities infeasible. In this article, we aim to leverage the representational power of neural networks to approximate the solutions of the CME and propose a framework for the Neural Estimation of Stochastic Simulations for Inference and Exploration (Nessie). Our approach is based on training neural networks to learn the distributions predicted by the CME from relatively few stochastic simulations. We show on biologically relevant examples that simple neural networks with one hidden layer can capture highly complex distributions across parameter space, thereby accelerating computationally intensive tasks such as parameter exploration and inference.

Mean-field theory accurately captures the variation of copy number distributions across the mRNA life cycle

We consider a stochastic model where a gene switches between two states, an mRNA transcript is released in the active state, and subsequently it undergoes an arbitrary number of sequential unimolecular steps before being degraded. The reactions effectively describe various stages of the mRNA life cycle such as initiation, elongation, termination, splicing, export, and degradation. We construct a mean-field approach that leads to closed-form steady-state distributions for the number of transcript molecules at each stage of the mRNA life cycle. By comparison with stochastic simulations, we show that the approximation is highly accurate over all the parameter space, independent of the type of expression (constitutive or bursty) and of the shape of the distribution (unimodal, bimodal, and nearly bimodal). The theory predicts that in a population of identical cells, any bimodality is gradually washed away as the mRNA progresses through its life cycle.

Kinetic Study of 17α-Estradiol Activity in Comparison with 17β-Estradiol and 17α-Ethynylestradiol

17α-estradiol (αE2), an endogenous stereoisomer of the hormone 17β-estradiol (E2), is capable of binding to estrogen receptors (ER). We aimed to mathematically describe, using experimental data, the possible interactions between αE2 and sperm ER during the process of sperm capacitation and to develop a kinetic model. The goal was to compare the suggested kinetic model with previously published results of ER interactions with E2 and 17α-ethynylestradiol (EE2). The HPLC-MS/MS method was developed to monitor the changes of αE2 concentration during capacitation. The calculated relative concentrations Bt were used for kinetic analysis. Rate constants k and molar ratio n were optimized and used for the construction of theoretical B(t) curves. Modifications in αE2–ER interactions were discovered during comparison with models for E2 and EE2. These new interactions displayed autocatalytic formation of an unstable adduct between the hormone and the cytoplasmic receptors. αE2 accumulates between the plasma membrane lipid bilayer with increasing potential, and when the critical level is reached, αE2 penetrates through the inner layer of the plasma membrane into the cytoplasm. It then rapidly reacts with the ER and creates an unstable adduct. The revealed dynamics of αE2–ER action may contribute to understanding tissue rejuvenation and the cancer-related physiology of αE2 signaling.

Poly(A) tail dynamics: Measuring polyadenylation, deadenylation and poly(A) tail length

Transcription of mRNAs culminates in RNA cleavage and a coordinated polyadenylation event at the 3' end. In its journey to be translated, the resulting transcript is under constant regulation by cap-binding proteins, miRNAs, and RNA binding proteins, including poly(A) binding proteins (PABPs). The interplay between all these factors determines whether nuclear or cytoplasmic exoribonucleases will gain access to and remove the poly(A) tail, which is so critical to the stability and translation capacity of the mRNA. In this chapter, we present an overview of two of the key features of the mRNA life-cycle: cleavage/polyadenylation and deadenylation, and describe biochemical assays that have been generated to study the activity of each of these enzymatic reactions. Finally, we also provide protocols to investigate mRNA's poly(A) length. The importance of these assays is highlighted by the dynamic and essential role the poly(A) tail length plays in controlling gene expression.

Tob2 phosphorylation regulates global mRNA turnover to reshape transcriptome and impact cell proliferation

Tob2, an anti-proliferative protein, promotes deadenylation through recruiting Caf1 deadenylase to the mRNA poly(A) tail by simultaneously interacting with both Caf1 and poly(A)-binding protein (PABP). Previously, we found that changes in Tob2 phosphorylation can alter its PABP-binding ability and deadenylation-promoting function. However, it remained unknown regarding the relevant kinase(s). Moreover, it was unclear whether Tob2 phosphorylation modulates the transcriptome and whether the phosphorylation is linked to Tob2's anti-proliferative function. In this study, we found that c-Jun amino-terminal kinase (JNK) increases phosphorylation of Tob2 at many Ser/Thr sites in the intrinsically disordered region (IDR) that contains two separate PABP-interacting PAM2 motifs. JNK-induced phosphorylation or phosphomimetic mutations at these sites weaken the Tob2-PABP interaction. In contrast, JNK-independent phosphorylation of Tob2 at serine 254 (S254) greatly enhances Tob2 interaction with PABP and its ability to promote deadenylation. We discovered that both PAM2 motifs are required for Tob2 to display these features. Combining mass spectrometry analysis, poly(A) size-distribution profiling, transcriptome-wide mRNA turnover analyses, and cell proliferation assays, we found that the phosphomimetic mutation at S254 (S254D) enhances Tob2's association with PABP, leading to accelerated deadenylation and decay of mRNAs globally. Moreover, the Tob2-S254D mutant accelerates the decay of many transcripts coding for cell cycle related proteins and enhances anti-proliferation function. Our findings reveal a novel mechanism by which Ccr4-Not complex is recruited by Tob2 to the mRNA 3' poly(A)-PABP complex in a phosphorylation dependent manner to promote rapid deadenylation and decay across the transcriptome, eliciting transcriptome reprogramming and suppressed cell proliferation.

The Dynamics of Cytoplasmic mRNA Metabolism

For all but a few mRNAs, the dynamics of metabolism are unknown. Here, we developed an experimental and analytical framework for examining these dynamics for mRNAs from thousands of genes. mRNAs of mouse fibroblasts exit the nucleus with diverse intragenic and intergenic poly(A)-tail lengths. Once in the cytoplasm, they have a broad (1000-fold) range of deadenylation rate constants, which correspond to cytoplasmic lifetimes. Indeed, with few exceptions, degradation appears to occur primarily through deadenylation-linked mechanisms, with little contribution from either endonucleolytic cleavage or deadenylation-independent decapping. Most mRNA molecules degrade only after their tail lengths fall below 25 nt. Decay rate constants of short-tailed mRNAs vary broadly (1000-fold) and are larger for short-tailed mRNAs that have previously undergone more rapid deadenylation. This coupling helps clear rapidly deadenylated mRNAs, enabling the large range in deadenylation rate constants to impart a similarly large range in stabilities.

Translational control of gene expression via interacting feedback loops

Translation is a key step in the synthesis of proteins. Accordingly, cells have evolved an intricate array of control mechanisms to regulate this process. By constructing a multicomponent mathematical framework we uncover how translation may be controlled via interacting feedback loops. Our results reveal that this interplay gives rise to a remarkable range of protein synthesis dynamics, including oscillations, step change, and bistability. This suggests that cells may have recourse to a much richer set of control mechanisms than was previously understood.

C3P3-G1: first generation of a eukaryotic artificial cytoplasmic expression system

Most eukaryotic expression systems make use of host-cell nuclear transcriptional and post-transcriptional machineries. Here, we present the first generation of the chimeric cytoplasmic capping-prone phage polymerase (C3P3-G1) expression system developed by biological engineering, which generates capped and polyadenylated transcripts in host-cell cytoplasm by means of two components. First, an artificial single-unit chimeric enzyme made by fusing an mRNA capping enzyme and a DNA-dependent RNA polymerase. Second, specific DNA templates designed to operate with the C3P3-G1 enzyme, which encode for the transcripts and their artificial polyadenylation. This system, which can potentially be adapted to any in cellulo or in vivo eukaryotic expression applications, was optimized for transient expression in mammalian cells. C3P3-G1 shows promising results for protein production in Chinese Hamster Ovary (CHO-K1) cells. This work also provides avenues for enhancing the performances for next generation C3P3 systems.

Inverted translational control of eukaryotic gene expression by ribosome collisions

The canonical model of eukaryotic translation posits that efficient translation initiation increases protein expression and mRNA stability. Contrary to this model, we find that increasing initiation rate can decrease both protein expression and stability of certain mRNAs in the budding yeast Saccharomyces cerevisiae. These mRNAs encode a stretch of polybasic residues that cause ribosome stalling. Our computational modeling predicts that the observed decrease in gene expression at high initiation rates occurs when ribosome collisions at stalls stimulate abortive termination of the leading ribosome or cause endonucleolytic mRNA cleavage. Consistent with this prediction, the collision-associated quality-control factors Asc1 and Hel2 (orthologs of human RACK1 and ZNF598, respectively) decrease gene expression from stall-containing mRNAs only at high initiation rates. Remarkably, hundreds of S. cerevisiae mRNAs that contain ribosome stall sequences also exhibit lower translation efficiency. We propose that inefficient translation initiation allows these stall-containing endogenous mRNAs to escape collision-stimulated reduction in gene expression.

Higher rates of translation counterintuitively lead to lower protein levels from eukaryotic mRNAs that encode ribosome stalls; modelling suggests that this occurs when ribosome collisions at stalls trigger abortive termination of the leading ribosome or cause endonucleolytic mRNA cleavage.

Destabilization of Eukaryote mRNAs by 5' Proximal Stop Codons Can Occur Independently of the Nonsense-Mediated mRNA Decay Pathway

In eukaryotes, the binding of poly(A) binding protein (PAB) to the poly(A) tail is central to maintaining mRNA stability. PABP interacts with the translation termination apparatus, and with eIF4G to maintain 3′–5′ mRNA interactions as part of an mRNA closed loop. It is however unclear how ribosome recycling on a closed loop mRNA is influenced by the proximity of the stop codon to the poly(A) tail, and how post-termination ribosome recycling affects mRNA stability. We show that in a yeast disabled for nonsense mediated mRNA decay (NMD), a PGK1 mRNA with an early stop codon at codon 22 of the reading frame is still highly unstable, and that this instability cannot be significantly countered even when 50% stop codon readthrough is triggered. In an NMD-deficient mutant yeast, stable reporter alleles with more 3′ proximal stop codons could not be rendered unstable through Rli1-depletion, inferring defective Rli1 ribosome recycling is insufficient in itself to trigger mRNA instability. Mathematical modelling of a translation system including the effect of ribosome recycling and poly(A) tail shortening supports the hypothesis that impaired ribosome recycling from 5′ proximal stop codons may compromise initiation processes and thus destabilize the mRNA. A model is proposed wherein ribosomes undergo a maturation process during early elongation steps, and acquire competency to re-initiate on the same mRNA as translation elongation progresses beyond the very 5′ proximal regions of the mRNA.

An mRNA decapping mutant deficient in P body assembly limits mRNA stabilization in response to osmotic stress

Yeast is exposed to changing environmental conditions and must adapt its genetic program to provide a homeostatic intracellular environment. An important stress for yeast in the wild is high osmolarity. A key response to this stress is increased mRNA stability primarily by the inhibition of deadenylation. We previously demonstrated that mutations in decapping activators (edc3∆ lsm4∆C), which result in defects in P body assembly, can destabilize mRNA under unstressed conditions. We wished to examine whether mRNA would be destabilized in the edc3∆ lsm4∆C mutant as compared to the wild-type in response to osmotic stress, when P bodies are intense and numerous. Our results show that the edc3∆ lsm4∆C mutant limits the mRNA stability in response to osmotic stress, while the magnitude of stabilization was similar as compared to the wild-type. The reduced mRNA stability in the edc3∆ lsm4∆C mutant was correlated with a shorter PGK1 poly(A) tail. Similarly, the MFA2 mRNA was more rapidly deadenylated as well as significantly stabilized in the ccr4∆ deadenylation mutant in the edc3∆ lsm4∆C background. These results suggest a role for these decapping factors in stabilizing mRNA and may implicate P bodies as sites of reduced mRNA degradation.

The decapping activator Edc3 and the Q/N-rich domain of Lsm4 function together to enhance mRNA stability and alter mRNA decay pathway dependence in Saccharomyces cerevisiae

The rate and regulation of mRNA decay are major elements in the proper control of gene expression. Edc3 and Lsm4 are two decapping activator proteins that have previously been shown to function in the assembly of RNA granules termed P bodies. Here, we show that deletion of edc3, when combined with a removal of the glutamine/asparagine rich region of Lsm4 (edc3Δ lsm4ΔC) reduces mRNA stability and alters pathways of mRNA degradation. Multiple tested mRNAs exhibited reduced stability in the edc3Δ lsm4ΔC mutant. The destabilization was linked to an increased dependence on Ccr4-mediated deadenylation and mRNA decapping. Unlike characterized mutations in decapping factors that either are neutral or are able to stabilize mRNA, the combined edc3Δ lsm4ΔC mutant reduced mRNA stability. We characterized the growth and activity of the major mRNA decay systems and translation in double mutant and wild-type yeast. In the edc3Δ lsm4ΔC mutant, we observed alterations in the levels of specific mRNA decay factors as well as nuclear accumulation of the catalytic subunit of the decapping enzyme Dcp2. Hence, we suggest that the effects on mRNA stability in the edc3Δ lsm4ΔC mutant may originate from mRNA decay protein abundance or changes in mRNPs, or alternatively may imply a role for P bodies in mRNA stabilization.

Summary: A strain mutated in two decapping activators, previously implicated in P body assembly, has reduced mRNA stability and increased dependence on decapping and Ccr4-dependent deadenylation for mRNA degradation.

Stochastic modeling of biochemical systems with multistep reactions using state-dependent time delay

To deal with the growing scale of molecular systems, sophisticated modelling techniques have been designed in recent years to reduce the complexity of mathematical models. Among them, a widely used approach is delayed reaction for simplifying multistep reactions. However, recent research results suggest that a delayed reaction with constant time delay is unable to describe multistep reactions accurately. To address this issue, we propose a novel approach using state-dependent time delay to approximate multistep reactions. We first use stochastic simulations to calculate time delay arising from multistep reactions exactly. Then we design algorithms to calculate time delay based on system dynamics precisely. To demonstrate the power of proposed method, two processes of mRNA degradation are used to investigate the function of time delay in determining system dynamics. In addition, a multistep pathway of metabolic synthesis is used to explore the potential of the proposed method to simplify multistep reactions with nonlinear reaction rates. Simulation results suggest that the state-dependent time delay is a promising and accurate approach to reduce model complexity and decrease the number of unknown parameters in the models.

Modeling and gene knockdown to assess the contribution of nonsense-mediated decay, premature termination, and selenocysteine insertion to the selenoprotein hierarchy

The expression of selenoproteins, a specific group of proteins that incorporates selenocysteine, is hierarchically regulated by the availability of Se, with some, but not all selenoprotein mRNA transcripts decreasing in abundance with decreasing Se. Selenocysteine insertion into the peptide chain occurs during translation following recoding of an internal UGA stop codon. There is increasing evidence that this UGA recoding competes with premature translation termination, which is followed by nonsense-mediated decay (NMD) of the transcript. In this study, we tested the hypothesis that the susceptibility of different selenoprotein mRNAs to premature termination during translation and differential sensitivity of selenoprotein transcripts to NMD are major factors in the selenoprotein hierarchy. Selenoprotein transcript abundance was measured in Caco-2 cells using real-time PCR under different Se conditions and the data obtained fitted to mathematical models of selenoprotein translation. A calibrated model that included a combination of differential sensitivity of selenoprotein transcripts to NMD and different frequency of non-NMD related premature translation termination was able to fit all the measurements. The model predictions were tested using SiRNA to knock down expression of the crucial NMD factor UPF1 (up-frameshift protein 1) and selenoprotein mRNA expression. The calibrated model was able to predict the effect of UPF1 knockdown on gene expression for all tested selenoproteins, except SPS2 (selenophosphate synthetase), which itself is essential for selenoprotein synthesis. These results indicate an important role for NMD in the hierarchical regulation of selenoprotein mRNAs, with the exception of SPS2 whose expression is likely regulated by a different mechanism.

Nonsense-Mediated mRNA Decay: Degradation of Defective Transcripts Is Only Part of the Story

Nonsense-mediated mRNA decay (NMD) is a eukaryotic surveillance mechanism that monitors cytoplasmic mRNA translation and targets mRNAs undergoing premature translation termination for rapid degradation. From yeasts to humans, activation of NMD requires the function of the three conserved Upf factors: Upf1, Upf2, and Upf3. Here, we summarize the progress in our understanding of the molecular mechanisms of NMD in several model systems and discuss recent experiments that address the roles of Upf1, the principal regulator of NMD, in the initial targeting and final degradation of NMD-susceptible mRNAs. We propose a unified model for NMD in which the Upf factors provide several functions during premature termination, including the stimulation of release factor activity and the dissociation and recycling of ribosomal subunits. In this model, the ultimate degradation of the mRNA is the last step in a complex premature termination process.

Kinetics effects and modeling of mRNA turnover

Broader comprehension of gene expression regulatory mechanisms can be gained from a global analysis of how transcription and degradation are coordinated to orchestrate complex cell responses. The role of messenger RNA (mRNA) turnover modulation in gene expression levels has become increasingly recognized. From such perspective, in this review we briefly illustrate how a simple but effective mathematical model of mRNA turnover and some experimental findings, may together shed light on the molecular mechanisms underpinning the major role of mRNA decay rates in shaping the kinetics of gene activation and repression.

Systematic analysis of the role of RNA-binding proteins in the regulation of RNA stability

mRNA half-lives are transcript-specific and vary over a range of more than 100-fold in eukaryotic cells. mRNA stabilities can be regulated by sequence-specific RNA-binding proteins (RBPs), which bind to regulatory sequence elements and modulate the interaction of the mRNA with the cellular RNA degradation machinery. However, it is unclear if this kind of regulation is sufficient to explain the large range of mRNA stabilities. To address this question, we examined the transcriptome of 74 Schizosaccharomyces pombe strains carrying deletions in non-essential genes encoding predicted RBPs (86% of all such genes). We identified 25 strains that displayed changes in the levels of between 4 and 104 mRNAs. The putative targets of these RBPs formed biologically coherent groups, defining regulons involved in cell separation, ribosome biogenesis, meiotic progression, stress responses and mitochondrial function. Moreover, mRNAs in these groups were enriched in specific sequence motifs in their coding sequences and untranslated regions, suggesting that they are coregulated at the posttranscriptional level. We performed genome-wide RNA stability measurements for several RBP mutants, and confirmed that the altered mRNA levels were caused by changes in their stabilities. Although RBPs regulate the decay rates of multiple regulons, only 16% of all S. pombe mRNAs were affected in any of the 74 deletion strains. This suggests that other players or mechanisms are required to generate the observed range of RNA half-lives of a eukaryotic transcriptome.

Virus-mediated compartmentalization of the host translational machinery

Viruses require the host translational apparatus to synthesize viral proteins. Host stress response mechanisms that suppress translation, therefore, represent a significant obstacle that viruses must overcome. Here, we report a strategy whereby the mammalian orthoreoviruses compartmentalize the translational machinery within virus-induced inclusions known as viral factories (VF). VF are the sites of reovirus replication and assembly but were thought not to contain ribosomes. It was assumed viral mRNAs exited the VF to undergo translation by the cellular machinery, and proteins reentered the factory to participate in assembly. Here, we used ribopuromycylation to visualize active translation in infected cells. These studies revealed that active translation occurs within VF and that ribosomal subunits and proteins required for translation initiation, elongation, termination, and recycling localize to the factory. Interestingly, we observed components of the 43S preinitiation complex (PIC) concentrating primarily at factory margins, suggesting a spatial and/or dynamic organization of translation within the VF. Similarly, the viral single-stranded RNA binding protein σNS localized to the factory margins and had a tubulovesicular staining pattern that extended a short distance from the margins of the factories and colocalized with endoplasmic reticulum (ER) markers. Consistent with these colocalization studies, σNS was found to associate with both eukaryotic translation initiation factor 3 subunit A (eIF3A) and the ribosomal subunit pS6R. Together, these findings indicate that σNS functions to recruit 43S PIC machinery to the primary site of viral translation within the viral factory. Pathogen-mediated compartmentalization of the translational apparatus provides a novel mechanism by which viruses might avoid host translational suppression.

Viruses lack biosynthetic capabilities and depend upon the host for protein synthesis. This dependence requires viruses to evolve mechanisms to coerce the host translational machinery into synthesizing viral proteins in the face of ongoing cellular stress responses that suppress global protein synthesis. Reoviruses replicate and assemble within cytoplasmic inclusions called viral factories. However, synthesis of viral proteins was thought to occur in the cytosol. To identify the site(s) of viral translation, we undertook a microscopy-based approach using ribopuromycylation to detect active translation. Here, we report that active translation occurs within viral factories and that translational factors are compartmentalized within factories. Furthermore, we find that the reovirus nonstructural protein σNS associates with 43S preinitiation complexes at the factory margins, suggesting a role for σNS in translation. Together, virus-induced compartmentalization of the host translational machinery represents a strategy for viruses to spatiotemporally couple viral protein synthesis with viral replication and assembly.

A halt in poly(A) shortening during S-adenosyl-L-methionine-induced translation arrest in CGS1 mRNA of Arabidopsis thaliana

Cystathionine γ-synthase (CGS) catalyzes the first committed step of methionine (Met) biosynthesis in plants. Expression of the Arabidopsis thaliana CGS1 gene is negatively feedback-regulated in response to the direct Met metabolite S-adenosyl-L-methionine (AdoMet). This regulation occurs at the step of mRNA stability during translation and is coupled with AdoMet-induced CGS1-specific translation arrest. In general, mRNA decay is initiated by a shortening of the poly(A) tail. However, this process has not been studied in detail in cases where regulatory events, such as programmed translation arrest, are involved. Here, we report that the poly(A) tail of the full-length CGS1 mRNA showed an apparent increase from 50-80 nucleotides (nt) to 140-150 nt after the induction of CGS1 mRNA degradation. This finding was unexpected because mRNAs that are destined for degradation harbored longer poly(A) tail than mRNAs that were not targeted for degradation. The results suggest that poly(A) shortening is inhibited or delayed during AdoMet-induced translation arrest of CGS1 mRNA. We propose an explanation for this phenomenon that remains consistent with the recent model of actively translating mRNA. We also found that CGS1 mRNA degradation intermediates, which are 5'-truncated forms of CGS1 mRNA, had a short poly(A) tail of 10-30 nt. This suggests that poly(A) shortening occurs rapidly on the degradation intermediates. The present study highlights CGS1 mRNA degradation as a useful system to understand the dynamic features of poly(A) shortening.

RNA biology and the adaptation of Cryptococcus neoformans to host temperature and stress

Cryptococcus neoformans is an environmental fungus that can cause severe disease in humans. C. neoformans encounters a multitude of stresses within the human host to which it must adapt in order to survive and proliferate. Upon stressful changes in the external milieu, C. neoformans must reprogram its gene expression to properly respond to and combat stress in order to maintain homeostasis. Several studies have investigated the changes that occur in response to these stresses to begin to unravel the mechanisms of adaptation in this organism. Here, we review studies that have explored stress-induced changes in gene expression with a focus on host temperature adaptation. We compare global messenger RNA (mRNA) expression data compiled from several studies and identify patterns that suggest that orchestrated, transient responses occur. We also utilize the available expression data to explore the possibility of a common stress response that may contribute to cellular protection against a variety of stresses in C. neoformans. In addition, we review studies that have revealed the significance of post-transcriptional mechanisms of mRNA regulation in response to stress, and discuss how these processes may contribute to adaptation and virulence.

Lsm2 and Lsm3 bridge the interaction of the Lsm1-7 complex with Pat1 for decapping activation

The evolutionarily conserved Lsm1-7-Pat1 complex is the most critical activator of mRNA decapping in eukaryotic cells and plays many roles in normal decay, AU-rich element-mediated decay, and miRNA silencing, yet how Pat1 interacts with the Lsm1-7 complex is unknown. Here, we show that Lsm2 and Lsm3 bridge the interaction between the C-terminus of Pat1 (Pat1C) and the Lsm1-7 complex. The Lsm2-3-Pat1C complex and the Lsm1-7-Pat1C complex stimulate decapping in vitro to a similar extent and exhibit similar RNA-binding preference. The crystal structure of the Lsm2-3-Pat1C complex shows that Pat1C binds to Lsm2-3 to form an asymmetric complex with three Pat1C molecules surrounding a heptameric ring formed by Lsm2-3. Structure-based mutagenesis revealed the importance of Lsm2-3-Pat1C interactions in decapping activation in vivo. Based on the structure of Lsm2-3-Pat1C, a model of Lsm1-7-Pat1 complex is constructed and how RNA binds to this complex is discussed.

Stochastic modelling of biochemical systems of multi-step reactions using a simplified two-variable model

Background

A fundamental issue in systems biology is how to design simplified mathematical models for describing the dynamics of complex biochemical reaction systems. Among them, a key question is how to use simplified reactions to describe the chemical events of multi-step reactions that are ubiquitous in biochemistry and biophysics. To address this issue, a widely used approach in literature is to use one-step reaction to represent the multi-step chemical events. In recent years, a number of modelling methods have been designed to improve the accuracy of the one-step reaction method, including the use of reactions with time delay. However, our recent research results suggested that there are still deviations between the dynamics of delayed reactions and that of the multi-step reactions. Therefore, more sophisticated modelling methods are needed to accurately describe the complex biological systems in an efficient way.

Results

This work designs a two-variable model to simplify chemical events of multi-step reactions. In addition to the total molecule number of a species, we first introduce a new concept regarding the location of molecules in the multi-step reactions, which is the second variable to represent the system dynamics. Then we propose a simulation algorithm to compute the probability for the firing of the last step reaction in the multi-step events. This probability function is evaluated using a deterministic model of ordinary differential equations and a stochastic model in the framework of the stochastic simulation algorithm. The efficiency of the proposed two-variable model is demonstrated by the realization of mRNA degradation process based on the experimentally measured data.

Conclusions

Numerical results suggest that the proposed new two-variable model produces predictions that match the multi-step chemical reactions very well. The successful realization of the mRNA degradation dynamics indicates that the proposed method is a promising approach to reduce the complexity of biological systems.

The RRM1 domain of the poly(A)-binding protein from Saccharomyces cerevisiae is critical to control of mRNA deadenylation

The poly(A)-binding protein PAB1 from the yeast Saccharomyces cerevisiae plays an important role in controlling mRNA deadenylation rates. Deletion of either its RRM1 or proline-rich domain (P domain) severely restricts deadenylation and slows mRNA degradation. Because these large deletions could be having unknown effects on the structure of PAB1, different strategies were used to determine the importance of the RRM1 and P domains to deadenylation. Since the P domain is quite variable in size and sequence among eukaryotes, P domains from two human PABPCs and from Xenopus were substituted for that of PAB1. The resultant PAB1 hybrid proteins, however, displayed limited or no difference in mRNA deadenylation as compared with PAB1. In contrast to the P domain, the RRM1 domain is highly conserved across species, and a systematic mutagenesis of the RRM1 domain was undertaken to identify its functional regions. Several mutations along the RNA-binding surface of RRM1 inhibited deadenylation, whereas one set of mutations on its exterior non-RNA binding surface shifted deadenylation from a slow distributive process to a rapid processive deadenylation. These results suggest that the RRM1 domain is the more critical region of PAB1 for controlling deadenylation and consists of at least two distinguishable functional regions.

Eukaryotic mRNA decay: methodologies, pathways, and links to other stages of gene expression

mRNA concentration depends on the balance between transcription and degradation rates. On both sides of the equilibrium, synthesis and degradation show, however, interesting differences that have conditioned the evolution of gene regulatory mechanisms. Here, we discuss recent genome-wide methods for determining mRNA half-lives in eukaryotes. We also review pre- and posttranscriptional regulons that coordinate the fate of functionally related mRNAs by using protein- or RNA-based trans factors. Some of these factors can regulate both transcription and decay rates, thereby maintaining proper mRNA homeostasis during eukaryotic cell life.

Circadian control of mRNA polyadenylation dynamics regulates rhythmic protein expression

Poly(A) tails are 3' modifications of eukaryotic mRNAs that are important in the control of translation and mRNA stability. We identified hundreds of mouse liver mRNAs that exhibit robust circadian rhythms in the length of their poly(A) tails. Approximately 80% of these are primarily the result of nuclear adenylation coupled with rhythmic transcription. However, unique decay kinetics distinguish these mRNAs from other mRNAs that are transcribed rhythmically but do not exhibit poly(A) tail rhythms. The remaining 20% are uncoupled from transcription and exhibit poly(A) tail rhythms even though the steady-state mRNA levels are not rhythmic. These are under the control of rhythmic cytoplasmic polyadenylation, regulated at least in some cases by cytoplasmic polyadenylation element-binding proteins (CPEBs). Importantly, we found that the rhythmicity in poly(A) tail length is closely correlated with rhythmic protein expression, with a several-hour delay between the time of longest tail and the time of highest protein level. Our study demonstrates that the circadian clock regulates the dynamic polyadenylation status of mRNAs, which can result in rhythmic protein expression independent of the steady-state levels of the message.

Deadenylation and P-bodies

Deadenylation is the major step in triggering mRNA decay and results in mRNA translation inhibition in eukaryotic cells. Therefore, it is plausible that deadenylation also induces the mRNP remodeling required for formation of GW bodies or RNA processing bodies (P-bodies), which harbor translationally silenced mRNPs. In this chapter, we discuss several examples to illustrate the roles of deadenylation in regulating gene expression. We highlight several lines of evidence indicating that even though non-translatable mRNPs may be prepared and/or assembled into P-bodies in different ways, deadenylation is always a necessary, and perhaps the earliest, step in mRNA decay pathways that enable mRNP remodeling required for P-body formation. Thus, deadenylation and the participating deadenylases are not simply required for preparing mRNA substrates; they play an indispensable role both structurally and functionally in P-body formation and regulation.

Regulation of release factor expression using a translational negative feedback loop: a systems analysis

The essential eukaryote release factor eRF1, encoded by the yeast SUP45 gene, recognizes stop codons during ribosomal translation. SUP45 nonsense alleles are, however, viable due to the establishment of feedback-regulated readthrough of the premature termination codon; reductions in full-length eRF1 promote tRNA-mediated stop codon readthrough, which, in turn, drives partial production of full-length eRF1. A deterministic mathematical model of this eRF1 feedback loop was developed using a staged increase in model complexity. Model predictions matched the experimental observation that strains carrying the mutant SUQ5 tRNA (a weak UAA suppressor) in combination with any of the tested sup45(UAA) nonsense alleles exhibit threefold more stop codon readthrough than that of an SUQ5 yeast strain. The model also successfully predicted that eRF1 feedback control in an SUQ5 sup45(UAA) mutant would resist, but not completely prevent, imposed changes in eRF1 expression. In these experiments, the introduction of a plasmid-borne SUQ5 copy into a sup45(UAA) SUQ5 mutant directed additional readthrough and full-length eRF1 expression, despite feedback. Secondly, induction of additional sup45(UAA) mRNA expression in a sup45(UAA) SUQ5 strain also directed increased full-length eRF1 expression. The autogenous sup45 control mechanism therefore acts not to precisely control eRF1 expression, but rather as a damping mechanism that only partially resists changes in release factor expression level. The validated model predicts that the degree of feedback damping (i.e., control precision) is proportional to eRF1 affinity for the premature stop codon. The validated model represents an important tool to analyze this and other translational negative feedback loops.

Architecture of the nuclease module of the yeast Ccr4-not complex: the Not1-Caf1-Ccr4 interaction

Shortening eukaryotic poly(A) tails represses mRNA translation and induces mRNA turnover. The major cytoplasmic deadenylase, the Ccr4-Not complex, is a conserved multisubunit assembly. Ccr4-Not is organized around Not1, a large scaffold protein that recruits two 3'-5' exoribonucleases, Caf1 and Ccr4. We report structural studies showing that the N-terminal arm of yeast Not1 has a HEAT-repeat structure with domains related to the MIF4G fold. A MIF4G domain positioned centrally within the Not1 protein recognizes Caf1, which in turn binds the LRR domain of Ccr4 and tethers the Ccr4 nuclease domain. The interactions that form the nuclease core of the Ccr4-Not complex are evolutionarily conserved. Their specific disruption affects cell growth and mRNA deadenylation and decay in vivo in yeast. Thus, the N-terminal arm of Not1 forms an extended platform reminiscent of scaffolding proteins like eIF4G and CBP80, and places the two nucleases in a pivotal position within the Ccr4-Not complex.

RNA degradation in Saccharomyces cerevisae

All RNA species in yeast cells are subject to turnover. Work over the past 20 years has defined degradation mechanisms for messenger RNAs, transfer RNAs, ribosomal RNAs, and noncoding RNAs. In addition, numerous quality control mechanisms that target aberrant RNAs have been identified. Generally, each decay mechanism contains factors that funnel RNA substrates to abundant exo- and/or endonucleases. Key issues for future work include determining the mechanisms that control the specificity of RNA degradation and how RNA degradation processes interact with translation, RNA transport, and other cellular processes.

Degradation of mRNAs that lack a stop codon: a decade of nonstop progress

Nonstop decay is the mechanism of identifying and disposing aberrant transcripts that lack in-frame stop codons. It is hypothesized that these transcripts are identified during translation when the ribosome arrives at the 3' end of the mRNA and stalls. Presumably, the ribosome stalling recruits additional cofactors, Ski7 and the exosome complex. The exosome degrades the transcript using either one of its ribonucleolytic activities, and the ribosome and the peptide are both released. Additional precautionary measures by the nonstop decay pathway may include translational repression of the nonstop transcript after translation, and proteolysis of the released peptide by the proteasome. This surveillance mechanism protects the cells from potentially harmful truncated proteins, but it may also be involved in mediating critical cellular functions of transcripts that are prone to stop codon read-through. Important advances have been made in the past decade as we learn that nonstop decay may have implications in human disease.

Scd6 targets eIF4G to repress translation: RGG motif proteins as a class of eIF4G-binding proteins

The formation of mRNPs controls the interaction of the translation and degradation machinery with individual mRNAs. The yeast Scd6 protein and its orthologs regulate translation and mRNA degradation in yeast, C. elegans, D. melanogaster, and humans by an unknown mechanism. We demonstrate that Scd6 represses translation by binding the eIF4G subunit of eIF4F in a manner dependent on its RGG domain, thereby forming an mRNP repressed for translation initiation. Strikingly, several other RGG domain-containing proteins in yeast copurify with eIF4E/G and we demonstrate that two such proteins, Npl3 and Sbp1, also directly bind eIF4G and repress translation in a manner dependent on their RGG motifs. These observations identify the mechanism of Scd6 function through its RGG motif and indicate that eIF4G plays an important role as a scaffolding protein for the recruitment of translation repressors.

Promoter elements regulate cytoplasmic mRNA decay

Promoters are DNA elements that enable transcription and its regulation by trans-acting factors. Here, we demonstrate that yeast promoters can also regulate mRNA decay after the mRNA leaves the nucleus. A conventional yeast promoter consists of a core element and an upstream activating sequence (UAS). We find that changing UASs of a reporter gene without altering the transcript sequence affects the transcript's decay kinetics. A short cis element, comprising two Rap1p-binding sites, and Rap1p itself, are necessary and sufficient to induce enhanced decay of the reporter mRNA. Furthermore, Rap1p stimulates both the synthesis and the decay of a specific population of endogenous mRNAs. We propose that Rap1p association with target promoter in the nucleus affects the composition of the exported mRNP, which in turn regulates mRNA decay in the cytoplasm. Thus, promoters can play key roles in determining mRNA levels and have the capacity to coordinate rates of mRNA synthesis and decay.

SPT5 affects the rate of mRNA degradation and physically interacts with CCR4 but does not control mRNA deadenylation

The CCR4-NOT complex has been shown to have multiple roles in mRNA metabolism, including that of transcriptional elongation, mRNA transport, and nuclear exosome function, but the primary function of CCR4 and CAF1 is in the deadenylation and degradation of cytoplasmic mRNA. As previous genetic analysis supported an interaction between SPT5, known to be involved in transcriptional elongation, and that of CCR4, the physical association of SPT5 with CCR4 was examined. A two-hybrid screen utilizing the deadenylase domain of CCR4 as a bait identified SPT5 as a potential interacting protein. SPT5 at its physiological concentration was shown to immunoprecipitate CCR4 and CAF1, and in vitro purified SPT5 specifically could bind to CAF1 and the deadenylase domain of CCR4. We additionally demonstrated that mutations in SPT5 or an spt4 deletion slowed the rate of mRNA degradation, a phenotype associated with defects in the CCR4 mRNA deadenylase complex. Yet, unlike ccr4 and caf1 deletions, spt5 and spt4 defects displayed little effect on the rate of deadenylation. They also did not affect decapping or 5' - 3' degradation of mRNA. These results suggest that the interactions between SPT5/SPT4 and the CCR4-NOT complex are probably the consequences of effects involving nuclear events and do not involve the primary role of CCR4 in mRNA deadenylation and turnover.

The circadian deadenylase Nocturnin is necessary for stabilization of the iNOS mRNA in mice

Nocturnin is a member of the CCR4 deadenylase family, and its expression is under circadian control with peak levels at night. Because it can remove poly(A) tails from mRNAs, it is presumed to play a role in post-transcriptional control of circadian gene expression, but its target mRNAs are not known. Here we demonstrate that Nocturnin expression is acutely induced by the endotoxin lipopolysaccharide (LPS). Mouse embryo fibroblasts (MEFs) lacking Nocturnin exhibit normal patterns of acute induction of TNFα and iNOS mRNAs during the first three hours following LPS treatment, but by 24 hours, while TNFα mRNA levels are indistinguishable from WT cells, iNOS message is significantly reduced 20-fold. Accordingly, analysis of the stability of the mRNAs showed that loss of Nocturnin causes a significant decrease in the half-life of the iNOS mRNA (t(1/2) = 3.3 hours in Nocturnin knockout MEFs vs. 12.4 hours in wild type MEFs), while having no effect on the TNFα message. Furthermore, mice lacking Nocturnin lose the normal nighttime peak of hepatic iNOS mRNA, and have improved survival following LPS injection. These data suggest that Nocturnin has a novel stabilizing activity that plays an important role in the circadian response to inflammatory signals.

A systems view of the protein expression process

Many biological processes are regulated by changing the concentration and activity of proteins. The presence of a protein at a given subcellular location at a given time with a certain conformation is the result of an apparently sequential process. The rate of protein formation is influenced by chromatin state, and the rates of transcription, translation, and degradation. There is an exquisite control system where each stage of the process is controlled both by seemingly unregulated proteins as well as through feedbacks mediated by RNA and protein products. Here we review the biological facts and mathematical models for each stage of the protein production process. We conclude that advances in experimental techniques leading to a detailed description of the process have not been matched by mathematical models that represent the details of the process and facilitate analysis. Such an exercise is the first step towards development of a framework for a systems biology analysis of the protein production process.

ELECTRONIC SUPPLEMENTARY MATERIAL

The online version of this article (doi:10.1007/s11693-011-9088-1) contains supplementary material, which is available to authorized users.

A quantitative assay for measuring mRNA decapping by splinted ligation reverse transcription polymerase chain reaction: qSL-RT-PCR

The degradation of messenger RNA is a critical node of gene regulation. A major pathway of mRNA decay is initiated by shortening of the poly(A) tail, followed by removal of the 5' cap structure (decapping) and subsequent degradation. Decapping is an important determinate in the destruction of many transcripts. Detailed kinetic analysis of in vivo decapping rates is necessary to understand how this step is regulated. Importantly, the product of decapping is recalcitrant for investigation, in part due to its transient nature. As such, little in vivo kinetic information is available. Here we report the development of an assay that measures decapping of mRNAs by combining splinted ligation and quantitative RT-PCR (qSL-RT-PCR). We apply this method to determine the decapping rate constant for a natural mRNA in vivo for the first time. The qSL-RT-PCR assay may be adapted for use on any mRNA, providing a new tool to study regulation of mRNA decay.

Different nuclease requirements for exosome-mediated degradation of normal and nonstop mRNAs

Two general pathways of mRNA decay have been characterized in yeast. In one pathway, the mRNA is degraded by the cytoplasmic form of the exosome. The exosome has both 3' to 5' exoribonuclease and endoribonuclease activity, and the available evidence suggests that the exonuclease activity is required for the degradation of mRNAs. We confirm here that this is true for normal mRNAs, but that aberrant mRNAs that lack a stop codon can be efficiently degraded in the absence of the exonuclease activity of the exosome. Specifically, we show that the endo- and exonuclease activities of the exosome are both capable of rapidly degrading nonstop mRNAs and ribozyme-cleaved mRNAs. Additionally, the endonuclease activity of the exosome is not required for endonucleolytic cleavage in no-go decay. In vitro, the endonuclease domain of the exosome is active only under nonphysiological conditions, but our findings show that the in vivo activity is sufficient for the rapid degradation of nonstop mRNAs. Thus, whereas normal mRNAs are degraded by two exonucleases (Xrn1p and Rrp44p), several endonucleases contribute to the decay of many aberrant mRNAs, including transcripts subject to nonstop and no-go decay. Our findings suggest that the nuclease requirements for general and nonstop mRNA decay are different, and describe a molecular function of the core exosome that is not disrupted by inactivating its exonuclease activity.