Beta Phase: Square45 is currently in beta testing. Expect some features or content to be incomplete or missing.
45

Cell Signaling

Field: Cell Biology

The ability of a cell to receive, process, and transmit signals with its environment and with itself.

Sequence of Expressions

Let [S][S] be the concentration of the second messenger (e.g., Ca2+Ca^{2+} or cAMPcAMP). Its dynamics are governed by the balance between production (catalyzed by enzymes EE) and degradation: \n\nd[S]dt=(Vmax[L]KM+[L])−kdeg[S]\frac{d[S]}{dt} = \left( \frac{V_{max} [L]}{K_M + [L]} \right) - k_{deg} [S] \n\nwhere VmaxV_{max} is the maximum production rate (e.g., adenylyl cyclase activity), KMK_M is the Michaelis constant for the activating ligand [L][L], and kdegk_{deg} is the first-order degradation rate constant. The system exhibits signal amplification when the production term is non-linear and significantly exceeds the degradation term over a short time scale.
Consider the binding equilibrium between a receptor RR, a ligand LL, and a drug DD. The fractional occupancy θ\theta of the receptor sites is defined by the generalized binding equation: \n\nθ=[L]+KA[D]KD+[L]+KA[D]\theta = \frac{[L] + K_{A} [D]}{K_{D} + [L] + K_{A} [D]} \n\nwhere [L][L] is the concentration of the natural ligand, [D][D] is the concentration of the drug, KDK_{D} is the dissociation constant for the natural ligand, and KAK_{A} is the apparent affinity constant for the drug, derived from the ratio of the binding rate constants KA=kon,D/koff,DK_{A} = k_{on, D} / k_{off, D}. For an antagonist, the binding term is modeled such that the drug binding does not induce a conformational change leading to signal activation.
The relationship between the measured response RR (e.g., enzyme activity or binding) and the ligand concentration [L][L] is described by the Hill equation: \n\nR([L])=[L]λK0.5λ+[L]λR([L]) = \frac{[L]^{\lambda}}{K_{0.5}^{\lambda} + [L]^{\lambda}} \n\nwhere K0.5K_{0.5} is the concentration of ligand required to achieve half-maximal response (R/2R/2), and λ\lambda is the Hill coefficient. The coefficient λ\lambda quantifies the degree of cooperativity: λ=1\lambda = 1 indicates non-cooperative binding (hyperbolic), λ>1\lambda > 1 indicates positive cooperativity (sigmoidal), and λ<1\lambda < 1 indicates negative cooperativity.
Define the binding interaction between a receptor RR and a ligand LL as an equilibrium process: R+L⇌RL R + L \rightleftharpoons RL The binding affinity is characterized by the dissociation constant KD=[R][L][RL]K_D = \frac{[R][L]}{[RL]}. The geometric complementarity constraint dictates that the potential energy of interaction, U(r)U(r), must exhibit a minimum at the equilibrium distance r0r_0, and this minimum must be highly specific, such that the binding free energy ΔGbind=−RTln(1/KD)\Delta G_{bind} = -RT \text{ln}(1/K_D) is maximized only when the spatial coordinates rR\mathbf{r}_R and rL\mathbf{r}_L satisfy a geometric matching function G(rR,rL)=0\mathcal{G}(\mathbf{r}_R, \mathbf{r}_L) = 0.
Let [Mi][M_i] be the concentration of molecular species MiM_i at time tt. A cascade of NN sequential reactions is modeled by the system of coupled ODEs: d[M1]dt=k1[S]−k2[M1]\frac{d[M_1]}{dt} = k_1 [S] - k_2 [M_1] d[Mi]dt=ki−1[Mi−1]−ki[Mi], for i=2, to N\frac{d[M_i]}{dt} = k_{i-1} [M_{i-1}] - k_i [M_i], \text{ for } i=2, \text{ to } N d[MN+1]dt=kN[MN]−kN+1[MN+1]\frac{d[M_{N+1}]}{dt} = k_N [M_N] - k_{N+1} [M_{N+1}] where [S][S] is the initial signal concentration, and kik_i are the rate constants for phosphorylation/activation steps. The overall response is determined by the steady-state concentration of the final effector, [MN+1]ss[M_{N+1}]_{ss}.
The relationship between the input signal concentration [S][S] and the resulting cellular activity AA is modeled by the Hill equation, which incorporates a cooperative binding mechanism: A([S])=Amax[S]hKDh+[S]h A([S]) = A_{max} \frac{[S]^h}{K_D^h + [S]^h} where AmaxA_{max} is the maximum achievable activity, KDK_D is the signal concentration required for half-maximal response, and hh is the Hill coefficient. A threshold effect is mathematically manifested when h>1h > 1, indicating a cooperative, non-linear switch-like response.
Let C(t)=([M1],[M2],…,[MN])\mathbf{C}(t) = \left( [M_1], [M_2], \dots, [M_N] \right) be the vector of concentrations of NN molecular species at time tt. The dynamics of the pathway are governed by the system of Ordinary Differential Equations (ODEs): \n\nd[Mi]dt=∑j(kon,ij[Mj][L]−koff,ij[Mi])+∑k(vk−kdeg,k[Mi])\frac{d[M_i]}{dt} = \sum_{j} \left( k_{on, ij} [M_j] [L] - k_{off, ij} [M_i] \right) + \sum_{k} \left( v_{k} - k_{deg, k} [M_i] \right) \n\nwhere kon,ijk_{on, ij} and koff,ijk_{off, ij} are the association and dissociation rate constants for the interaction between MjM_j and LL (ligand/precursor), and vkv_k represents the production rate from a preceding step kk. The overall pathway response RR is a function of the final activated species: R=f(C(t))R = f(\mathbf{C}(t)).
The output signal RR resulting from an input stimulus II is modeled by a threshold function Θ\Theta: \n\nR(I)={Rmaxif I≥Ithreshold0if I<IthresholdR(I) = \begin{cases} R_{max} & \text{if } I \ge I_{threshold} \\ 0 & \text{if } I < I_{threshold} \end{cases} \n\nAlternatively, using a continuous approximation for mathematical analysis, the response can be modeled by a steep sigmoidal function, such as the Fermi function, centered at the threshold IthresholdI_{threshold}: \n\nR(I)=Rmax1+e−k(I−Ithreshold)R(I) = \frac{R_{max}}{1 + e^{-k(I - I_{threshold})}} \n\nwhere kk controls the steepness of the transition.
Let S(t)S(t) be the concentration of the initial signal molecule and R(t)R(t) be the concentration of the final cellular response product. If the reaction cascade follows a rate law r=k[S]nr = k[S]^n, where n>1n > 1 is the effective reaction order, the relationship between the initial input S0S_0 and the final output RfinalR_{final} can be modeled by the ratio of rates: d[R]dt=keff[A][S]nKMn+[S]n \frac{d[R]}{dt} = k_{eff} [A] \frac{[S]^n}{K_M^n + [S]^n} where keffk_{eff} is the effective rate constant and nn quantifies the amplification exponent. For significant amplification, we require n>1n > 1 and keff[S]nKMn+[S]n→large constantk_{eff} \frac{[S]^n}{K_M^n + [S]^n} \to \text{large constant} as [S][S] increases.
Consider a system of two molecular species, AA and BB, governed by the following coupled ODE system: d[A]dt=k1[S]−k2[A][B]\frac{d[A]}{dt} = k_1 [S] - k_2 [A][B] d[B]dt=k3[A][B]−k4[B]\frac{d[B]}{dt} = k_3 [A][B] - k_4 [B] For negative feedback, the rate constant k3k_3 must be dependent on the concentration of the product BB, such that k3=k3,max[B]KB+[B]k_3 = k_{3,max} \frac{[B]}{K_B + [B]}, leading to a self-limiting production rate for BB that stabilizes the system around a steady state [B]ss[B]_{ss}.