Nonlatching positive feedback enables robust bimodality by decoupling expression noise from the mean

Mathematical modeling and mechanisms of HIV latency for personalized anti latency therapies

Combination antiretroviral therapy controls human immunodeficiency virus-1 (HIV) but cannot eradicate latent proviruses in immune cells, which reactivate upon treatment interruption. Anti-latency therapies like "shock-and-kill" are being developed but are yet to succeed due to the complexity of latency mechanisms. This review discusses recent advances in understanding HIV latency via mathematical modeling, covering key regulatory factors and models to predict latency reversal, highlighting gaps to guide future therapeutic approaches.

Origin and functional role of antisense transcription in endogenous and exogenous retroviruses

Most proteins expressed by endogenous and exogenous retroviruses are encoded in the sense (positive) strand of the genome and are under the control of regulatory elements within the 5' long terminal repeat (LTR). A number of retroviral genomes also encode genes in the antisense (negative) strand and their expression is under the control of negative sense promoters within the 3' LTR. In the case of the Human T-cell Lymphotropic Virus 1 (HTLV-1), the antisense protein HBZ has been shown to play a critical role in the virus lifecycle and in the pathogenic process, while the function of the Human Immunodeficiency Virus 1 (HIV-1) antisense protein ASP remains unknown. However, the expression of 3' LTR-driven antisense transcripts is not always demonstrably associated with the presence of an antisense open reading frame encoding a viral protein. Moreover, even in the case of retroviruses that do express an antisense protein, such as HTLV-1 and the pandemic strains of HIV-1, the 3' LTR-driven antisense transcript shows both protein-coding and noncoding activities. Indeed, the ability to express antisense transcripts appears to be phylogenetically more widespread among endogenous and exogenous retroviruses than the presence of a functional antisense open reading frame within these transcripts. This suggests that retroviral antisense transcripts may have originated as noncoding molecules with regulatory activity that in some cases later acquired protein-coding function. Here, we will review examples of endogenous and exogenous retroviral antisense transcripts, and the ways through which they benefit viral persistence in the host.

Enhanced Transcriptional Strength of HIV-1 Subtype C Minimizes Gene Expression Noise and Confers Stability to the Viral Latent State

Stochastic fluctuations in gene expression emanating from the HIV-1 long terminal repeat (LTR), amplified by the Tat positive feedback circuit, determine the choice between viral infection fates: active transcription (ON) or transcriptional silence (OFF). The emergence of several transcription factor binding site (TFBS) variant strains in HIV-1 subtype C (HIV-1C), especially those containing the duplication of the NF-κB motif, mandates the evaluation of the effect of enhanced transcriptional strength on gene expression noise and its influence on viral fate selection switch. Using a panel of subgenomic LTR-variant strains containing different copy numbers of the NF-κB motif (ranging from 0 to 4), we used flow cytometry, mRNA quantification, and pharmacological perturbations to demonstrate an inverse correlation between promoter strength and gene expression noise in Jurkat T cells and primary CD4⁺ T cells. The inverse correlation is consistent in clonal cell populations at constant intracellular concentrations of Tat and when NF-κB levels were regulated pharmacologically. Further, we show that strong LTRs containing at least two copies of the NF-κB motif in the enhancer establish a more stable latent state and demonstrate more rapid latency reversal than weak LTRs containing fewer motifs. We also demonstrate a cooperative binding of NF-κB to the motif cluster in HIV-1C LTRs containing two, three, or four NF-κB motifs (Hill coefficient [H] = 2.61, 3.56, and 3.75, respectively). The present work alludes to a possible evolution of the HIV-1C LTR toward gaining transcriptional strength associated with attenuated gene expression noise with implications for viral latency. IMPORTANCE Over the past two consecutive decades, HIV-1 subtype C (HIV-1C) has been undergoing directional evolution toward augmenting the transcriptional strength of the long terminal repeat (LTR) by adding more copies of the existing transcription factor binding site (TFBS) by sequence duplication. Additionally, the duplicated elements are genetically diverse, suggesting broader-range signal receptivity by variant LTRs. The HIV-1 promoter is inherently noisy, and the stochastic fluctuations in gene expression of variant LTRs may influence the active transcription (ON)/transcriptional silence (OFF) latency decisions. The evolving NF-κB motif variations of HIV-1C offer a powerful opportunity to examine how the transcriptional strength of the LTR might influence gene expression noise. Our work here shows that the augmented transcriptional strength of the HIV-1C LTR leads to concomitantly reduced gene expression noise, consequently leading to stabler latency maintenance and rapid latency reversal. The present work offers a novel lead toward appreciating the molecular mechanisms governing HIV-1 latency.

