Deriving Pulse Responses In Incoherent Feedforward Loops

From Branches to a Transient Output

An incoherent feedforward loop sends an input along two paths with opposing effects on an output: one branch activates it and the other inhibits it. If activation builds quickly while inhibition accumulates more slowly, a sustained input can produce a temporary rise followed by a decline. The key question here is quantitative: which branch gains and response rates create a peak, and when does it occur?

We will use a deliberately simple model to make those dependencies explicit. The two branches respond independently to a step input, and their effects combine linearly at the output. This makes it possible to derive the full time course rather than infer pulse shape from a diagram alone.

With this convention, each branch approaches a steady level proportional to the input:

The rates \(k_a\) and \(k_i\) have units of inverse time; the gains convert input units to branch-state units. We call activation the fast branch and inhibition the slow branch when \(k_a>k_i>0\). A larger rate means a shorter response time scale, approximately \(1/k\), not a larger steady-state effect.

Step Responses and the Steady State

For a constant input, each branch approaches its own target exponentially. Solving the first-order equations gives \(a(t)=g_aU(1-e^{-k_at})\) and \(i(t)=g_iU(1-e^{-k_it})\). Subtracting the inhibitory contribution gives the output trajectory:

At \(t=0\), both branches are zero, so \(y(0)=0\). As time tends to infinity, both exponentials vanish and the output settles to

Thus, equal gains produce exact return to the original baseline in this model, even though the output can rise substantially in between. If the gains differ, the output settles above or below baseline. A transient peak is therefore not automatically the same thing as exact adaptation to baseline.

When Does the Output Have a Peak?

Differentiate the output to find whether it is rising or falling:

Suppose \(U>0\) and \(k_a>k_i>0\). At late times the inhibitory term decays more slowly, so the derivative eventually becomes negative if \(g_i>0\). For the response to start by rising, its initial slope must be positive: \(U(g_ak_a-g_ik_i)>0\). When this condition holds, the derivative crosses zero once. Setting it to zero and rearranging gives the peak time:

The condition makes \(t_*>0\). It also makes the crossing a maximum: after dividing the derivative equation by the positive quantity \(Ue^{-k_it}\), the derivative has the sign of \(g_ak_ae^{-(k_a-k_i)t}-g_ik_i\), which decreases strictly with time. It is positive before \(t_*\) and negative after it. Because \(y(0)=0\) and the output initially rises, this maximum is positive.

Worked Calculations

What the Derivation Tells Us—and What It Does Not

In the stated model, the peak time depends on the gains and rates but not on the magnitude \(U\) of the step: multiplying \(U\) scales the whole trajectory without changing where its derivative is zero. This conclusion relies on linear branch dynamics and a linear output. If branch activation saturates or the output combines the branches nonlinearly, changing the input magnitude can change both peak height and timing.

The model also makes the role of rates distinct from the role of gains. Rates determine how quickly each branch responds; gains determine the eventual branch levels per unit input. The peak criterion compares initial branch drives, \(g_ak_a\) and \(g_ik_i\), whereas the final level depends on the gain difference, \(g_a-g_i\). A large transient can coexist with a zero, positive, or negative final output.

This closed-form calculation is useful for interpreting a mechanism or checking a simulation, but it is not a substitute for checking whether the mechanism is plausible. In a biological pathway, branch dynamics may not be first-order, the branches may share components, and the measured output may not be proportional to a simple difference of states. The derivation identifies what follows from this idealized motif; data or a more detailed model are needed to assess whether its assumptions fit a particular system.

Diagnostics and Common Mistakes

• Confusing a large rate with a large effect. The rate controls how quickly a branch approaches its target; the gain controls the target level. Check units and compare \(k_a\) with \(k_i\) separately from \(g_a\) with \(g_i\).
• Calling every transient maximum exact adaptation. A maximum followed by a decline does not imply return to the starting baseline. Calculate \(y_{\infty}=U(g_a-g_i)\) explicitly.
• Using the peak formula when its conditions fail. A negative logarithm argument, a nonpositive time, or a parameter set that violates \(g_ak_a>g_ik_i\) signals that the assumed positive-time maximum is not present. Do not report the formula as a valid peak time in that case.
• Ignoring initial conditions. The derivation assumes both branches start at zero relative to their pre-input targets. If the system is already displaced from baseline, the exponentials have different coefficients and the peak calculation must be redone for those initial states.
• Trusting a simulation without checking limiting behavior. A numerical trajectory should start at zero under these initial conditions and approach \(U(g_a-g_i)\). A peak, if predicted, should occur near \(t_*\); disagreement can indicate a sign error, incorrect units, or an implementation error.