Search

Polyhedral models of metabolic networks are computationally tractable and can predict some cellular functions. A longstanding challenge is incorporating metabolites without losing tractability. In this paper, we do so using a new second-order cone representation of the Michaelis–Menten kinetics. The resulting model consists of linear stoichiometric constraints alongside second-order cone constraints that couple the reaction fluxes to metabolite concentrations. We formulate several new problems around this model: conic flux balance analysis, which augments flux balance analysis with metabolite concentrations; dynamic conic flux balance analysis; and finding minimal cut sets of networks with both reactions and metabolites. Solving these problems yields information about both fluxes and metabolite concentrations. They are second-order cone or mixed-integer second-order cone programs, which, while not as tractable as their linear counterparts, can nonetheless be solved at practical scales using existing software. [ABSTRACT FROM AUTHOR]

We study the dynamics of a model of membrane vesicle transport into dendritic spines, which are bulbous intracellular compartments in neurons driven by molecular motors. We reduce the lubrication model proposed in Fai et al. (Phys Rev Fluids 2:113601, 2017) to a fast–slow system, yielding an analytically and numerically tractable equation equivalent to the original model in the overdamped limit. The model’s key parameters include: (1) the ratio of motors that prefer to push toward the head of the dendritic spine to the motors that prefer to push in the opposite direction, and (2) the viscous drag exerted on the vesicle by the spine constriction. We perform a numerical bifurcation analysis in these parameters and find that steady-state vesicle velocities appear and disappear through several saddle-node bifurcations. This process allows us to identify the region of parameter space in which multiple stable velocities exist. We show by direct calculations that there can only be unidirectional motion for sufficiently close vesicle-to-spine diameter ratios. Our analysis predicts the critical vesicle-to-spine diameter ratio, at which there is a transition from unidirectional to bidirectional motion, consistent with experimental observations of vesicle trajectories in the literature. [ABSTRACT FROM AUTHOR]

We develop a mathematical model to study the immediate effect of low-dose radiation on the G2 checkpoint and the G2/M transition of the cell cycle via a radiation pathway (the ATM–Chk2 pathway) of an individual mammalian cell. The model consists of a system of nonlinear differential equations describing the dynamics of a network of regulatory proteins that play key roles in the G2/M transition, cell cycle oscillations, and the radiation pathway. We simulate the application of a single pulse of low-dose radiation at different intensities (∼ 0–0.4 Gy) and times during the latter part of the G2-phase. We use bifurcation analysis to characterize the effect of radiation on the G2/M transition via the ATM–Chk2 pathway. We show that radiation between 0.1 and 0.3 Gy can delay the G2/M transition, and radiation higher than 0.3 Gy can fully activate the G2 checkpoint. Also, our results show that radiation can be low enough to neither delay the G2/M transition nor activate the G2 checkpoint (∼ 0.1 Gy). Our model supports the idea that the cell response to radiation during G2-phase explains hyper-radiosensitivity and increased radioresistance (HRS/IRR) observed at low dose. [ABSTRACT FROM AUTHOR]