Live-Cell Single RNA Imaging Reveals Bursts of Translational Frameshifting.

Ribosomal frameshifting during the translation of RNA is implicated in human disease and viral infection. While previous work has uncovered many details about single RNA frameshifting kinetics in vitro, little is known about how single RNA frameshift in living systems. To confront this problem, we have developed technology to quantify live-cell single RNA translation dynamics in frameshifted open reading frames. Applying this technology to RNA encoding the HIV-1 frameshift sequence reveals a small subset (∼8%) of the translating pool robustly frameshift. Frameshifting RNA are translated at similar rates as non-frameshifting RNA (∼3 aa/s) and can continuously frameshift for more than four rounds of translation. Fits to a bursty model of frameshifting constrain frameshifting kinetic rates and demonstrate how ribosomal traffic jams contribute to the persistence of the frameshifting state. These data provide insight into retroviral frameshifting and could lead to alternative strategies to perturb the process in living cells.

Ten simple rules for managing laboratory information.

Information is the cornerstone of research, from experimental (meta)data and computational processes to complex inventories of reagents and equipment. These 10 simple rules discuss best practices for leveraging laboratory information management systems to transform this large information load into useful scientific findings.

Stochastic dynamics of Type-I interferon responses.

Interferon (IFN) activates the transcription of several hundred of IFN stimulated genes (ISGs) that constitute a highly effective antiviral defense program. Cell-to-cell variability in the induction of ISGs is well documented, but its source and effects are not completely understood. The molecular mechanisms behind this heterogeneity have been related to randomness in molecular events taking place during the JAK-STAT signaling pathway. Here, we study the sources of variability in the induction of the IFN-alpha response by using MxA and IFIT1 activation as read-out. To this end, we integrate time-resolved flow cytometry data and stochastic modeling of the JAK-STAT signaling pathway. The complexity of the IFN response was matched by fitting probability distributions to time-course flow cytometry snapshots. Both, experimental data and simulations confirmed that the MxA and IFIT1 induction circuits generate graded responses rather than all-or-none responses. Subsequently, we quantify the size of the intrinsic variability at different steps in the pathway. We found that stochastic effects are transiently strong during the ligand-receptor activation steps and the formation of the ISGF3 complex, but negligible for the final induction of the studied ISGs. We conclude that the JAK-STAT signaling pathway is a robust biological circuit that efficiently transmits information under stochastic environments.

Modeling the intracellular dynamics of the dengue viral infection and the innate immune response.

The interplay between the dengue virus and the innate immune response is not fully understood. Here, we use deterministic and stochastic approaches to investigate the dynamics of the interaction between the interferon-mediated innate immune response and the dengue virus. We aim to develop a quantitative representation of these complex interactions and predict their system-level dynamics. Our simulation results predict bimodal and bistable dynamics that represent viral clearance and virus-producing states. Under normal conditions, we determined that the viral infection outcome is modulated by the innate immune response and the positive-strand viral RNA concentration. Additionally, we tested system perturbations by external stimulation, such as the direct induction of the innate immune response by interferon, and a therapeutic intervention consisting of the direct application of mRNA encoding for several interferon-stimulated genes. Our simulation results suggest optimal regimes for the studied intervention approaches.

Computational design and interpretation of single-RNA translation experiments.

Advances in fluorescence microscopy have introduced new assays to quantify live-cell translation dynamics at single-RNA resolution. We introduce a detailed, yet efficient sequence-based stochastic model that generates realistic synthetic data for several such assays, including Fluorescence Correlation Spectroscopy (FCS), ribosome Run-Off Assays (ROA) after Harringtonine application, and Fluorescence Recovery After Photobleaching (FRAP). We simulate these experiments under multiple imaging conditions and for thousands of human genes, and we evaluate through simulations which experiments are most likely to provide accurate estimates of elongation kinetics. Finding that FCS analyses are optimal for both short and long length genes, we integrate our model with experimental FCS data to capture the nascent protein statistics and temporal dynamics for three human genes: KDM5B, ß-actin, and H2B. Finally, we introduce a new open-source software package, RNA Sequence to NAscent Protein Simulator (rSNAPsim), to easily simulate the single-molecule translation dynamics of any gene sequence for any of these assays and for different assumptions regarding synonymous codon usage, tRNA level modifications, or ribosome pauses. rSNAPsim is implemented in Python and is available at: https://github.com/MunskyGroup/rSNAPsim.git.

