Step Changes in Transcription and Protein Response

Why a Transcription Step Does Not Produce an Instant Protein Step

A change in transcription changes the supply of mRNA, but protein abundance depends on translation from the mRNA that is present and on protein degradation. If mRNA turns over quickly while protein turns over slowly, the transcript can approach its new level well before the protein does. This separation in time scales is a basic reason that a change in gene expression may have a delayed protein-level response.

We consider a constant transcription rate before time zero, followed by an abrupt switch to a different constant rate at $t=0$. The aim is to derive both trajectories and relate their shapes to the mRNA and protein half-lives.

Model and assumptions: Let $M(t)$ be mRNA abundance and $P(t)$ be protein abundance. After the step, transcription occurs at constant rate $r_1$, mRNA is degraded by first-order turnover rate $k_m>0$, translation occurs at constant rate $k_{tr}>0$ per unit mRNA, and protein is degraded by first-order turnover rate $k_p>0$. Thus, for $t\geq 0$, $$ \frac{dM}{dt}=r_1-k_mM,\qquad \frac{dP}{dt}=k_{tr}M-k_pP. $$ The pre-step system is assumed to have been at steady state, and the step changes transcription only. The rates and translation coefficient remain constant; there are no delays, feedbacks, or additional compartments in this model.

The rate $r_1$ has units of mRNA abundance per time, while $k_m$ and $k_p$ have units of inverse time. The coefficient $k_{tr}$ converts mRNA abundance to protein production per time. For a first-order turnover process, the half-life is $\ln(2)/k$, so the mRNA and protein half-lives are $\ln(2)/k_m$ and $\ln(2)/k_p$, respectively.

mRNA Response to a Step

Let $r_0$ be the transcription rate before the step. Since the pre-step system is at steady state, its initial mRNA abundance is $M_0=r_0/k_m$. The post-step steady-state abundance is $M_1=r_1/k_m$. Solving the mRNA equation with initial condition $M(0)=M_0$ gives:

M(t)=M_1+(M_0-M_1)e^{-k_mt} =M_0+(M_1-M_0)(1-e^{-k_mt}).

The response approaches $M_1$ exponentially. At one mRNA half-life after the step, the difference between $M(t)$ and $M_1$ is half its initial value. The formula applies to an increase or a decrease in transcription: when $M_1<M_0$, the same exponential describes a decline toward the lower steady state.

Deriving the Protein Trajectory

Before the step, steady state requires $k_{tr}M_0=k_pP_0$, so $P_0=k_{tr}M_0/k_p$. The post-step steady state is $P_1=k_{tr}M_1/k_p$. To isolate the response shape, define $F(t)$ as the fraction of the total change in protein abundance completed by time $t$, so that $P(t)=P_0+(P_1-P_0)F(t)$.

The mRNA fraction of its total change is $1-e^{-k_mt}$. Substitution into the protein equation yields a first-order equation for the normalized protein response:

\frac{dF}{dt}=k_p\left[(1-e^{-k_mt})-F(t)\right],\qquad F(0)=0.

Multiplying by the integrating factor $e^{k_pt}$ and integrating from zero to $t$ gives the response explicitly. For unequal turnover rates, the result is:

F(t)=1-\frac{k_pe^{-k_mt}-k_me^{-k_pt}}{k_p-k_m}, \qquad k_p\ne k_m.

This derivation also explains why the response depends on both turnover rates. Protein production is driven by the mRNA response, while protein degradation removes protein throughout the process. The final protein trajectory is therefore not simply the mRNA trajectory multiplied by a constant, except at steady state.

Key result: Under the stated assumptions, $P(t)=P_0+(P_1-P_0)F(t)$, with $F(t)$ given above. The same $F(t)$ describes an upward or downward transcription step. For a step up, $P_1-P_0>0$; for a step down, $P_1-P_0<0$. In either case, $F(t)$ is the fraction of the change from the old protein level toward the new one.

If the turnover rates are equal, the displayed expression for $F(t)$ has a zero denominator, but the response itself is well-defined. Taking the limit as $k_p$ approaches $k_m=k$ gives:

F(t)=1-(1+kt)e^{-kt},\qquad k_p=k_m=k.

This limit is useful as a check on the derivation and prevents a numerical implementation from dividing by a value close to zero when the two rates are nearly equal.

Worked Examples

Worked Example: A Transcription Doubling with Long-Lived Protein

Suppose a hypothetical cell has an mRNA half-life of 1 hour and a protein half-life of 4 hours. Transcription doubles at time zero, and the system was at steady state beforehand. Find the protein abundance after 1 hour as a fraction of its initial abundance.

The turnover rates are $k_m=\ln(2)/1=0.6931\ \text{h}^{-1}$ and $k_p=\ln(2)/4=0.1733\ \text{h}^{-1}$. At $t=1$ hour, $e^{-k_mt}=0.5000$ and $e^{-k_pt}=0.8409$. Substitution into the normalized response gives:

F(1)=1-\frac{(0.1733)(0.5000)-(0.6931)(0.8409)} {0.1733-0.6931} \approx 0.0455.

Because the new steady-state protein abundance is twice the initial abundance, $P_1-P_0=P_0$. Therefore $P(1)/P_0=1+F(1)\approx 1.0455$. The protein level has increased by about 4.55% after 1 hour, even though the eventual increase is 100%. The small early response is consistent with a transcript that is changing and a protein with a longer half-life.