A transcriptional cycling model recapitulates chromatin-dependent features of noisy inducible transcription

Activation of gene expression in response to environmental cues results in substantial phenotypic heterogeneity between cells that can impact a wide range of outcomes including differentiation, viral activation, and drug resistance. An important source of gene expression noise is transcriptional bursting, or the process by which transcripts are produced during infrequent bursts of promoter activity. Chromatin accessibility impacts transcriptional bursting by regulating the assembly of transcription factor and polymerase complexes on promoters, suggesting that the effect of an activating signal on transcriptional noise will depend on the initial chromatin state at the promoter. To explore this possibility, we simulated transcriptional activation using a transcriptional cycling model with three promoter states that represent chromatin remodeling, polymerase binding and pause release. We initiated this model over a large parameter range representing target genes with different chromatin environments, and found that, upon increasing the polymerase pause release rate to activate transcription, changes in gene expression noise varied significantly across initial promoter states. This model captured phenotypic differences in activation of latent HIV viruses integrated at different chromatin locations and mediated by the transcription factor NF-κB. Activating transcription in the model via increasing one or more of the transcript production rates, as occurs following NF-κB activation, reproduced experimentally measured transcript distributions for four different latent HIV viruses, as well as the bimodal pattern of HIV protein expression that leads to a subset of reactivated virus. Importantly, the parameter ‘activation path’ differentially affected gene expression noise, and ultimately viral activation, in line with experimental observations. This work demonstrates how upstream signaling pathways can be connected to biological processes that underlie transcriptional bursting, resulting in target gene-specific noise profiles following stimulation of a single upstream pathway.

Many genes are transcribed in infrequent bursts of mRNA production through a process called transcriptional bursting, which contributes to variability in responses between cells. Heterogeneity in cell responses can have important biological impacts, such as whether a cell supports viral replication or responds to a drug, and thus there is an effort to describe this process with mathematical models to predict biological outcomes. Previous models described bursting as a transition between an “OFF” state or an “ON” state, an elegant and simple mathematical representation of complex molecular mechanisms, but one which failed to capture how upstream activation signals affected bursting. To address this, we added an additional promoter state to better reflect biological mechanisms underlying bursting. By fitting this model to variable activation of quiescent HIV infections in T cells, we showed that our model more accurately described viral expression variability across cells in response to an upstream stimulus. Our work highlights how mathematical models can be further developed to understand complex biological mechanisms and suggests ways to connect transcriptional bursting to upstream activation pathways.

T cell stimulation remodels the latently HIV-1 infected cell population by differential activation of proviral chromatin

The reservoir of latently HIV-1 infected cells is heterogeneous. To achieve an HIV-1 cure, the reservoir of activatable proviruses must be eliminated while permanently silenced proviruses may be tolerated. We have developed a method to assess the proviral nuclear microenvironment in single cells. In latently HIV-1 infected cells, a zinc finger protein tethered to the HIV-1 promoter produced a fluorescent signal as a protein of interest came in its proximity, such as the viral transactivator Tat when recruited to the nascent RNA. Tat is essential for viral replication. In these cells we assessed the proviral activation and chromatin composition. By linking Tat recruitment to proviral activity, we dissected the mechanisms of HIV-1 latency reversal and the consequences of HIV-1 production. A pulse of promoter-associated Tat was identified that contrasted to the continuous production of viral proteins. As expected, promoter H3K4me3 led to substantial expression of the provirus following T cell stimulation. However, the activation-induced cell cycle arrest and death led to a surviving cell fraction with proviruses encapsulated in repressive chromatin. Further, this cellular model was used to reveal mechanisms of action of small molecules. In a proof-of-concept study we determined the effect of modifying enhancer chromatin on HIV-1 latency reversal. Only proviruses resembling active enhancers, associated with H3K4me1 and H3K27ac and subsequentially recognized by BRD4, efficiently recruited Tat upon cell stimulation. Tat-independent HIV-1 latency reversal of unknown significance still occurred. We present a method for single cell assessment of the microenvironment of the latent HIV-1 proviruses, used here to reveal how T cell stimulation modulates the proviral activity and how the subsequent fate of the infected cell depends on the chromatin context.

