Cellular signaling, crucial for biological processes like immune response and homeostasis, relies on specificity and fidelity in signal transduction to accurately respond to stimuli amidst biological ...noise. Kinetic proofreading (KPR) is a key mechanism enhancing signaling specificity through time-delayed steps, although its effectiveness is debated due to intrinsic noise potentially reducing signal fidelity. In this study, we reformulate the theory of kinetic proofreading (KPR) by convolving multiple intermediate states into a single state and then define an overall “processing” time required to traverse these states. This simplification allows us to succinctly describe kinetic proofreading in terms of a single waiting time parameter, facilitating a more direct evaluation and comparison of KPR performance across different biological contexts such as DNA replication and T cell receptor (TCR) signaling. We find that loss of fidelity for longer proofreading steps relies on the specific strategy of information extraction and show that in the first-passage time (FPT) discrimination strategy, longer proofreading steps can exponentially improve the accuracy of KPR at the cost of speed. Thus, KPR can still be an effective discrimination mechanism in the high noise regime. However, in a product concentration-based discrimination strategy, longer proofreading steps do not necessarily lead to an increase in performance. However, by introducing activation thresholds on product concentrations, can we decompose the product-based strategy into a series of FPT-based strategies to better resolve the subtleties of KPR-mediated product discrimination. Our findings underscore the importance of understanding KPR in the context of how information is extracted and processed in the cell.
Nanopores spanning synthetic membranes have been used as key components in proof-of-principle nanofluidic applications, particularly those involving manipulation of biomolecules or sequencing of DNA. ...The only practical way of manipulating charged macromolecules near nanopores is through a voltage difference applied across the nanopore-spanning membrane. However, recent experiments have shown that salt concentration gradients applied across nanopores can also dramatically enhance charged particle capture from a low concentration reservoir of charged molecules at one end of the nanopore. This puzzling effect has hitherto eluded a physically consistent theoretical explanation. Here, we propose an electrokinetic mechanism of this enhanced capture that relies on the electrostatic potential near the pore mouth. For long pores with diameter much greater than the local screening length, we obtain accurate analytic expressions showing how salt gradients control the local conductivity which can lead to increased local electrostatic potentials and charged analyte capture rates. We also find that the attractive electrostatic potential may be balanced by an outward, repulsive electro-osmotic flow that can in certain cases conspire with the salt gradient to further enhance the analyte capture rate.
We develop and analyze mathematical models for receptor-mediated transcytosis of monoclonal antibodies (MAb) targeting the transferrin receptor (TfR) or the insulin receptor (IR), which are expressed ...at the blood-brain barrier (BBB). The mass-action kinetic model for both the TfR and IR antibodies were solved numerically to generate predictions for the concentrations of all species in all compartments considered. Using these models, we estimated the rates of MAb endocytosis into brain capillary endothelium, which forms the BBB in vivo, the rates of MAb exocytosis from the intra-endothelial compartment into brain extracellular space, and the rates of receptor recycling from the endothelial space back to the luminal endothelial plasma membrane. Our analysis highlights the optimal rates of MAb association with the targeted receptor. An important role of the endogenous ligand, transferrin (Tf) or insulin, in receptor-mediated-transport (RMT) of the associated MAb was found and was attributed to the five order magnitude difference between plasma concentrations of Tf (25,000 nM) and insulin (0.3 nM). Our modeling shows that the very high plasma concentration of Tf leads to only 5% of the endothelial TfR expressed on the luminal endothelial membrane.
Efficient testing and vaccination protocols are critical aspects of epidemic management. To study the optimal allocation of limited testing and vaccination resources in a heterogeneous contact ...network of interacting susceptible, infected, and recovered individuals, we present a degree-based testing and vaccination model for which we derive optimal policies using control-theoretic methods. Within our framework, we find that optimal intervention policies first target high-degree nodes before shifting to lower-degree nodes in a time-dependent manner. Using such optimal policies, it is possible to delay outbreaks and reduce incidence rates to a greater extent than uniform and reinforcement-learning-based interventions, particularly on certain scale-free networks.
We consider age-structured models with an imposed refractory period between births. These models can be used to formulate alternative population control strategies to China's one-child policy. By ...allowing any number of births, but with an imposed delay between births, we show how the total population can be decreased and how a relatively older age distribution can be generated. This delay represents a more 'continuous' form of population management for which the strict one-child policy is a limiting case. Such a policy approach could be more easily accepted by society. Our analyses provide an initial framework for studying demographics and how social constraints influence population structure.
Abstract
Solving analytically intractable partial differential equations (PDEs) that involve at least one variable defined on an unbounded domain arises in numerous physical applications. Accurately ...solving unbounded domain PDEs requires efficient numerical methods that can resolve the dependence of the PDE on the unbounded variable over at least several orders of magnitude. We propose a solution to such problems by combining two classes of numerical methods: (i) adaptive spectral methods and (ii) physics-informed neural networks (PINNs). The numerical approach that we develop takes advantage of the ability of PINNs to easily implement high-order numerical schemes to efficiently solve PDEs and extrapolate numerical solutions at any point in space and time. We then show how recently introduced adaptive techniques for spectral methods can be integrated into PINN-based PDE solvers to obtain numerical solutions of unbounded domain problems that cannot be efficiently approximated by standard PINNs. Through a number of examples, we demonstrate the advantages of the proposed spectrally adapted PINNs in solving PDEs and estimating model parameters from noisy observations in unbounded domains.
The specificity of T cells is that each T cell has only one T cell receptor (TCR). A T cell clone represents a collection of T cells with the same TCR sequence. Thus, the number of different T cell ...clones in an organism reflects the number of different T cell receptors (TCRs) that arise from recombination of the V(D)J gene segments during T cell development in the thymus. TCR diversity and more specifically, the clone abundance distribution, are important factors in immune functions. Specific recombination patterns occur more frequently than others while subsequent interactions between TCRs and self-antigens are known to trigger proliferation and sustain naive T cell survival. These processes are TCR-dependent, leading to clone-dependent thymic export and naive T cell proliferation rates. We describe the heterogeneous steady-state population of naive T cells (those that have not yet been antigenically triggered) by using a mean-field model of a regulated birth-death-immigration process. After accounting for random sampling, we investigate how TCR-dependent heterogeneities in immigration and proliferation rates affect the shape of clone abundance distributions (the number of different clones that are represented by a specific number of cells, or "clone counts"). By using reasonable physiological parameter values and fitting predicted clone counts to experimentally sampled clone abundances, we show that realistic levels of heterogeneity in immigration rates cause very little change to predicted clone-counts, but that modest heterogeneity in proliferation rates can generate the observed clone abundances. Our analysis provides constraints among physiological parameters that are necessary to yield predictions that qualitatively match the data. Assumptions of the model and potentially other important mechanistic factors are discussed.
In a recent clone-tracking experiment, millions of uniquely tagged hematopoietic stem cells (HSCs) and progenitor cells were autologously transplanted into rhesus macaques and peripheral blood ...containing thousands of tags were sampled and sequenced over 14 years to quantify the abundance of hundreds to thousands of tags or "clones." Two major puzzles of the data have been observed: consistent differences and massive temporal fluctuations of clone populations. The large sample-to-sample variability can lead clones to occasionally go "extinct" but "resurrect" themselves in subsequent samples. Although heterogeneity in HSC differentiation rates, potentially due to tagging, and random sampling of the animals' blood and cellular demographic stochasticity might be invoked to explain these features, we show that random sampling cannot explain the magnitude of the temporal fluctuations. Moreover, we show through simpler neutral mechanistic and statistical models of hematopoiesis of tagged cells that a broad distribution in clone sizes can arise from stochastic HSC self-renewal instead of tag-induced heterogeneity. The very large clone population fluctuations that often lead to extinctions and resurrections can be naturally explained by a generation-limited proliferation constraint on the progenitor cells. This constraint leads to bursty cell population dynamics underlying the large temporal fluctuations. We analyzed experimental clone abundance data using a new statistic that counts clonal disappearances and provided least-squares estimates of two key model parameters in our model, the total HSC differentiation rate and the maximum number of progenitor-cell divisions.
To analyze the influence of ocular magnification on the peripapillary retinal nerve fiber layer (RNFL) thickness measurement and its performance as acquired with spectral-domain optical coherence ...tomography (SD-OCT).
Spectral domain OCT measurements from 108 normal eyes (59 subjects) and 72 glaucoma eyes (58 patients) were exported and custom software was used to correct RNFL measurements for ocular magnification. Retinal nerve fiber layer prediction limits in normal subjects, structure-function relationships, and RNFL performance for detection of glaucoma were compared before and after correction for ocular magnification (Bennett's formula). Association of disc area with cross-sectional RNFL area was explored.
The median (interquartile range, IQR) visual field mean deviation and scaling factor were 0 (-0.85 to 0.73) dB and 0.96 (0.93-1.00) in normal eyes and -4.0 (-6.0 to -2.2) dB and 0.99 (0.95-1.03) in the glaucoma group (P < 0.001 and P = 0.003, respectively; average correction 3%). Correction for ocular magnification caused a reversal of the negative relationship between the cross-sectional RNFL area and axial length (slope = -0.022 mm(2)/mm, P = 0.015 vs. = 0.22 mm(2)/mm, P = 0.007). However, such correction did not change RNFL prediction limits (except in superior and nasal quadrants), improve global or regional structure-function relationships, or enhance the ability of RNFL measurements to discriminate glaucoma from normal eyes (P > 0.05). The cross-sectional RNFL area was not correlated with optic disc area (P = 0.325).
Correction of RNFL measurements for ocular magnification did not improve prediction limits in normal subjects or enhance the performance of SD-OCT in this group of eyes in which the axial length did not deviate significantly from average values. The cross-sectional area of the RNFL was not related to the optic disc area.
Motivated by recent efforts by the criminal justice system to treat and rehabilitate nonviolent offenders rather than focusing solely on their punishment, we introduce an evolutionary game theoretic ...model to study the effects of "carrot and stick" intervention programs on criminal recidivism. We use stochastic simulations to study the evolution of a population where individuals may commit crimes depending on their past history, surrounding environment and, in the case of recidivists, on any counseling, educational or training programs available to them after being punished for their previous crimes. These sociological factors are embodied by effective parameters that determine the decision making probabilities. Players may decide to permanently reform or continue engaging in criminal activity, eventually reaching a state where they are considered incorrigible. Depending on parameter choices, the outcome of the game is a society with a majority of virtuous, rehabilitated citizens or incorrigibles. Since total resources may be limited, we constrain the combined punishment and rehabilitation costs per crime to be fixed, so that increasing one effort will necessarily decrease the other. We find that the most successful strategy in reducing crime is to optimally allocate resources so that after being punished, criminals experience impactful intervention programs, especially during the first stages of their return to society. Excessively harsh or lenient punishments are less effective. We also develop a system of coupled ordinary differential equations with memory effects to give a qualitative description of our simulated societal dynamics. We discuss our findings and sociological implications.