Despite high-resolution structural studies on the USP14-proteasome-substrate complexes, time-resolved cryo-electron microscopy (cryo-EM) results on USP14-regulated allostery of the 26S proteasome are still very limited and a quantitative understanding of substrate degradation dynamics remains elusive. In this study, we propose a mean field model of ordinary differential equations (ODEs) for USP14 regulated substrate degradation in 26S proteasome. The kinetic model incorporates recent cryo-EM findings on the allostery of 26S proteasome and generates results in good agreement with time-resolved experimental observations. The model elucidates that USP14 typically reduces the substrate degradation rate and reveals the functional dependence of this rate on the concentrations of substrate and adenosine triphosphate (ATP). The half-maximal effective concentration (EC50) of the substrate for different ATP concentrations is predicted. When multiple substrates are present, the model suggests that substrates that are easier to insert into the OB-ring and disengage from the proteasome, or less likely to undergo deubiquitination would be more favored to be degraded by the USP14-bound proteasome. The mean field model proposed here quantitatively considers the process of proteasomal substrate degradation from the perspective of chemical kinetics, and provides a quantitative framework to decode the dynamic interplay between USP14 and the proteasome.

Once challenged by the SARS-CoV-2 virus, the human host immune system triggers a dynamic process against infection. We constructed a mathematical model to describe host innate and adaptive immune response to viral challenge. Based on the dynamic properties of viral load and immune response, we classified the resulting dynamics into four modes, reflecting increasing severity of COVID-19 disease. We found the numerical product of immune system's ability to clear the virus and to kill the infected cells, namely immune efficacy, to be predictive of disease severity. We also investigated vaccine-induced protection against SARS-CoV-2 infection. Results suggested that immune efficacy based on memory T cells and neutralizing antibody titers could be used to predict population vaccine protection rates. Finally, we analyzed infection dynamics of SARS-CoV-2 variants within the construct of our mathematical model. Overall, our results provide a systematic framework for understanding the dynamics of host response upon challenge by SARS-CoV-2 infection, and this framework can be used to predict vaccine protection and perform clinical diagnosis.

Motivation: One of the long-expected goals of genome-scale metabolic modelling is to evaluate the influence of the perturbed enzymes on flux distribution. Both ordinary differential equation (ODE) models and constraint-based models, like Flux balance analysis (FBA), lack the capacity to perform metabolic control analysis (MCA) for large-scale networks. Results: In this study, we developed a hyper-cube shrink algorithm (HCSA) to incorporate the enzymatic properties into the FBA model by introducing a pseudo reaction V constrained by enzymatic parameters. Our algorithm uses the enzymatic information quantitatively rather than qualitatively. We first demonstrate the concept by applying HCSA to a simple three-node network, whereby we obtained a good correlation between flux and enzyme abundance. We then validate its prediction by comparison with ODE and with a synthetic network producing voilacein and analogues in Saccharomyces cerevisiae. We show that HCSA can mimic the state-state results of ODE. Finally, we show its capability of predicting the flux distribution in genome-scale networks by applying it to sporulation in yeast. We show the ability of HCSA to operate without biomass flux and perform MCA to determine rate-limiting reactions. Availability and implementation: Algorithm was implemented by Matlab and C ++. The code is available at https://github.com/kekegg/HCSA. Contact: [email protected] or [email protected]. Supplementary information: Supplementary data are available at Bioinformatics online.

BACKGROUND: Single-particle cryo-electron microscopy (cryo-EM) has become a mainstream tool for the structural determination of biological macromolecular complexes. However, high-resolution cryo-EM reconstruction often requires hundreds of thousands of single-particle images. Particle extraction from experimental micrographs thus can be laborious and presents a major practical bottleneck in cryo-EM structural determination. Existing computational methods for particle picking often use low-resolution templates for particle matching, making them susceptible to reference-dependent bias. It is critical to develop a highly efficient template-free method for the automatic recognition of particle images from cryo-EM micrographs. RESULTS: We developed a deep learning-based algorithmic framework, DeepEM, for single-particle recognition from noisy cryo-EM micrographs, enabling automated particle picking, selection and verification in an integrated fashion. The kernel of DeepEM is built upon a convolutional neural network (CNN) composed of eight layers, which can be recursively trained to be highly "knowledgeable". Our approach exhibits an improved performance and accuracy when tested on the standard KLH dataset. Application of DeepEM to several challenging experimental cryo-EM datasets demonstrated its ability to avoid the selection of un-wanted particles and non-particles even when true particles contain fewer features. CONCLUSIONS: The DeepEM methodology, derived from a deep CNN, allows automated particle extraction from raw cryo-EM micrographs in the absence of a template. It demonstrates an improved performance, objectivity and accuracy. Application of this novel method is expected to free the labor involved in single-particle verification, significantly improving the efficiency of cryo-EM data processing.

