Fluctuations in Tat copy number when it counts the most: a possible mechanism to battle the HIV latency
© Konkoli and Jesorka; licensee BioMed Central Ltd. 2013
Received: 26 October 2012
Accepted: 12 February 2013
Published: 5 March 2013
The HIV-1 virus can enter a dormant state and become inactive, which reduces accessibility by antiviral drugs. We approach this latency problem from an unconventional point of view, with the focus on understanding how intrinsic chemical noise (copy number fluctuations of the Tat protein) can be used to assist the activation process of the latent virus. Several phase diagrams have been constructed in order to visualize in which regions of the parameter space noise can drive the activation process. Essential to the study is the use of a hyperbolic coordinate system, which greatly facilitates quantification of how the various reaction rate combinations shape the noise behavior of the Tat protein feedback system. We have designed a mathematical manual of how to approach the problem of activation quantitatively, and introduce the notion of an “operating point” of the virus. For both noise-free and noise-based strategies we show how operating point off-sets induce changes in the number of Tat molecules. The major result of the analysis is that for every noise-free strategy there is a noise-based strategy that requires lower dosage, but achieves the same anti-latency effect. It appears that the noise-based activation is advantageous for every operating point.
HIV-1 latency is a serious problem that prevents eradication of the virus
Human immunodeficiency virus 1 (HIV-1), initially reported nearly 30 years ago, represents a major global health problem with millions of people infected [1, 2]. One of the biggest problems with HIV-1 is that the virus can enter a dormant state and effectively “hide” from drug cocktails that are therapeutically administered. Several reviews have been written on the subject [3–13]. Latent viral reservoirs can persist for many years. Once a therapy is interrupted, the latent reservoirs remain as the source of eventually renewed infection. This behaviour has been identified as the key problem in eradicating HIV-1.
There are many reservoirs, e.g., cell types, in the body that can harbor the latent virus. CD4+ T-cells have been identified as one of the largest pools of the viral DNA. This is also one of the best characterized reservoirs . Once a T-cell is infected, it can either become activated and produce new virus particles that will infect other cells, or it can enter an inactive state. In the inactive (latent) state the transcription of the viral DNA is silenced, despite the fact that the viral RNA has been reversely transcribed and inserted into the host DNA.
The process of entry into the latent state is rather complex since it is controlled by a sizeable number of processes which need to occur nearly at the same time [8–10, 13]. There is a relatively low percentage of latently infected T-cells, being roughly one per 106, though the frequency can be lower . This also suggests that once the latent state is established it remains very stable, and spontaneous activation events are very rare. In order to activate a latent cell, the inactive cellular processes need to be reactivated, which is not likely to occur simultaneously. A large number of studies have been performed in the past with the goal to find a way of reliably activating the latent virus, as reviewed in [3–13].
The mechanisms responsible for the maintenance of the HIV-1 latency work mostly at the molecular level. Processes such as chromatin control, a shortage of host transcription factors that initiate transcription, the presence of molecules that slow down (or even block) polymerase elongation, transcriptional interference, DNA methylation, and insufficient transport of viral mRNA from the nucleus to the cytoplasm are some prominent examples. All known mechanisms have been reviewed in detail, for example in [8–10, 13].
A small number of studies focused on the role of chemical (intrinsic) noise in establishing the latent state [14–17]. A generic conclusion extracted from these studies is that noise is detrimental for the latency decision. Noise comes from fluctuations in copy numbers of the proteins that are involved in the latency control. The level of noise in the system is partially modulated by host factors (proteins that are present in healthy cells) and partially with virus specific factors (proteins encoded by viral genetic machinery). Noise makes the latency decision stochastic (unpredictable) and the literature suggests that the noise driven inactivation occurs spontaneously.
In the present study the focus is placed on understanding how noise can be manipulated to reverse the latency decision, i.e., to drive the activation process. This idea is fully in line with several previous findings and suggestions in the literature. For example, it has been clearly appreciated that noise plays a role in various cellular decision making processes [18, 19], but the idea to use noise to steer cell fate decisions has never been seriously explored. It has been argued, merely on a general basis, that one should strive to control noise better in order to steer cell fate decisions . In the context of the HIV activation it has been found that altering noise of the HIV genetic machinery can bias virus decision making towards productive replication or latency .
In our study the focus is on finding ways to manipulate noise in order to specifically facilitate the productive activation, which could become the foundation of medical treatment strategies in the future. To our knowledge, this way of approaching the eradication of the latent virus has not yet been much discussed in the literature.
The first goal of this study is to identify regions in parameter space where noise greatly influences the dynamics. We investigate fluctuations of copy numbers of several key proteins that are involved in the maintenance of the latent state of infected cells. Based on these insights, the second goal is to identify suitable medical strategies that can be used to harvest noise in order to achieve more efficient activation. These intuitive considerations will be addressed in a quantitative way by means of a rigorous mathematics.
The Tat feedback loop is a main source of noise driving the activation of the latent virus
The Tat feedback loop has been identified as an important part of the HIV gene expression machinery, reviewed in [3–5, 7] and others. The main biological function of the Tat loop is to accelerate viral RNA production which further increases the number of viral particles in the infected cell. The presence of Tat molecules in the cell increases the transcription rate by roughly two orders of magnitude. This normally leads to cell lysis and continued infection of other cells. A lack of Tat molecules in latently infected cells has been identified as a strong barrier to activation of the latent cell [8–10, 13].
Interestingly, the loop appears to be important for noise driven entry into the latent state [14–17]. If viral DNA integrates into regions of high basal transcription then noise does not play such a big role. However, if it integrates in regions of low basal transcription, then the positive feedback loop can amplify noise effectively and produce bursts of activity. In such a way the fate of the cell is determined in a stochastic manner, as shown in [14, 15]. Such behavior splits the typical low transcription isogenic cell population into two groups. In the first group the Tat feedback loop produces new Tat particles (the “on” state). In the second group the Tat feedback loop is inactivated (the “off” state). This phenomenon was termed the phenotypic bifurcation (PheB). A series of carefully designed experiments were performed  to show that such behavior is indeed a product of intracellular noise, and is independent from extracellular perturbations. The possibility of a spontaneous latency decision caused by uncontrollable fluctuations of the Tat copy number has been clearly demonstrated.
Current anti-latency treatments need to be improved
Even though several anti-latency drugs have been investigated, it is not clear how to administer them efficiently. For example, there is no consensus whether they should be administered aggressively (all at once), by repeated injection of smaller amounts, or constantly administered over a longer time period [4, 5, 21, 22].
It is generally agreed upon that due to potential side effects of the treatment (e.g. toxicity), minimal dosages are preferable. A convincing argument for administering Minimal dosages is that too large amounts of the activation drug might release the virus beyond control, such that the usual highly active antiretroviral therapy (HAART) cannot contain it any longer . Moreover, there are viral reservoirs, e.g., the brain, that are not easily accessible by HAART. In such organs it is very important to activate the virus gradually.
Clinical studies show that a global T-cell activation cannot eradicate the virus. Instead, it causes unwanted side effects . An additional problem is that many of the activating agents are generic to a wide range of gene regulation processes; administering them is expected to show severe toxicity in the host. Thus the reduction of the dosage of anti-latency drugs appears desirable. In this context, noise-driven activation could be highly beneficial.
Intrinsic fluctuations in protein copy numbers could be used to achieve more efficient activation
A therapy can be envisioned where several types of drugs work in synergy with the specific aim of moving the operating point of the virus into a noisy region, where the frequency of spontaneous virus activation would increase. Such therapy could be sustained for a longer time since, presumably, lower dosages of anti-latency drugs could be used. Targeting fluctuations in Tat copy number is very natural in this context, since the Tat protein is an essential part of the gene transcription machinery which produces new viral particles.
Fluctuations in Tat copy number have the potential to drive the activation. The up to date understanding of the feedback loop is that Tat controls transcriptional elongation rather than initiation. However, it seems that these two processes cannot be clearly separated, as it has been shown that an exogenous injection of Tat can activate the latent cell [14, 23–27]. We note that a study on mice  showed no obvious side effects.
It is somewhat surprising that exogenous administration of Tat can in fact have a positive effect on the activation of infected cells. The lack of Tat molecules is just one of the many barriers to activation. If the transcription of the viral DNA is also blocked by other mechanisms that are not controlled by the Tat protein, the injection of exogenous Tat should not speed up the transcription process automatically. In fact, there are also cases where Tat by itself cannot reactivate the virus .
This supports a notion that in latently infected cells the transcription system is in a rather labile state, and that the factors affecting the transcription process do not work strictly in a binary on-off fashion. Such a system could be activated by spontaneous fluctuations in Tat copy number. Based on this insight, in the next section we will construct a mathematical machinery to identify useful strategies to achieve noise-assisted activation of the virus.
A mathematical manual for the design of noise-based activation strategies
The central idea in the subsequent discussions is the concept of the “operating point” of the virus. The term will refer to a particular choice of the parameters that define the dynamics of the system, i.e., the virus and the cell it infects. In fact, one of the difficulties in combating HIV is that the integration of the virus into different sites gives rise to proviruses that have variable gene expression properties. It is conceivable that this results in different operating points of the virus. Therefore it is not immediately clear how a drug designed to achieve noise-based activation at a specific operation point would work consistently at other operating points. It is important to address this concern.
The manual we envision must contain a description of how (i) each anti-viral treatment affects an operating point, (ii) how the amounts of administered agents affect the magnitude of the off-set, and (iii) which anti-latency effect the induced off-set has. If the operating point of the virus is moved to regions where the latent state is less stable, the effect of antiviral drugs might be enhanced, which implies reduced dosage and shortened treatment. Below, each of these key considerations is mathematically formalized.
Mathematical description of therapies and dosage
Assumption I. The main function of each anti-latency agent is to alter the operating point of a virus: Assume that the virus is at an operating point , being a list of relevant parameters that describe the reaction system (the latent cell). If a certain amount of anti-latency agents is administered then this will off-set the operating point of the virus by and move it to a new operating point . Thus mathematically, with every activation strategy one needs to associate the related displacement of the operating point it induces. It is possible that a treatment generates always the same displacement regardless of the operating point of the virus, but this does not necessarily need to be the case.
should be always satisfied, i.e., with nothing administered no operating point off-set is induced.
being the standard Cartesian norm of the vector .
To avoid working with complicated metric spaces and projection procedures, the simplest possible relationship in the form of a proportionality law will be assumed.
Probability of activation and related observables
The central quantity of interest is the probability of the activation of a latent cell in an observation interval t0<t′<t, to be denoted by P(t0,t).
where n i (t′) for i=1,2,⋯,N are copy numbers of relevant reactants (e.g. proteins) at time t′.
where Φ∗ is a constant.
The noise-free and noise-based activation concepts
Workings and effects of noise-free therapies
is the gradient computed at the operating point .
Such a treatment will be referred to as optimal noise-free.
Workings and effects of noise-based treatments
The precise value for κ depends on the character of the particle number distribution function and the confidence interval, but in here it will be taken as κ∼1 for simplicity reasons. Note that the effects of noise can be shut-off by taking κ=0.
Such a treatment will be referred to as optimal noise-based.
The dose reduction coefficient can be used to quantify at which operating points the noise-based treatment is useful
If the size of the vector is smaller than the size of the vector then a noise based therapy with lower dosage is possible.
will be used to quantify the comparison between the two strategies. This quantity indicates the degree of the dose reduction for the treatment associated with a vector . It does not depend on the size of the displacement vector , only on its direction relative to the gradient of the mean.
to be referred to as the dose reduction coefficient. Ideally, but the opposite is perfectly possible. When noise-based therapy simply would not work as well as its noise-free counterpart.
The space of operating points where a noise-based activation strategy can work
Now the mathematical manual will be used to find regions in the space of operating points of the latent virus where a noise-based activation is effective. To do this we focus on the Tat feedback loop. It is essential to model the noise of the Tat feedback loop, and to choose a relevant observable. The loop itself will be modelled in the simplest possible way as shown in the next subsection. The features of the system we wish to describe and the related observables are discussed subsequently.
The resistor model of HIV latency
The “resistor model” of the HIV dormancy control has been suggested to explain how the lack of the Tat molecules maintains (stabilizes) the latent state [14, 15]. This simplified description of the transcription machinery will be used for the theoretical analysis. The biochemical model  consists of entities (or particle types) TatA and TatD which denote the acetylated and deacetylated form of the Tat molecule respectively. To simply the notation TatA and TatD will be abbreviated as A and D respectively.
It is clear that the model defined above neglects many biochemical details and is far from complete. For example, the mechanisms listed above are just a subset of all modifications of the Tat protein that occur in the cell . Furthermore, a weak basal expression of TatD is continuously occurring  but this process will be neglected in the same way as in . The model does not discriminate between mRNAs and proteins, either.
A class of noise-based activation strategies
This condition is a fundamental property of the system. It defines the balance between the processes that produce and destroy TatA molecules. The system is stable if the production is slower than the destruction.
The key behavior of the resistor model we wish to exploit is related to pulses of activity. We investigate a hypothetical situation where Tat molecules are administered weakly at a constant rate in combination with an anti-viral drug and one or more operation point off-setting agents. The injection of Tat should trigger repeated pulses of activity. The anti-viral drug should reduce the copy number of viruses, while the off-setting agents move the operating point of the virus towards a nosier regime. The role of the off-setting agents is to boost fluctuations in the Tat expression levels so that the dosage of the anti-viral drug can be reduced.
Two pulse characteristics as key observables
We now apply the mathematical machinery to compute the dose reduction coefficient for a wide range of reaction rate parameters (operating points) of the resistor model. It is necessary to identify an observable upon which the activation probability strongly depends. It will be assumed that the determining observable that governs the activation probability is the amount of viral mRNAs (proteins) produced by the feedback loop for the duration of a single pulse triggered by a single injection of Tat. If there are many viral particles that are ready for packaging and transport, the latent cell should activate.
is the pulse surface (PS), i.e., the surface under the n A (t) curve (see Figure 2). This scenario will be referred to as the pulse surface-dependent threshold scenario, and when appropriate quantities computed in this context will be labelled by PS as in the example above.
There is no explicit experimental evidence that the surface is the one determining factor that most strongly influences activation. However, this idea is supported by several experimental studies. For example, in the experiment published in  it was argued that the duration of the activity pulse is likely to be important. The strength of the feedback loop regulates the duration of activity bursts. For less stable states the duration of such bursts increases dramatically, until they become equally long, and eventually longer than the life-time of the cell. Furthermore, in  it was experimentally shown that a stronger feedback implies more likely activation. Stronger feedback should correlate with the size of the surface . Thus the first observable that will be used to implement the mathematical manual is the area .
where n A (t) is the number of A at time t (see Figure 2). The quantities computed in this context will be labeled by PH as above.
While ΦPH will be primarily used for illustrative purposes to gain qualitative understanding, the concentration-dependent threshold scenario does have practical significance. One expects that the number of Tat molecules in the cell correlates with the number of other viral proteins since their expression is encoded in the same gene. In that sense, high copy numbers of Tat protein should correlate with large activation probability. There are also some experimental indications in favour of the scenario.
Although there is no direct evidence that Tat on its own ensures activation, there is experimental evidence that exogenous injection of Tat can activate the latent cell . Interestingly, there is also a report of the opposite . Moreover, it has been suggested in several publications that Tat operates in a threshold dependent manner [6, 8, 31, 32], but there is no direct experimental evidence for that. Interestingly, it is known that the Rev protein which controls export of viral mRNAs operates in a threshold dependent manner, and Tat drives the concentration of Rev above the threshold (See  for a discussion).
Computation of the dose reduction coefficients for the concentration-dependent threshold scenario (the pulse height scenario)
The functions and will be defined later on.
Computation of the dose reduction coefficients for the pulse-surface dependent threshold scenario
where teq denotes the average equilibration time of the system. This quantity has to be a well-defined function of the operating point .
where is a typical fluctuation of the number of A at t=tmax and is the size of a typical fluctuation of the pulse duration. To avoid the mathematical complexity involved in computing , it will be assumed that this quantity does not fluctuate.
Thus, the required key quantities are the average number of Tat molecules at this peak, its standard deviation , and the typical time of the pulse duration scaled with the expression rate .
The mathematical description of noise: An overview of computing the mean and the standard deviation
Of particular interest will be the values of μ A (t) and σ A (t) at . The equations of motion for the factorial moments needed to compute these quantities are discussed in the next section.
The equation system for factorial moments
In principle, such equation system forms an infinite hierarchy which for the present system decouples automatically. It is known that when all propensity functions in the stochastic formulation are linear with respect to the population count the computation of statistical moments is simple. For example, the equations for ρ1,0 and ρ0,1 form a closed system, and likewise the equations for ρ2,0, ρ0,2, and ρ1,1. It can be seen that the equations for ρx,y with x+y≤ξ form a closed system for every ξ=1,2,3,⋯. This analysis implies that the solution to a particular equation system is exact. Also, since the equation system for factorial moments up to order ξ=2 is exact, the values for the mean and the variance will be exact.
Direct numerical integration should be avoided
In order to solve the equations it is necessary to integrate them from some initial condition from time t=0 until t=tmax. There are several reasons why a direct numerical integration should be avoided. First, the computational cost of a single time integration scales linearly with the inverse of a typical time step size used. A direct numerical integration is highly impractical, since to construct the phase diagrams, the noise measures have to be computed at many points in the reaction rate space.
Second, there will be a need to compute derivatives of the mean and the standard deviation with respect to the reaction rates. If obtained numerically, his would multiply the computational effort by a factor two or more, depending on which technique is used to numerically calculate the derivatives.
where the functions f μ and f σ are detailed in the next section. When convenient the subscript “max” will be omitted. To simplify notation we will use only μ and σ instead of μmax and σmax.
Details of computing the mean and the standard deviation
In this section it will be shown how to obtain the equations of motion for factorial moments and the functional forms for the mean and the standard deviation.
A procedure for avoiding numerical integration
The variable φ A D (t) is a generalization of the noise strength concept for a pair of particle types.
This can be used to obtain approximations for the noise strengths at t=tmax.
This can be seen from the fact that for the Poisson initial condition and which after using the noise strength definitions (66) and (67) automatically leads to the above result. Accordingly, given that the particle number distribution function is Poisson-like at t≈0, one can expect that φ A (t) and φ D (t) vary within a finite interval.
A single cell can be infected by more than one virus particle. The number of viral particles that infect a given cell (the multiplicity of infection) varies randomly and is usually Poisson-distributed. Accordingly, in what follows it will be assumed that the initial particle number distribution function is Poisson-like.
This is the approximation that will be used to compute noise strengths at the time instance where the number of the acetylated Tat particles reaches maximum.
There is no a priori reason why the noise strengths should vary linearly with time. We have inspected several curves where noise strengths were computed numerically to see whether the time dependence is linear. Interestingly, while the time dependence is not strictly linear it seems that the approximation used is qualitatively correct. We performed more rigorous tests of such an approximation by comparing it with the results of a numerical integration for wide range of parameters and found reasonable agreement.
The equation system for noise strengths
It will be shown later that the terms proportional to a(t) and d(t) also drop out in the asymptotic limit when the stability condition is satisfied.
Locating the peak region
In order to see what happens when relatively few deacetylated Tat molecules are injected the equation above will be used with d0=1.
Computing the asymptotic noise strengths
Note that there are no restrictions on the reaction rates, only that they are positive real numbers. This implies that in this model the ratio of the number of acetylated and deacetylated Tat molecules approaches a constant value regardless on whether the system is stable or not.
When β δ−α k>0, z(t) grows exponentially fast, implying that the means approach zero. This condition is fully equivalent to the stability condition (32). On the other hand, when β δ−α k 0, z(t) approaches zero exponentially fast, which implies that the average copy numbers become infinite. In both cases the ratio of the means becomes constant and is given by u∗.
where a(tmax) is given in (88) and is obtained by solving (93-94). All algebraic equations have been solved by using the Mathematica package (not shown).
A hyperbolic coordinate system facilitates understanding
In fact, this set of coordinates is very intuitive which can be appreciated by analysing the meaning of the inverse transformation.
is particularly important for several reasons. First, it measures how intensive the transcription process (the production of acetylated Tat molecules) is. This coordinate also measures how far in the ility region the operating point of the virus is placed. For example, ϵ=0 at the border of the stab stability region when α k=β δ. Furthermore, ϵ≈1 when the stability condition is strongly satisfied, i.e. when β δ≫k α. Thus, in the stable region ϵ attains values in the interval between zero and one.
This variable could potentially be used to experimentally quantify the degree of the latency for a given cell. For example, the literature suggests that for a latent virus its operating point has to lie in the region of stability where ϵ>0. Otherwise, for ϵ 0, an infected cell would lyse relatively fast. Thus, it seems necessary that a treatment meant to activate a latent virus has to move its operating point into the unstable region where ϵ 0 or at least sufficiently close to the instability boundary ϵ=0.
The symbol “R” emphasizes the resistor model operating point.
Regions of parameter space where noise is large
A case study: A strategy for improving noise-free histone deacetylase inhibitor treatment by moving the operating point of the resistor system to a noisier region
A typical anti-latency strategy in the HIV therapy context is to increase the transcription rate k, for example in treatments based on the use of histone deacetylase inhibitors . These molecules open up the chromatin environment such that the transcription factors needed for viral expression can attach easier to their respective binging sites. Several molecules have been suggested as drug candidates. Valproic acid, trichostatin A, vorinostat (SAHA), and many more have been reviewed in , and references therein. Some of these molecules are rather toxic. In the following we investigate how toxicity could be reduced by lowering their dosage at a typical operation point, e.g. . In order to analyse the effects of various treatments it is useful to explicitly compute the gradients of the mean and the standard deviation at this operating point. We will focus on the conceptually simpler concentration-dependent threshold scenario.
Different changes in the reaction rates have to be used in order to generate these two vastly different offsets. For example, let us investigate which changes in the transcription rate k are needed to achieve the changes in (108) and (109).
By comparing (117) and (120) one sees that the dose of the primary drug can be reduced roughly twenty times (assuming the linear relationship between the dose change and the transcription rate change). The price one has to pay is that additional drugs have to be administered which reduce β and increase δ (note that Δ ϵ is negative). It is not unrealistic that this can be eventually verified experimentally, e.g. in the context of the kinetics experiments on LTR [14, 15].
Even if the optimal noise-free treatment would be applied in this case, noise-based treatment would still result in dose reduction. An immediate use of Eq. (27) shows that .
Note that since the SAHA therapy is not optimal. As the optimal noise free (gradient based) therapy is more efficient than the SAHA treatment, the corresponding noise based therapy would yield a smaller dose reduction.
This case study shows that it is sufficient to inspect the value of the dosage reduction coefficient in order to test whether a noise-based treatment would be beneficial at an operating point .
Regions of parameter space where noise-driven activation is possible
A striking result of this visual analysis is that regions where could not be found within the hyperplanes. For both scenarios, we computed the dose reduction coefficient for 10000 randomly selected operating points, and found always Θ 1. Moreover, for the PH scenario, at all points the angle between and is always in the ±π/2 interval, i.e., the vectors point in roughly the same direction. It appears that the noise-based activation is advantageous everywhere, regardless which value for κ is chosen. The degree of the dose reduction is clearly controlled by κ. Note that the related graphs for (ΘPH), (ΘPS) are roughly the same, which implies that the conclusion is generic. There are some differences that are worth discussing.
The contours in panels (a) and (d) differ in the middle region, and likewise for panels (c) and (f). There are even some regions where a relatively large dose reduction is possible. Figure 4a indicates that an efficient dosage reduction can be achieved in the regions around u≈0 where ϵ is either very large or very small. This does not hold for panel (d). However, both (a) and (d) panels suggest that the noise based therapy is advantageous in the region where the absolute value of u is very large.
Figures 4b and 4e show that for fixed u the largest dose reduction can be reached for very small values of ϵ and very large values of v. Figures 4c and 4f indicate that the noise-based activation can be useful in the regions where v is large and where either u→∞ or u→−∞. The panels do not agree in the middle region.
The linear theory yields qualitative predictions when operating point off-sets are not small
Several contour lines are shown: (NF1) μmax=μ1, (NF2) μmax=μ1+Δ N, (NB1) μmax+σmax=μ1+σ1, and (NB2) μmax+σmax=μ1+σ1+Δ N where μ1=0.05378, μ1+σ1=0.064536, and Δ N=0.010756.
The long black arrow denotes the linearized noise-free off-set that needs to be induced in order to shift the mean from μmax to μmax+Δ N. The long grey arrow is the shortest non-infinitesimal noise-free off-set that needs to be induced to move from the operating point to the NF2 contour line. The large grey circle (only in part visible) starts touching NF2 exactly at the point where the arrow meets the contour line. The short black arrow denotes the noise-based off-set that needs to be induced in order to shift μ+σ by Δ N. The small grey circle denotes the minimal circle that touches NB2.
The length ratio between the black vectors is specified exactly by the dose reduction coefficient . The grey arrow (the exact off-set) does not match its black counterpart neither in length nor direction. Clearly, the size of the noise-free offset needed to achieve the shift in μ is too large for the linearized version of the problem. This approximation is more satisfactory for the noise-based treatment since the related noise-based off-set is still relatively small. The linearized theory provides qualitative results for non-infinitesimal off-sets and that can be used in practise to guide future experiments.
A few comments on the robustness of the results with regard to the model extension
The model used in here cannot capture all phenomena that might be important. The first feature of the model that can be clearly improved is its complexity. For example, we made no distinction between mRNAs and proteins. It has been found that after integration into the genome the HIV promoter makes mRNAs in random bursts of transcriptional activity [37, 38]. Moreover, each mRNA makes proteins in translational bursts. This suggests that there might be other sources of noise in the system that were not considered in this work.
Also, a rather phenomenological model was used for the probability of activation and the related observable. There is clearly a need for refinement. It is possible that the probability of the virus activation depends on other features of the system, e.g., the amount of relative fluctuations. This is rather speculative, based upon indications in the literature that Tat operates also outside the Tat feedback loop. Tat is a multi-functional protein that is involved in other intracellular processes [32, 39, 40] and there is certainly the possibility that it influences cell homeostasis in a more complex manner than discussed here. For example, in addition to acetylation, Tat undergoes several other post-translational modifications and interacts with several other proteins. It is currently not entirely clear how this might affect gene expression noise on one hand, and the activation probability on the other. Post-translational modifications or other interactions could buffer, or propagate noise in levels of Tat into noise in gene expression. The question remains whether such effects could be beneficial in bringing the operating point closer to the stability boundary. Out of all parameters in the model, such effects will very likely exert strong influence on the decay constant δ of the deacetylated Tat and cause it to fluctuate in time. One could speculate that in such a situation both u in v will start fluctuating, though u, being essentially the natural logarithm of the ratio between β and δ, will fluctuate much less. This implies that v fluctuates towards smaller values, away from the noise-rich regions, while u and ϵ are kept essentially constant (e.g. see the phase diagrams). Such changes in v could reduce the activation probability, and need to be further investigated.
We now briefly discuss possible effects of chromatin on the operating point of a virus. The transcription rate k will likely be influenced the most. This implies that instead of analysing the effect of a drug on a single operating point one should investigated a set of operating points. Such points should be distributed around the region where v≈v R , u≈u R and where ϵ∈[0,1]. Accordingly, a noise-based treatment designed to move the operating point into a noisier region will work uniformly on all points in this set. This would always result in a relatively large dose reduction; Figure 4, panels (a), (b), (d) and (e).
We have designed a generic mathematical manual of how to approach the problem of the HIV latency in a quantitative manner, accompanied by an example of how to use this manual. We formalised several concepts that are vaguely defined or understood only intuitively. The first key concept is the notion of an “operating point” of the virus and how it is affected by the action of a drug. The second key concepts of the mathematical formalism is the notion of a particular observable that strongly affects the activation probability. We suggested a mathematical way of describing how each therapy affects an operating point and how this in turn influences the observable that controls the activation probability. The third key ingredient is a dose reduction coefficient, which can be used to quantitatively compare various therapies.
We have suggested and investigated two rather general strategies for the virus activation, the noise-free and the noise-based strategy. In the first approach, anti-latency agents are administered such that the activation happens with almost absolute certainty. In order to achieve such certainty, possibly unreasonable quantities of drugs need to be administered. In the second approach, drugs are administered in such a way that the activation is less certain but still happens with a relatively large probability. The idea behind the noise-based strategy is to reduce the quantity of drugs that need to be administered in order to achieve activation.
The mathematical manual is rather generic. To demonstrate how the mathematical manual can be used, we have focused on the simplest possible model of the Tat protein feedback loop, the most important part of the the HIV latency control. Based on the structure of the loop we suggested a class of noise-based activation strategies. We envision such an activation strategy in a procedure where one constantly supplies exogenous Tat at a very small rate, and adds a combination of anti-latency drugs that would off-set the operating point of the virus towards noisier regions.
In this context we considered two observables and could compute the dose reduction coefficient for both cases to answer the fundamental question: for which operating points the noise-base therapy is advantageous over the noise-free therapy in the sense of possible dosage reduction? This addresses the practical problem of reducing effects of toxicity during the anti-latent treatment. Three phase diagrams were constructed to explain what controls the noise, and how this control can be used to battle latency. Our analysis of the phase diagrams indicates that the noise-based therapy is always advantageous, no matter which operating point the virus adopts. This results holds regardless of which observable is targeted. This is the major result of our analysis.
The mathematical manual is currently based on rather qualitative assumptions, but it is very generic and can be easily extended and refined. For example, we have used the assumption of small operating point offsets in order to linearize the theory. One can easily consider non-infinitesimal off-sets, as demonstrated in the example. We discussed several possible extensions which are left for future work.
In summary, we suggest an activation principle where intrinsic noise is considered a feature benefiting treatment. We showed that such strategy should be efficient for any latent cell.
This work was partly funded by the Swedish Research Council (VR).
- Greene WC: A history of AIDS: looking back to see ahead. Eur J Immunol. 2007, 37: S94-S102. 10.1002/eji.200737441.View ArticlePubMed
- Le Douce V, Janossy A, Hallay H, Ali S, Riclet R, Rohr O, Schwartz C: Achieving a cure for HIV infection: do we have reasons to be optimistic?. J Antimicrob Chemother. 2012, 10.1093/jac/dkr599.
- Coiras M, Lopez-Huertas MR, Perez-Olmeda M, Alcami J: Understanding HIV-1 latency provides clues for the eradication of long-term reservoirs. Nat Rev Microbiol. 2009, 7: 798-812. 10.1038/nrmicro2223.View ArticlePubMed
- Contreras X, Lenasi T, Peterlin BM: HIV latency: present knowledge and future directions. Future Virol. 2006, 1: 733-745. 10.2217/174607220.127.116.113.PubMed CentralView ArticlePubMed
- Dahl V, Josefsson L, Palmer S: HIV reservoirs, latency, and reactivation: prospects for eradication. Antiviral Res. 2010, 85: 286-294. 10.1016/j.antiviral.2009.09.016.View ArticlePubMed
- Graci JD, Colacino JM, Peltz SW, Dougherty JP, Gu Z: HIV type-1 latency: targeted induction of proviral reservoirs. Antivir Chem Chemother. 2009, 19: 177-87.View ArticlePubMed
- Karn J: The molecular biology of HIV latency: breaking and restoring the Tat-dependent transcriptional circuit. Curr Opin HIV AIDS. 2011, 6: 4-11. 10.1097/COH.0b013e328340ffbb.PubMed CentralView ArticlePubMed
- Mbonye U, Karn J: Control of HIV latency by epigenetic and non-epigenetic mechanisms. Curr Hiv Res. 2011, 9: 554-567. 10.2174/157016211798998736.PubMed CentralView ArticlePubMed
- Schiralli Lester GM, Henderson AJ: Mechanisms of HIV transcriptional regulation and their contribution to latency. Mol Biol Int. 2012, 2012: 614120-PubMed CentralView ArticlePubMed
- Siliciano RF, Greene WC: HIV latency. Cold Spring Harb Perspect Med. 2011, 1: a007096-10.1101/cshperspect.a007096.PubMed CentralView ArticlePubMed
- Soudeyns H, Pantaleo G: The moving target: mechanisms of HIV persistence during primary infection. Immunol Today. 1999, 20: 446-450. 10.1016/S0167-5699(99)01504-2.View ArticlePubMed
- Williams SA, Greene WC: Regulation of HIV-1 latency by T-cell activation. Cytokine. 2007, 39: 63-74. 10.1016/j.cyto.2007.05.017.PubMed CentralView ArticlePubMed
- Richman DD, Margolis DM, Delaney M, Greene WC, Hazuda D, Pomerantz RJ: The challenge of finding a cure for HIV infection. Science. 2009, 323: 1304-1307. 10.1126/science.1165706.View ArticlePubMed
- Weinberger LS, Burnett JC, Toettcher JE, Arkin AP, Schaffer DV: Stochastic gene expression in a lentiviral positive-feedback loop: HIV-1 Tat fluctuations drive phenotypic diversity. Cell. 2005, 122: 169-182. 10.1016/j.cell.2005.06.006.View ArticlePubMed
- Weinberger LS, Shenk T: An HIV feedback resistor: auto-regulatory circuit deactivator and noise buffer. Plos Biol. 2007, 5: 67-81. 10.1371/journal.pbio.0050067.View Article
- Burnett JC, Miller-Jensen K, Shah PS, Arkin AP, Schaffer DV: Control of stochastic gene expression by host factors at the HIV promoter. Plos Pathog. 2009, 5: e1000260-10.1371/journal.ppat.1000260.PubMed CentralView ArticlePubMed
- Razooky BS, Weinberger LS: Mapping the architecture of the HIV-1 Tat circuit: a decision-making circuit that lacks bistability and exploits stochastic noise. Methods. 2011, 53: 68-77. 10.1016/j.ymeth.2010.12.006.PubMed CentralView ArticlePubMed
- Raj A, van Oudenaarden A: Nature, Nurture, or Chance: stochastic gene expression and its consequences. Cell. 2008, 135: 216-226. 10.1016/j.cell.2008.09.050.PubMed CentralView ArticlePubMed
- Balazsi G, van Oudenaarden A, Collins JJ: Cellular decision making and biological noise: from microbes to mammals. Cell. 2011, 144: 910-925. 10.1016/j.cell.2011.01.030.PubMed CentralView ArticlePubMed
- Singh A, Weinberger LS: Stochastic gene expression as a molecular switch for viral latency. Curr Opin Microbiol. 2009, 12: 460-466. 10.1016/j.mib.2009.06.016.PubMed CentralView ArticlePubMed
- Margolis DM: Eradication therapies for HIV infection: time to begin again. AIDS Res Hum Retroviruses. 2011, 27: 347-353. 10.1089/aid.2011.0017.PubMed CentralView ArticlePubMed
- Shan L, Deng K, Shroff NS, Durand CM, Rabi SA, Yang HC, Zhang H, Margolick JB, Blankson JN, Siliciano RF: Stimulation of HIV-1-Specific cytolytic T lymphocytes facilitates elimination of latent viral reservoir after virus reactivation. Immunity. 2012, 36: 491-501. 10.1016/j.immuni.2012.01.014.PubMed CentralView ArticlePubMed
- Donahue DA, Kuhl BD, Sloan RD, Wainberg MA: The viral protein Tat can inhibit the establishment of HIV-1 latency. J Virol. 2012, 86: 3253-3263. 10.1128/JVI.06648-11.PubMed CentralView ArticlePubMed
- Gonzalez-Hernandez MJ, Swanson MD, Contreras-Galindo R, Cookinham S, King SR, Noel JRJ, Kaplan MH, Markovitz DM: Expression of human endogenous retrovirus type-K (HML-2) is activated by the Tat protein of HIV-1. J Virol. 2012, 86 (15): 7790-7805. 10.1128/JVI.07215-11.PubMed CentralView ArticlePubMed
- Lin X, Irwin D, Kanazawa S, Huang L, Romeo J, Yen TS, Peterlin BM: Transcriptional profiles of latent human immunodeficiency virus in infected individuals: effects of Tat on the host and reservoir. J Virol. 2003, 77: 8227-36. 10.1128/JVI.77.15.8227-8236.2003.PubMed CentralView ArticlePubMed
- Emiliani S, Fischle W, Ott M, Van Lint C, Amella CA, Verdin E: Mutations in the tat gene are responsible for human immunodeficiency virus type 1 postintegration latency in the U1 cell line. J Virol. 1998, 72: 1666-1670.PubMed CentralPubMed
- Bettaccini AA, Baj A, Accolla RS, Basolo F, Toniolo AQ: Proliferative activity of extracellular HIV-1 Tat protein in human epithelial cells: expression profile of pathogenetically relevant genes. Bmc Microbiol. 2005, 5: 20-10.1186/1471-2180-5-20.PubMed CentralView ArticlePubMed
- Pearson R, Kim YK, Hokello J, Lassen K, Friedman J, Tyagi M, Karn J: Epigenetic silencing of human immunodeficiency virus (HIV) transcription by formation of restrictive chromatin structures at the viral long terminal repeat drives the progressive entry of HIV into latency. J Virol. 2008, 82: 12291-12303. 10.1128/JVI.01383-08.PubMed CentralView ArticlePubMed
- Ott M, Geyer M, Zhou Q: The control of HIV transcription: keeping RNA polymerase II on track. Cell Host Microbe. 2011, 10: 426-435. 10.1016/j.chom.2011.11.002.PubMed CentralView ArticlePubMed
- Weinberger LS, Dar RD, Simpson ML: Transient-mediated fate determination in a transcriptional circuit of HIV. Nat Genet. 2008, 40: 466-470. 10.1038/ng.116.View ArticlePubMed
- Leonard JN, Schaffer DV: Computational design of antiviral RNA interference strategies that resist human immunodeficiency virus escape. J Virol. 2005, 79: 1645-1654. 10.1128/JVI.79.3.1645-1654.2005.PubMed CentralView ArticlePubMed
- Romani B, Engelbrecht S, Glashoff RH: Functions of Tat: the versatile protein of human immunodeficiency virus type 1. J Gen Virol. 2010, 91: 1-12. 10.1099/vir.0.016303-0.View ArticlePubMed
- Konkoli Z: Mathematical explanation of the predictive power of the X-level approach reaction noise estimator method. Theor Biol Med Model. 2012, 9: 12-10.1186/1742-4682-9-12.PubMed CentralView ArticlePubMed
- Konkoli Z: Modeling reaction noise with a desired accuracy by using the X level approach reaction noise estimator (XARNES) method. J Theor Biol. 2012, 305: 1-14.View ArticlePubMed
- Konkoli Z, Johannesson H, Lee BP: Fluctuation effects in steric reaction-diffusion systems. Phys Rev E. 1999, 59: R3787-R3790. 10.1103/PhysRevE.59.R3787.View Article
- Wightman F, Ellenberg P, Churchill M, Lewin SR: HDAC inhibitors in HIV. Immunol Cell Biol. 2012, 90: 47-54. 10.1038/icb.2011.95.View ArticlePubMed
- Skupsky R, Burnett JC, Foley JE, Schaffer DV, Arkin AP: HIV promoter integration site primarily modulates transcriptional burst size rather than frequency. Plos Comput Biol. 2010, 6: e1000952-10.1371/journal.pcbi.1000952.PubMed CentralView ArticlePubMed
- Dar RD, Razooky BS, Singh A, Trimeloni TV, McCollum JM, Cox CD, Simpson ML, Weinberger LS: Transcriptional burst frequency and burst size are equally modulated across the human genome. Proc Natl Acad Sci U S A. 2012, 109: 17454-17459. 10.1073/pnas.1213530109.PubMed CentralView ArticlePubMed
- Das AT, Harwig A, Berkhout B: The HIV-1 Tat protein has a versatile role in activating viral transcription. J Virol. 2011, 85: 9506-9516. 10.1128/JVI.00650-11.PubMed CentralView ArticlePubMed
- Peruzzi F: The multiple functions of HIV-1 Tat: proliferation versus apoptosis. Front Biosci. 2006, 11: 708-717. 10.2741/1829.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.