Complex problems to challenge AI models
Protein aggregation problem
We are investigating aggregation of Amyloid Beta ($A\beta$) proteins in neuronal cells under fluctuating physiological and environmental conditions. The rate of change of $A\beta$ production is affected by the levels of misfolded $A\beta$ proteins $M(t)$ in the cell and reactive oxygen species (ROS) $R_t$:
\[
\text{Rate of change A}\beta\,\text{production }=\frac{{57.9}}{1+e^{0.6\cdot\,M(t)}+0.8\cdot\,R(t)}
\]
Rate of change of protein degradation increases with oxidative stress and is defined as:
\[
\text{Rate of change in A}\beta \text{ depletion}=\frac{1}{1+e^{0.4\cdot(4-R(t))}}\cdot\,A(t)
\]
The rate of change of misfolded proteins increases with elevated ROS, disturbed pH and temperature stress:
\[
\frac{dM}{dt}=\frac{A(t)}{1+e^{4-R(t)}}\cdot\,{1-e^{-(pH(t)-7.4)^2}}\cdot\,{1-e^{-0.05\cdot\,(T(t)-37)^2}}
\]
The rate of change of misfolded protein is infu
ROS level $R(t)$ fluctuates with metabolic cycles defined as:
\[
R(t)=0.3\cdot\,\left[3.6-0.7\cdot\,\sin\left({\frac{2\pi\,t}{9}}\right)\right]
\]
Change in the Chaperone protein $H(t)$ levels is influenced by the oxidative stress and misfolded protein in the cell by following relation:
\[
\frac{dH(t)}{dt}=\frac{0.08\cdot\, (R(t)+0.07\cdot\,M(t))}{1+e^{2\cdot\,(pH(t)-7)}}-0.2\cdot\,H(t)
\]
The turnover rate $k_{cat}$ for activated chaperone is found to be 0.56 $min^{-1}$. Chaperone-mediated clearance rate ($V_c$) is defined
\[
V_c=\frac{V_{max}\cdot\,M(t)}{0.3+M(t)}
\]
where
$V_{max}=k_{cat}\cdot\,H(t)$
Due to incubator malfunction, temperature fluctuates sinusoidally over time:
\[
T(t) = T_{\text{avg}} + A \cdot \cos\left( \frac{2\pi}{60}(t - 20) \right)
\]
$T(t)$ : Temperature at time $t$
$T_{\text{avg}}$ : Average temperature of the day
$A$ : Amplitude (half of the daily temperature range)
$t$ : Time in minutes
$t_{\text{peak}}$ : Time of the peak temperature
We noted that the average temperature of the incubator is $35^\circ\text{C}$ and maximum temperature was noted to be $43^\circ\text{C}$.
The misfolded proteins in the cell have tendency to stick together and form clusters following Smoluchowski coagulation equation:
For each cluster of size $k>1$
\[
\frac{dC_k(t)}{dt}=\frac{1}{2}\sum_{i+j=k} K(i,j)\cdot,c_i(t)\cdot\,c_j(t)-c_k(t)\sum_{j} K(k,j)\cdot,c_j(t)
\]
Where,
$C_k$ is cluster of size $k$
$K(i,j)$: aggregation kernel between clusters of size i and j
The aggregation kernel is a mathematical function that describes the rate at which two clusters of sizes stick together and is modeled by following relationship:
\[
K(i,j,t)=(i^{1/3}+j^{1/3})\cdot\,\frac{T(t)}{T_0}\cdot\,e^{-0.35\cdot(pH(t)-7.4)^2}, T_0 = 37
\]
As the cluster size and number grow, they contribute to the acidification of the cellular environment.
\[
pH(t)=7+ 0.5\cdot\sin\left({\frac{2\pi\, t}{60}}\right)-10\cdot\,\sum_{k=5}^{N}c_k(t)
\]
Initial conditions at \( t = 0 \):
\[
\begin{aligned}
A(0) &= 5 \\
M(0) &= 1 \\
H(0) &= 0.05 \\
C_k(0) &= 0
\end{aligned}
\]
What is the concentration of $C_3$, the cluster of size 3, present in the cell after 2 hour?
Quorum sensing problem
We are studying a mixed population of two bacterial strains M and N. When cultured together, tthese strains compete for limited resources following the Lotka-Volterra competition model:
$$\frac{dN_m}{dt}=r_m\cdot N_m\left(1-\frac{N_m+\gamma N_n}{K}\right), \frac{dN_n}{dt}=r_n\cdot N_n\left(1-\frac{N_n+\delta N_m}{K}\right) $$
with parameters:
$r_m =0.15, \gamma=0.3, r_n=0.5, \delta=0.8, K=1\times10^6$
Each bacterial strain expresses surface transporters that mediate glucose uptake via receptor-ligand binding:
$$\frac{dRL_m}{dt}=N_m\cdot 0.6\times10^{-5}\cdot\,L(t)-0.3\cdot RL_m$$, $$\frac{dRL_n}{dt}=N_n\cdot 0.9\times10^{-5}\cdot\,L(t)-0.4\cdot RL_n$$
where,
$RL_m$ and $RL_n$ are the receptor-ligand complexes of strain M and N respectively, and
$L(t)=(5+ 2\cdot\ sin(\frac{t}{10}))\cdot 10^6$ is the fluctuating glucose concentration over time (in minutes).
Glucose-bound receptors trigger production of protein B via first-order kinetics::
$$\frac{dB_m}{dt}=0.025\cdot RL_m-0.3\cdot B_m, \frac{dB_n}{dt}=0.03\cdot RL_n-0.4\cdot B_n$$
Strains M and N synthesize autoinducers $A_m$ and $A_n$, regulated by their respective B protein levels through a Hill function with Hill coefficient n=2:
$$\frac{dA_m}{dt}=2.02\times10^6\cdot\,\frac{B_m^n}{(2\times10^6)^n+B_m^n}-0.4\cdot A_m$$
$$\frac{dA_n}{dt}=2.6\times10^6\cdot\,\frac{B_n^n}{(2.3\times10^6)^n+B_n^n}-0.2\cdot A_n$$
Once $A_m$ and $A_n$ cross thresholds of $4 \times 10^6$ and $2.5 \times 10^6$ respectively, a fraction of the bacterial populations begins secreting virulence factors:
$$\text{Strain M: } 0.6\% \text{ of cells secrete } V_m$$
$$\text{Strain N: } 1.6\% \text{ of cells secrete } V_n$$
The virulence factors degrade at the following rates:
$$\frac{dV_m}{dt} = -0.2 \cdot V_m, \quad \frac{dV_n}{dt} = -0.3 \cdot V_n$$
The virulence factors collectively damage host HEK cells:
$$\frac{dC}{dt}=1\times10^{-7}\cdot\ V_t\cdot C$$
where C is the number of HEK cells and $V_t$ is $V_m+V_n$
Above a strain-specific toxicity threshold $\theta$, the virulence becomes self-toxic to the bacteria:
$$\frac{dN_m}{dt}=-\lambda_m\cdot\ (V_t-\theta_m)\cdot N_m, \frac{dN_n}{dt}=-\lambda_n\cdot\ (V_t-\theta_n)\cdot N_n$$
where,
$\lambda_m=1\times10^{-5}, \lambda_n=1.5\times10^{-5}, \theta_m=3\times10^4, \theta_n=3.2\times10^4$
We are investigating the efficacy of drug $I$ on these bacterial strains which reduces the effective concentration of $B_m$ and $B_n$ by the following equation:
$$B_m(effective)=\frac{B_m}{1+\frac{I}{K_m}}, B_n(effective)=\frac{B_n}{1+\frac{I}{K_n}}$$
We have maintained HEK cell culture seeded at the initial concentration of $1\times10^7$ cells. At $t=0$, we added 10 CFU of bacterial strain M. We let the bacterial strain M grow for an hour in the media. After that, we added 1000 CFU of strain N with and let both of them grow together for an hour. After that, we introduce the drug that follows sinusoidal release profile:
$$I(t)=1.5+1.2\cdot sin(\frac{t}{10})$$
We observed that CFU for HEK cells was measured at $7.52\times10^6$ after 1 hour of drug addition. Estimate the ratio of $\frac{K_m}{K_n}$
Protein synergy problem
You are tasked with estimating the degradation of fibrinogen in human plasma initiated by a snake venom serine protease (SVSP), which also activates prothrombin to thrombin. The system involves three key reactions: (R1) SVSP activates prothrombin to thrombin, (R2) thrombin degrades fibrinogen, and (R3) SVSP directly degrades fibrinogen. Initially, 0.5 mM SVSP and 1.7 g/mL fibrinogen are present at pH 7.4 and 37°C. Prothrombin is present at 1.4mM in the plasma.
The reaction environment changes dynamically: (1) pH decreases linearly by 0.01 units per minute due to accumulating acidic byproducts (2) Temperature remains at 37°C until 60 minutes, then increases at 0.05°C per minute. At 45 minutes, the serine protease inhibitor PMSF is added, inhibiting SVSP with a first-order rate constant of 0.03 min⁻¹.
\[
I_1(t) = 1 - e^{-0.03 \cdot (t - 45)} \quad \text{for } t \geq 45\text{ min}
\]
At 90 minutes, the thrombin-specific inhibitor Hirudin is introduced:
\[
I_2(t) = 1 - e^{-0.2 \cdot (t - 90)} \quad \text{for } t \geq 90\text{ min}
\]
Where,
$I_1$ and $I_2$ are the fraction of SVSP and Thrombin inhibited respectively.
SVSP follows Michaelis-Menten kinetics with substrate inhibition:
\[
\text{Rate}_{R3} = \frac{0.1 \cdot [E] \cdot [F]}{1.3 + [F] + \frac{[F]^2}{5}}
\]
Prothrombin Activation by SVSP (R1):
\[
\frac{d[P]}{dt} = -0.5 \cdot [E] \cdot [P] \cdot f(pH, T)
\]
$[P]$ : Prothrombin concentration (mM)
$[E]$ : SVSP concentration (mM)
$k_{\text{act}}$ : Rate constant for prothrombin activation by SVSP (mM$^{-1}$ min$^{-1}$)
$f(pH, T)$ : Environmental modulation factor accounting for pH and temperature effects on enzyme activity
Thrombin follows standard Michaelis-Menten kinetics:
\[
\text{Rate}_{R2} = \frac{0.1 \cdot [T] \cdot [F]}{1.0 + [F]}
\]
Where,
$[F]$ is fibrinogen amount
$[T]$ is thrombin concentration
$[E]$ is SVSP concentration
.
SVSP degradation follows the Arrhenius equation:
\[
k_{\text{deg}}(T) = 7.4\cdot e^{-\frac{E_a}{R \cdot (T + 273.15)}}
\]
where
$k_{deg}$ is the rate constant for SVSP degradation, $E_a$ is 17 KJ/mol, $R$ is 8.314 J/mol/K and T is in degree celsius.
SVSP activities are modulated by pH and temperature through the environmental factor:
\[
f(pH, T) = e^{-0.5 \cdot (pH - 7.4)^2} \cdot e^{\frac{-5000}{8.314 \cdot (T + 273.15)}}
\]
Estimate the concentration of fibrinogen remaining after 2 hours.