Clonal Expansion of Infected CD4+ T Cells in People Living with HIV

HIV infection is not curable with current antiretroviral therapy (ART) because a small fraction of CD4+ T cells infected prior to ART initiation persists. Understanding the nature of this latent reservoir and how it is created is essential to development of potentially curative strategies. The discovery that a large fraction of the persistently infected cells in individuals on suppressive ART are members of large clones greatly changed our view of the reservoir and how it arises. Rather than being the products of infection of resting cells, as was once thought, HIV persistence is largely or entirely a consequence of infection of cells that are either expanding or are destined to expand, primarily due to antigen-driven activation. Although most of the clones carry defective proviruses, some carry intact infectious proviruses; these clones comprise the majority of the reservoir. A large majority of both the defective and the intact infectious proviruses in clones of infected cells are transcriptionally silent; however, a small fraction expresses a few copies of unspliced HIV RNA. A much smaller fraction is responsible for production of low levels of infectious virus, which can rekindle infection when ART is stopped. Further understanding of the reservoir will be needed to clarify the mechanism(s) by which provirus expression is controlled in the clones of cells that constitute the reservoir.

Differences in Transcriptional Dynamics Between T-cells and Macrophages as Determined by a Three-State Mathematical Model

HIV-1 viral transcription persists in patients despite antiretroviral treatment, potentially due to intermittent HIV-1 LTR activation. While several mathematical models have been explored in the context of LTR-protein interactions, in this work for the first time HIV-1 LTR model featuring repressed, intermediate, and activated LTR states is integrated with generation of long (env) and short (TAR) RNAs and proteins (Tat, Pr55, and p24) in T-cells and macrophages using both cell lines and infected primary cells. This type of extended modeling framework allows us to compare and contrast behavior of these two cell types. We demonstrate that they exhibit unique LTR dynamics, which ultimately results in differences in the magnitude of viral products generated. One of the distinctive features of this work is that it relies on experimental data in reaction rate computations. Two RNA transcription rates from the activated promoter states are fit by comparison of experimental data to model predictions. Fitting to the data also provides estimates for the degradation/exit rates for long and short viral RNA. Our experimentally generated data is in reasonable agreement for the T-cell as well macrophage population and gives strong evidence in support of using the proposed integrated modeling paradigm. Sensitivity analysis performed using Latin hypercube sampling method confirms robustness of the model with respect to small parameter perturbations. Finally, incorporation of a transcription inhibitor (F07#13) into the governing equations demonstrates how the model can be used to assess drug efficacy. Collectively, our model indicates transcriptional differences between latently HIV-1 infected T-cells and macrophages and provides a novel platform to study various transcriptional dynamics leading to latency or activation in numerous cell types and physiological conditions.

Cell Size-Based Decision-Making of a Viral Gene Circuit

Latently infected T cells able to reinitiate viral propagation throughout the body remain a major barrier to curing HIV. Distinguishing between latently infected cells and uninfected cells will advance efforts for viral eradication. HIV decision-making between latency and active replication is stochastic, and drug cocktails that increase bursts of viral gene expression enhance reactivation from latency. Here, we show that a larger host-cell size provides a natural cellular mechanism for enhancing burst size of viral expression and is necessary to destabilize the latent state and bias viral decision-making. Latently infected Jurkat and primary CD4+ T cells reactivate exclusively in larger activated cells, while smaller cells remain silent. In addition, reactivation is cell-cycle dependent and can be modulated with cell-cycle-arresting compounds. Cell size and cell-cycle dependent decision-making of viral circuits may guide stochastic design strategies and applications in synthetic biology and may provide important determinants to advance diagnostics and therapies.

Bohn-Wippert et al. investigate reactivation of T cells latently infected with HIV. They discover that only larger cells exit latency, while smaller cells remain silent. Viral expression bursts are cell size and cell-cycle dependent, presenting dynamic cell states, capable of active control, as sources of viral fate determination.

Kinetics of HTLV-1 reactivation from latency quantified by single-molecule RNA FISH and stochastic modelling

The human T cell leukemia virus HTLV-1 establishes a persistent infection in vivo in which the viral sense-strand transcription is usually silent at a given time in each cell. However, cellular stress responses trigger the reactivation of HTLV-1, enabling the virus to transmit to a new host cell. Using single-molecule RNA FISH, we measured the kinetics of the HTLV-1 transcriptional reactivation in peripheral blood mononuclear cells (PBMCs) isolated from HTLV-1+ individuals. The abundance of the HTLV-1 sense and antisense transcripts was quantified hourly during incubation of the HTLV-1-infected PBMCs ex vivo. We found that, in each cell, the sense-strand transcription occurs in two distinct phases: the initial low-rate transcription is followed by a phase of rapid transcription. The onset of transcription peaked between 1 and 3 hours after the start of in vitro incubation. The variance in the transcription intensity was similar in polyclonal HTLV-1+ PBMCs (with tens of thousands of distinct provirus insertion sites), and in samples with a single dominant HTLV-1+ clone. A stochastic simulation model was developed to estimate the parameters of HTLV-1 proviral transcription kinetics. In PBMCs from a leukemic subject with one dominant T-cell clone, the model indicated that the average duration of HTLV-1 sense-strand activation by Tax (i.e. the rapid transcription) was less than one hour. HTLV-1 antisense transcription was stable during reactivation of the sense-strand. The antisense transcript HBZ was produced at an average rate of ~0.1 molecules per hour per HTLV-1+ cell; however, between 20% and 70% of HTLV-1-infected cells were HBZ-negative at a given time, the percentage depending on the individual subject. HTLV-1-infected cells are exposed to a range of stresses when they are drawn from the host, which initiate the viral reactivation. We conclude that whereas antisense-strand transcription is stable throughout the stress response, the HTLV-1 sense-strand reactivation is highly heterogeneous and occurs in short, self-terminating bursts.

Targeting HIV-1 proviral transcription

Despite the success of antiretroviral therapies, there is no cure for HIV-1 infection due to the establishment of a long-lived latent reservoir that fuels viral rebound upon treatment interruption. 'Shock-and-kill' strategies to diminish the latent reservoir have had modest impact on the reservoir leading to considerations of alternative approaches to target HIV-1 proviruses. This review explores approaches to target HIV-1 transcription as a way to block the provirus expression.

Post-Transcriptional Noise Control

Recent evidence indicates that transcriptional bursts are intrinsically amplified by messenger RNA cytoplasmic processing to generate large stochastic fluctuations in protein levels. These fluctuations can be exploited by cells to enable probabilistic bet-hedging decisions. But large fluctuations in gene expression can also destabilize cell-fate commitment. Thus, it is unclear if cells temporally switch from high to low noise, and what mechanisms enable this switch. Here, the discovery of a post-transcriptional mechanism that attenuates noise in HIV is reviewed. Early in its life cycle, HIV amplifies transcriptional fluctuations to probabilistically select alternate fates, whereas at late times, HIV utilizes a post-transcriptional feedback mechanism to commit to a specific fate. Reanalyzing various reported post-transcriptional negative feedback architectures reveals that they attenuate noise more efficiently than classic transcriptional autorepression, leading to the derivation of an assay to detect post-transcriptional motifs. It is hypothesized that coupling transcriptional and post-transcriptional autoregulation enables efficient temporal noise control to benefit developmental bet-hedging decisions.

Probabilistic control of HIV latency and transactivation by the Tat gene circuit

The reservoir of HIV latently infected cells is the major obstacle for eradication of HIV infection. The "shock-and-kill" strategy proposed earlier aims to reduce the reservoir by activating cells out of latency. While the intracellular HIV Tat gene circuit is known to play important roles in controlling latency and its transactivation in HIV-infected cells, the detailed control mechanisms are not well understood. Here we study the mechanism of probabilistic control of the latent and the transactivated cell phenotypes of HIV-infected cells. We reconstructed the probability landscape, which is the probability distribution of the Tat gene circuit states, by directly computing the exact solution of the underlying chemical master equation. Results show that the Tat circuit exhibits a clear bimodal probability landscape (i.e., there are two distinct probability peaks, one associated with the latent cell phenotype and the other with the transactivated cell phenotype). We explore potential modifications to reactions in the Tat gene circuit for more effective transactivation of latent cells (i.e., the shock-and-kill strategy). Our results suggest that enhancing Tat acetylation can dramatically increase Tat and viral production, while increasing the Tat-transactivation response binding affinity can transactivate latent cells more rapidly than other manipulations. Our results further explored the "block and lock" strategy toward a functional cure for HIV. Overall, our study demonstrates a general approach toward discovery of effective therapeutic strategies and druggable targets by examining control mechanisms of cell phenotype switching via exactly computed probability landscapes of reaction networks.

A Post-Transcriptional Feedback Mechanism for Noise Suppression and Fate Stabilization

Diverse biological systems utilize fluctuations ("noise") in gene expression to drive lineage-commitment decisions. However, once a commitment is made, noise becomes detrimental to reliable function, and the mechanisms enabling post-commitment noise suppression are unclear. Here, we find that architectural constraints on noise suppression are overcome to stabilize fate commitment. Using single-molecule and time-lapse imaging, we find that-after a noise-driven event-human immunodeficiency virus (HIV) strongly attenuates expression noise through a non-transcriptional negative-feedback circuit. Feedback is established through a serial cascade of post-transcriptional splicing, whereby proteins generated from spliced mRNAs auto-deplete their own precursor unspliced mRNAs. Strikingly, this auto-depletion circuitry minimizes noise to stabilize HIV's commitment decision, and a noise-suppression molecule promotes stabilization. This feedback mechanism for noise suppression suggests a functional role for delayed splicing in other systems and may represent a generalizable architecture of diverse homeostatic signaling circuits.

Relatively slow stochastic gene-state switching in the presence of positive feedback significantly broadens the region of bimodality through stabilizing the uninduced phenotypic state

Within an isogenic population, even in the same extracellular environment, individual cells can exhibit various phenotypic states. The exact role of stochastic gene-state switching regulating the transition among these phenotypic states in a single cell is not fully understood, especially in the presence of positive feedback. Recent high-precision single-cell measurements showed that, at least in bacteria, switching in gene states is slow relative to the typical rates of active transcription and translation. Hence using the lac operon as an archetype, in such a region of operon-state switching, we present a fluctuating-rate model for this classical gene regulation module, incorporating the more realistic operon-state switching mechanism that was recently elucidated. We found that the positive feedback mechanism induces bistability (referred to as deterministic bistability), and that the parameter range for its occurrence is significantly broadened by stochastic operon-state switching. We further show that in the absence of positive feedback, operon-state switching must be extremely slow to trigger bistability by itself. However, in the presence of positive feedback, which stabilizes the induced state, the relatively slow operon-state switching kinetics within the physiological region are sufficient to stabilize the uninduced state, together generating a broadened parameter region of bistability (referred to as stochastic bistability). We illustrate the opposite phenotype-transition rate dependence upon the operon-state switching rates in the two types of bistability, with the aid of a recently proposed rate formula for fluctuating-rate models. The rate formula also predicts a maximal transition rate in the intermediate region of operon-state switching, which is validated by numerical simulations in our model. Overall, our findings suggest a biological function of transcriptional “variations” among genetically identical cells, for the emergence of bistability and transition between phenotypic states.

Identifying the mechanism underlying the coexistence of multiple stable phenotypic states has been a challenging scientific problem for more than half a century, and an appropriate mathematical model at the single-cell level is also in high demand. Single-cell measurements conducted in the past ten years have shown that gene-state switching is slow relative to the typical rates of active transcription and translation; hence the recently proposed fluctuating-rate model is a good candidate for describing the single-cell dynamics. We use the classic gene regulation module of the lac operon as an archetype and build a specific fluctuating-rate model based on the recently identified operon-state switching mechanism. This model is analyzed to dissect the interplay between positive feedback and the stochastic switching of gene states in the emergence of bistability/multistablity and the transition between phenotypic states. We show that relatively slow operon-state switching stabilizes the uninduced state and that the positive feedback stabilizes the induced state. Thus, the parameter range for bistability is significantly broadened. In addition, recently proposed landscape theory and rate formula predict opposite phenotype-transition rate dependence on operon-state switching rates for the two types of bistability.