Modeling the effect of tat inhibitors on HIV latency.

Even in the presence of a successful combination therapy stalling the progress of AIDS, developing a cure for this disease is still an open question. One of the major steps towards a cure would be to be able to eradicate latent HIV reservoirs present in patients. During the last decade, multiple findings point to the dominant role of the viral protein Tat in the establishment of latency. Here we present a mathematical study to understand the potential role of Tat inhibitors as virus-suppressing agents. For this aim, we implemented a computational model that reproduces intracellular dynamics. Simulating an HIV-infected cell and its intracellular feedback we observed that removing Tat protein from the system via inhibitors resulted in a temporary and reversible viral suppression. In contrast, we observed that compounds that interact with Tat protein and disrupt the integrated viral genome produced a more permanent viral suppression.

Studying HIV latency by modeling the interaction between HIV proteins and the innate immune response.

HIV infection leads to two cell fates, the viral productive state or viral latency (a reversible non-productive state). HIV latency is relevant because infected active CD4+ T-lymphocytes can reach a resting memory state in which the provirus remains silent for long periods of time. Despite experimental and theoretical efforts, the causal molecular mechanisms responsible for HIV latency are only partially understood. Studies have determined that HIV latency is influenced by the innate immune response carried out by cell restriction factors that inhibit the postintegration steps in the virus replication cycle. In this study, we present a mathematical study that combines deterministic and stochastic approaches to analyze the interactions between HIV proteins and the innate immune response. Using wide ranges of parameter values, we observed the following: (1) a phenomenological description of the viral productive and latent cell phenotypes is obtained by bistable and bimodal dynamics, (2) biochemical noise reduces the probability that an infected cell adopts the latent state, (3) the effects of the innate immune response enhance the HIV latency state, (4) the conditions of the cell before infection affect the latent phenotype, i.e., the existing expression of cell restriction factors propitiates HIV latency, and existing expression of HIV proteins reduces HIV latency.

Analysis and design of single-cell experiments to harvest fluctuation information while rejecting measurement noise.

Introduction: Despite continued technological improvements, measurement errors always reduce or distort the information that any real experiment can provide to quantify cellular dynamics. This problem is particularly serious for cell signaling studies to quantify heterogeneity in single-cell gene regulation, where important RNA and protein copy numbers are themselves subject to the inherently random fluctuations of biochemical reactions. Until now, it has not been clear how measurement noise should be managed in addition to other experiment design variables (e.g., sampling size, measurement times, or perturbation levels) to ensure that collected data will provide useful insights on signaling or gene expression mechanisms of interest. Methods: We propose a computational framework that takes explicit consideration of measurement errors to analyze single-cell observations, and we derive Fisher Information Matrix (FIM)-based criteria to quantify the information value of distorted experiments. Results and Discussion: We apply this framework to analyze multiple models in the context of simulated and experimental single-cell data for a reporter gene controlled by an HIV promoter. We show that the proposed approach quantitatively predicts how different types of measurement distortions affect the accuracy and precision of model identification, and we demonstrate that the effects of these distortions can be mitigated through explicit consideration during model inference. We conclude that this reformulation of the FIM could be used effectively to design single-cell experiments to optimally harvest fluctuation information while mitigating the effects of image distortion.

Using Mechanistic Models and Machine Learning to Design Single-Color Multiplexed Nascent Chain Tracking Experiments.

mRNA translation is the ubiquitous cellular process of reading messenger-RNA strands into functional proteins. Over the past decade, large strides in microscopy techniques have allowed observation of mRNA translation at a single-molecule resolution for self-consistent time-series measurements in live cells. Dubbed Nascent chain tracking (NCT), these methods have explored many temporal dynamics in mRNA translation uncaptured by other experimental methods such as ribosomal profiling, smFISH, pSILAC, BONCAT, or FUNCAT-PLA. However, NCT is currently restricted to the observation of one or two mRNA species at a time due to limits in the number of resolvable fluorescent tags. In this work, we propose a hybrid computational pipeline, where detailed mechanistic simulations produce realistic NCT videos, and machine learning is used to assess potential experimental designs for their ability to resolve multiple mRNA species using a single fluorescent color for all species. Through simulation, we show that with careful application, this hybrid design strategy could in principle be used to extend the number of mRNA species that could be watched simultaneously within the same cell. We present a simulated example NCT experiment with seven different mRNA species within the same simulated cell and use our ML labeling to identify these spots with 90% accuracy using only two distinct fluorescent tags. The proposed extension to the NCT color palette should allow experimentalists to access a plethora of new experimental design possibilities, especially for cell signalling applications requiring simultaneous study of multiple mRNAs.