Worked Example: Transcription Shutoff with Slowly Turned-Over mRNA

In a hypothetical system, transcription is completely shut off at time zero. The mRNA half-life is 2 hours and the protein half-life is 8 hours. Find the remaining protein abundance after 8 hours as a fraction of its initial abundance.

Here $k_m=\ln(2)/2=0.3466\ \text{h}^{-1}$ and $k_p=\ln(2)/8=0.08664\ \text{h}^{-1}$. Shutoff means $M_1=P_1=0$, given the linear model and the pre-step steady state. At $t=8$ hours, $e^{-k_mt}=0.06250$ and $e^{-k_pt}=0.5000$. Thus:

F(8)=1-\frac{(0.08664)(0.06250)-(0.3466)(0.5000)} {0.08664-0.3466} \approx 0.3541.

The protein has completed about 35.41% of its total decrease toward zero. Its remaining fraction is $P(8)/P_0=1-F(8)\approx 0.6459$, or about 64.59%. Although mRNA has fallen to 6.25% of its initial abundance after four mRNA half-lives, protein remains because degradation is slower and protein continues to be produced from the remaining mRNA.

Worked Example: A Threefold Transcription Increase

Consider a hypothetical system with mRNA and protein half-lives of 3 hours and 6 hours, respectively. Transcription increases threefold at time zero. Calculate the protein abundance after 3 hours relative to its initial abundance.

The rates are $k_m=\ln(2)/3=0.2310\ \text{h}^{-1}$ and $k_p=\ln(2)/6=0.1155\ \text{h}^{-1}$. At 3 hours, $e^{-k_mt}=0.5000$ and $e^{-k_pt}=0.7071$. The response fraction is:

F(3)=1-\frac{(0.1155)(0.5000)-(0.2310)(0.7071)} {0.1155-0.2310} \approx 0.0858.

Since the new steady-state protein abundance is three times the initial abundance, $P_1-P_0=2P_0$. Hence $P(3)/P_0=1+2F(3)\approx 1.172$. Protein has increased by about 17.2% after 3 hours, not by threefold. The distinction between the eventual level and the transient response is essential when interpreting a measurement taken before the system reaches its new steady state.

Shape, Checks, and Practical Limitations

For positive turnover rates, the mRNA and protein responses move monotonically toward their new steady states in this model. The protein response starts with $F(0)=0$; because the initial protein is at steady state while mRNA has not yet changed, its initial rate of change is zero. This produces a lag. The new steady state is approached as both exponential terms decay. For positive rates and a fixed transcription step, the model does not predict overshoot.

A useful diagnostic is to compare the observed time course with these qualitative features. An apparent immediate protein jump, substantial overshoot, or persistent oscillation cannot be explained by this two-state, first-order step model alone. Such patterns may reflect measurement timing, an unmodelled mechanism, or a violation of the assumptions; they should not be “fixed” by misreading $F(t)$ as an instantaneous protein response.

The analytical formula is convenient when the rates and step timing are specified. Numerical integration of the two ODEs is a general alternative and is especially useful when inputs vary continuously or the model has additional states. A quasi-steady approximation for protein would treat protein as tracking its current production closely; that approximation is not appropriate when protein turnover is slow relative to the changing mRNA. Neither approach removes the need to check the model assumptions.

Common pitfall: After transcription doubles, it is incorrect to conclude that protein instantly doubles, or that protein follows the mRNA exponential with the same half-life. The twofold value is the new steady state, not the immediate response. Calculate the mRNA trajectory first, then use the coupled protein response; the protein turnover rate contributes its own lag.

The analytic solution relies on a pre-step steady state, an abrupt and sustained transcription change, first-order degradation, and constant translation and turnover parameters. If the initial state is not at steady state, the same ODEs can still be solved, but the initial conditions must be retained rather than replaced by $M_0=r_0/k_m$ and $P_0=k_{tr}M_0/k_p$. Likewise, a time-varying transcription input requires a time-varying input solution rather than this single-step formula. These limitations matter when interpreting time courses: a mismatch may indicate that the step-response assumptions do not describe the experiment.

Check Your Understanding

Answer each question, then compare with the solutions below.

If mRNA and protein have equal turnover rate $k$, what expression gives $F(t)$, and why is it needed?
A transcription step changes the eventual protein abundance from $P_0$ to $0.5P_0$. If $F(t)=0.4$, what is $P(t)$ relative to $P_0$?
Why can mRNA fall substantially while protein remains relatively abundant after transcription shutoff?
Name two assumptions whose failure could make the analytic step-response formula inappropriate.

Solutions. (1) $F(t)=1-(1+kt)e^{-kt}$. It is the limit of the unequal-rate expression when $k_p$ approaches $k_m$, so it avoids the apparent zero denominator. (2) $P(t)=P_0+(0.5P_0-P_0)(0.4)=0.8P_0$; 20% of the total decrease has occurred. (3) Protein degradation may be slower than mRNA degradation, and the residual mRNA can also continue to produce protein. (4) Examples include a non-steady initial state, a gradual rather than abrupt change in transcription, non-first-order degradation, or translation and turnover rates that change during the response.