Massive technological advances enabled high-throughput measurements of proteomic changes in biological processes. However, retrieving biological insights from large-scale protein dynamics data remains a challenging task. Here we used the mating differentiation in yeast Saccharomyces cerevisiae as a model and developed integrated experimental and computational approaches to analyze the proteomic dynamics during the process of cell fate determination. When exposed to a high dose of mating pheromone, the yeast cell undergoes growth arrest and forms a shmoo-like morphology; however, at intermediate doses, chemotropic elongated growth is initialized. To understand the gene regulatory networks that control this differentiation switch, we employed a high-throughput microfluidic imaging system that allows real-time and simultaneous measurements of cell growth and protein expression. Using kinetic modeling of protein dynamics, we classified the stimulus-dependent changes in protein abundance into two sources: global changes due to physiological alterations and gene-specific changes. A quantitative framework was proposed to decouple gene-specific regulatory modes from the growth-dependent global modulation of protein abundance. Based on the temporal patterns of gene-specific regulation, we established the network architectures underlying distinct cell fates using a reverse engineering method and uncovered the dose-dependent rewiring of gene regulatory network during mating differentiation. Furthermore, our results suggested a potential crosstalk between the pheromone response pathway and the target of rapamycin (TOR)-regulated ribosomal biogenesis pathway, which might underlie a cell differentiation switch in yeast mating response. In summary, our modeling approach addresses the distinct impacts of the global and gene-specific regulation on the control of protein dynamics and provides new insights into the mechanisms of cell fate determination. We anticipate that our integrated experimental and modeling strategies could be widely applicable to other biological systems.

Cell polarization toward an attractant is influenced by both physical and chemical factors. Most existing mathematical models are based on reaction-diffusion systems and only focus on the chemical process occurring during cell polarization. However, membrane tension has been shown to act as a long-range inhibitor of cell polarization. Here, we present a cell polarization model incorporating the interplay between Rac GTPase, filamentous actin (F-actin), and cell membrane tension. We further test the predictions of this model by performing single cell measurements of the spontaneous polarization of cancer stem cells (CSCs) and non-stem cancer cells (NSCCs), as the former have lower cell membrane tension. Based on both our model and the experimental results, cell polarization is more sensitive to stimuli under low membrane tension, and high membrane tension improves the robustness and stability of cell polarization such that polarization persists under random perturbations. Furthermore, our simulations are the first to recapitulate the experimental results described by Houk et al., revealing that aspiration (elevation of tension) and release (reduction of tension) result in a decrease in and recovery of the activity of Rac-GTP, respectively, and that the relaxation of tension induces new polarity of the cell body when a cell with the pseudopod-neck-body morphology is severed.

As significant growth in mobile data services and demands emerges, mobile Quality-of-Service (QoS) as well as mobile Quality-of-Experience (QoE) receive rapidly growing attentions from both academia and industry. This paper throws light upon various relations between mobile web QoS and QoE according to different kinds of web services, from a comprehensive point of view. We proposes an integrated mobile web QoE assessment method in a fine-grained manner. Works presented in previous literatures are first discussed, and then it is demonstrated that a dual-process theory in QoE assessment which separates the whole process into two rather independent ones, aligning closely with reality. Next, relations between system capability and user perception are studied through mathematical models, and these relations are demonstrated with functions of MOS and waiting time acquired in the dual-process model. The theoretical models proposed in this paper are then proved applicable for mobile web QoE assessment by comparisons conducted between results of our experiments and those of classic ones.

One of the major breakthroughs in oncogenesis research in recent years is the discovery that, in most patients, oncogenic mutations are concentrated in a few core biological functional pathways. This discovery indicates that oncogenic mechanisms are highly related to the dynamics of biologic regulatory networks, which govern the behaviour of functional pathways. Here, we propose that oncogenic mutations found in different biological functional pathways are closely related to parameter sensitivity of the corresponding networks. To test this hypothesis, we focus on the DNA damage-induced apoptotic pathway--the most important safeguard against oncogenesis. We first built the regulatory network that governs the apoptosis pathway, and then translated the network into dynamics equations. Using sensitivity analysis of the network parameters and comparing the results with cancer gene mutation spectra, we found that parameters that significantly affect the bifurcation point correspond to high-frequency oncogenic mutations. This result shows that the position of the bifurcation point is a better measure of the functionality of a biological network than gene expression levels of certain key proteins. It further demonstrates the suitability of applying systems-level analysis to biological networks as opposed to studying genes or proteins in isolation.