Using mechanistic models and machine learning to design single-color multiplexed nascent chain tracking experiments.

mRNA translation is the ubiquitous cellular process of reading messenger-RNA strands into functional proteins. Over the past decade, large strides in microscopy techniques have allowed observation of mRNA translation at a single-molecule resolution for self-consistent time-series measurements in live cells. Dubbed Nascent chain tracking (NCT), these methods have explored many temporal dynamics in mRNA translation uncaptured by other experimental methods such as ribosomal profiling, smFISH, pSILAC, BONCAT, or FUNCAT-PLA. However, NCT is currently restricted to the observation of one or two mRNA species at a time due to limits in the number of resolvable fluorescent tags. In this work, we propose a hybrid computational pipeline, where detailed mechanistic simulations produce realistic NCT videos, and machine learning is used to assess potential experimental designs for their ability to resolve multiple mRNA species using a single fluorescent color for all species. Our simulation results show that with careful application this hybrid design strategy could in principle be used to extend the number of mRNA species that could be watched simultaneously within the same cell. We present a simulated example NCT experiment with seven different mRNA species within the same simulated cell and use our ML labeling to identify these spots with 90% accuracy using only two distinct fluorescent tags. We conclude that the proposed extension to the NCT color palette should allow experimentalists to access a plethora of new experimental design possibilities, especially for cell Signaling applications requiring simultaneous study of multiple mRNAs.

What is the core oscillator in the speract-activated pathway of the Strongylocentrotus purpuratus sperm flagellum?

Sperm chemotaxis has an important role in fertilization. Most of our knowledge regarding this phenomenon comes from studies in organisms whose fertilization occurs externally, like sea urchins. Sea urchin spermatozoa respond to sperm-activating peptides, which diffuse from the egg jelly coat and interact with their receptor in the flagellum, triggering several physiological responses: changes in membrane potential, intracellular pH, cyclic nucleotide levels, and intracellular Ca2+ concentration ([Ca2+]). In particular, flagellar [Ca2+] has been shown to oscillate. These [Ca2+] oscillations are correlated with changes in the flagellar shape and so with the regulation of the sperm swimming paths. In this study, we demonstrate, from a mathematical modeling perspective, that the reported speract-activated signaling pathway in Strongylocentrotus purpuratus (speract being a sperm-activating peptide specific to this species) has the necessary elements to replicate the reported [Ca2+] oscillations. We further investigate which elements of this signaling pathway constitute the core oscillator.

A new efficient approach to fit stochastic models on the basis of high-throughput experimental data using a model of IRF7 gene expression as case study.

BACKGROUND: Mathematical models are used to gain an integrative understanding of biochemical processes and networks. Commonly the models are based on deterministic ordinary differential equations. When molecular counts are low, stochastic formalisms like Monte Carlo simulations are more appropriate and well established. However, compared to the wealth of computational methods used to fit and analyze deterministic models, there is only little available to quantify the exactness of the fit of stochastic models compared to experimental data or to analyze different aspects of the modeling results. RESULTS: Here, we developed a method to fit stochastic simulations to experimental high-throughput data, meaning data that exhibits distributions. The method uses a comparison of the probability density functions that are computed based on Monte Carlo simulations and the experimental data. Multiple parameter values are iteratively evaluated using optimization routines. The method improves its performance by selecting parameters values after comparing the similitude between the deterministic stability of the system and the modes in the experimental data distribution. As a case study we fitted a model of the IRF7 gene expression circuit to time-course experimental data obtained by flow cytometry. IRF7 shows bimodal dynamics upon IFN stimulation. This dynamics occurs due to the switching between active and basal states of the IRF7 promoter. However, the exact molecular mechanisms responsible for the bimodality of IRF7 is not fully understood. CONCLUSIONS: Our results allow us to conclude that the activation of the IRF7 promoter by the combination of IRF7 and ISGF3 is sufficient to explain the observed bimodal dynamics.