You are asking for the derivative of a non-linear 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 inputs---the 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 Y(a)\delta_a$$ \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 pre-image 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 analysis---one 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 distributions---hence the Dirac functions in the formula. I presume that this is the reason for the reference to distributions in your question).