Escherichia coli chemotactic motion in spatiotemporally varying environments is studied by using a computational model based on a coarse-grained description of the intracellular signaling pathway dynamics. We find that the cell's chemotaxis drift velocity v(d) is a constant in an exponential attractant concentration gradient [L] proportional, variantexp(Gx). v(d) depends linearly on the exponential gradient G before it saturates when G is larger than a critical value G(C). We find that G(C) is determined by the intracellular adaptation rate k(R) with a simple scaling law: G(C) infinity k(1/2)(R). The linear dependence of v(d) on G = d(ln[L])/dx directly demonstrates E. coli's ability in sensing the derivative of the logarithmic attractant concentration. The existence of the limiting gradient G(C) and its scaling with k(R) are explained by the underlying intracellular adaptation dynamics and the flagellar motor response characteristics. For individual cells, we find that the overall average run length in an exponential gradient is longer than that in a homogeneous environment, which is caused by the constant kinase activity shift (decrease). The forward runs (up the gradient) are longer than the backward runs, as expected; and depending on the exact gradient, the (shorter) backward runs can be comparable to runs in a spatially homogeneous environment, consistent with previous experiments. In (spatial) ligand gradients that also vary in time, the chemotaxis motion is damped as the frequency omega of the time-varying spatial gradient becomes faster than a critical value omega(c), which is controlled by the cell's chemotaxis adaptation rate k(R). Finally, our model, with no adjustable parameters, agrees quantitatively with the classical capillary assay experiments where the attractant concentration changes both in space and time. Our model can thus be used to study E. coli chemotaxis behavior in arbitrary spatiotemporally varying environments. Further experiments are suggested to test some of the model predictions.

Attribution of biological robustness to the specific structural properties of a regulatory network is an important yet unsolved problem in systems biology. It is widely believed that the topological characteristics of a biological control network largely determine its dynamic behavior, yet the actual mechanism is still poorly understood. Here, we define a novel structural feature of biological networks, termed 'regulation entropy', to quantitatively assess the influence of network topology on the robustness of the systems. Using the cell-cycle control networks of the budding yeast (Saccharomyces cerevisiae) and the fission yeast (Schizosaccharomyces pombe) as examples, we first demonstrate the correlation of this quantity with the dynamic stability of biological control networks, and then we establish a significant association between this quantity and the structural stability of the networks. And we further substantiate the generality of this approach with a broad spectrum of biological and random networks. We conclude that the regulation entropy is an effective order parameter in evaluating the robustness of biological control networks. Our work suggests a novel connection between the topological feature and the dynamic property of biological regulatory networks.

Drug molecules not only interact with specific targets, but also alter the state and function of the associated biological network. How to design drugs and evaluate their functions at the systems level becomes a key issue in highly efficient and low-side-effect drug design. The arachidonic acid metabolic network is the network that produces inflammatory mediators, in which several enzymes, including cyclooxygenase-2 (COX-2), have been used as targets for anti-inflammatory drugs. However, neither the century-old nonsteriodal anti-inflammatory drugs nor the recently revocatory Vioxx have provided completely successful anti-inflammatory treatment. To gain more insights into the anti-inflammatory drug design, the authors have studied the dynamic properties of arachidonic acid (AA) metabolic network in human polymorphous leukocytes. Metabolic flux, exogenous AA effects, and drug efficacy have been analyzed using ordinary differential equations. The flux balance in the AA network was found to be important for efficient and safe drug design. When only the 5-lipoxygenase (5-LOX) inhibitor was used, the flux of the COX-2 pathway was increased significantly, showing that a single functional inhibitor cannot effectively control the production of inflammatory mediators. When both COX-2 and 5-LOX were blocked, the production of inflammatory mediators could be completely shut off. The authors have also investigated the differences between a dual-functional COX-2 and 5-LOX inhibitor and a mixture of these two types of inhibitors. Their work provides an example for the integration of systems biology and drug discovery.