I am wondering what is the derivative of the following function with respect to $x(t)$ in sense of distributions. $$ I\left(\int_0^t x(\tau)d\tau \leq c\right) $$ where $I$ is the indicator function and $c$ is a constant.

You are asking for the derivative of a nonlinear function on an infinite dimensional space. (You do not specify the latter but the space locally integrable functions seems a natural candidate). The derivative can only exist in a very weak sense so it is natural to go for the directional derivative which has two inputsthe point $ x $ where the derivative is computed and the direction $ y $ in which the rate of change takes place. A back of the envelope calculation suggests that the derivative should then be $$\sum \frac {Y(a)}{x(a)}\delta_a$$ where we use capitals to indicate the primitives which occur in the formulation and the sum is taken over the $a$ in the preimage of $c$ under $X$. Some general remarks: the question was posed in such a vague manner that it is not really possible to give a precise, rigorous answer. Presumably you have some concrete application in view and I suggest that you check the above formula there to see if it leads to the expected result (ones I looked at were the functions $t$with primitive $t^2$and $\sin t$ for the function $x$). As is evident from the above formula, the derivative is in a much weaker sense than the usual concepts of functional analysisone requires special conditions on $x$ (beyond just smoothness) for the above expression to make sense and the limit of the difference quotients used to define the derivative does not take place in the underlying function space but in a larger space (of distributionshence the Dirac functions in the formula. I presume that this is the reason for the reference to distributions in your question). 


Let me try.
Your mapping is:
Simpler approach, but it might not answer your question:
Suppose $s\mapsto x(\tau,s)$ is a smooth variation of $x(\tau)$, i.e.,
$x\in C^\infty(\mathbb R^2)$.


