The size of a cell is central to many functions, including cellular communication and exchange of materials with the environment. This modeling and experimental study focused on understanding how the size of a cell determines its ability to uptake nanometer-scale extracellular materials from the environment. Several mechanisms in the cell plasma membrane mediate cellular uptake of nutrients, biomolecules, and particles. These mechanisms involve recognition and internalization of the extracellular molecules via endocytic components, such as clathrin-coated pits, vacuoles, and micropinocytic vesicles. Because the demand for an external resource could be different for cells of different sizes, the collective actions of these various endocytic routes should also vary based on the cell size. Here, using a reaction-diffusion model, we analyze single-cell data to interrogate the one/one mapping between the size of the MDA-MB 231 breast cancer cells and their ability to uptake nanoparticles. Our analysis indicates that under both reaction-and diffusion-controlled regimes, cellular uptake follows a linear relationship with the cell radius. Furthermore, this linear dependency is insensitive to particle size variation within 20-200 nm range. This result is counterintuitive because the general perception is that cellular uptake is proportional to the cell volume (mass) or surface area and hence follow a cubic or square relationship with the cell radius. A further analysis using our model reveals a potential mechanism underlying this linear relationship.

No abstract

Receptor tyrosine kinases (RTKs) typically contain multiple autophosphorylation sites in their cytoplasmic domains. Once activated, these autophosphorylation sites can recruit downstream signaling proteins containing Src homology 2 (SH2) and phosphotyrosine-binding (PTB) domains, which recognize phosphotyrosine-containing short linear motifs (SLiMs). These domains and SLiMs have polyspecific or promiscuous binding activities. Thus, multiple signaling proteins may compete for binding to a common SLiM and vice versa. To investigate the effects of competition on RTK signaling, we used a rule-based modeling approach to develop and analyze models for ligand-induced recruitment of SH2/PTB domain-containing proteins to autophosphorylation sites in the insulin-like growth factor 1 (IGF1) receptor (IGF1R). Models were parameterized using published datasets reporting protein copy numbers and site-specific binding affinities. Simulations were facilitated by a novel application of model restructuration, to reduce redundancy in rule-derived equations. We compare predictions obtained via numerical simulation of the model to those obtained through simple prediction methods, such as through an analytical approximation, or ranking by copy number and/or KD value, and find that the simple methods are unable to recapitulate the predictions of numerical simulations. We created 45 cell line-specific models that demonstrate how early events in IGF1R signaling depend on the protein abundance profile of a cell. Simulations, facilitated by model restructuration, identified pairs of IGF1R binding partners that are recruited in anti-correlated and correlated fashions, despite no inclusion of cooperativity in our models. This work shows that the outcome of competition depends on the physicochemical parameters that characterize pairwise interactions, as well as network properties, including network connectivity and the relative abundances of competitors.

This article introduces a multiscale framework for spatiotemporal modeling of protein-protein interaction. Cellular protein molecules represent multivalent species that contain modular features, such as binding domains and phosphorylation motifs. The binding and transformations of these features occur at a small time and spatial scale. On the other hand, space and time involved in protein diffusion, colocalization, and formation of complexes could be relatively large. Here, we present an agent-based framework integrated with a multiscale Brownian Dynamics (BD) simulation algorithm. The framework employs spatial graphs to describe multivalent molecules and complexes with their site-specific details. By implementing a time-adaptive feature, the BD algorithm enables efficient computation while capturing the site-specific interactions of the diffusing species at the sub-nanometer scale. We demonstrate these capabilities by modeling two multivalent molecules, one representing a ligand and the other a receptor, in a two-dimensional plane (cell membrane). Using the model, we show that the algorithm can accelerate computation by orders of magnitudes in both concentrated and dilute regimes. We also show that the algorithm enables robust model predictions against a wide range of selection of time step sizes.

Individual cells in a solution display variable uptake of nanomaterials, peptides, and nutrients. Such variability reflects their heterogeneity in endocytic capacity. In a recent work, we have shown that the endocytic capacity of a cell depends on its size and surface density of endocytic components (transporters). We also demonstrated that in MDA-MB-231 breast cancer cells, the cell-surface transporter density (n) may decay with cell radius (r) following the power rule n $ r a , where a z À1. In this work, we investigate how n and r may independently contribute to the endocytic heterogeneity of a cell population. Our analysis indicates that the smaller cells display more heterogeneity because of the higher stochastic variations in n. By contrast, the larger cells display a more uniform uptake, reflecting less-stochastic variations in n. We provide analyses of these dependencies by establishing a stochastic model. Our analysis reveals that the exponent a in the above relationship is not a constant; rather, it is a random variable whose distribution depends on cell size r. Using Bayesian analysis, we characterize the cell-size-dependent distributions of a that accurately capture the particle uptake heterogeneity of MDA-MB-231 cells.

Aggregation of cell surface receptor proteins by multivalent antigens is an essential early step for immune cell signalling. A number of experimental and modelling studies in the past have investigated multivalent ligand-mediated aggregation of IgE receptors (FcɛRI) in the plasma membrane of mast cells. However, understanding of the mechanisms of FcɛRI aggregation remains incomplete. Experimental reports indicate that FcɛRI forms relatively small and finite-sized clusters when stimulated by a multivalent ligand. By contrast, modelling studies have shown that receptor cross-linking by a trivalent ligand may lead to the formation of large receptor superaggregates that may potentially give rise to hyperactive cellular responses. In this work, we have developed a Brownian dynamics-based spatio-temporal model to analyse FcɛRI aggregation by a trivalent antigen. Unlike the existing models, which implemented non-spatial simulation approaches, our model explicitly accounts for the coarse-grained site-specific features of the multivalent species (molecules and complexes). The model incorporates membrane diffusion, steric collisions and sub-nanometre-scale site-specific interaction of the time-evolving species of arbitrary structures. Using the model, we investigated temporal evolution of the species and their diffusivities. Consistent with a recent experimental report, our model predicted sharp decay in species mobility in the plasma membrane in response receptor cross-linking by a multivalent antigen. We show that, due to such decay in the species mobility, post-stimulation receptor aggregation may become self-limiting. Our analysis reveals a potential regulatory mechanism suppressing hyperactivation of immune cells in response to multivalent antigens.

No abstract