A mathematical model of the coupled mechanisms of cell adhesion,contraction and spreading

Recent research has shown that cell spreading is highly dependent on the contractility of its cytoskeleton and the mechanical properties of the environment it is located in. The dynamics of such process is critical for the development of tissue engineering strategy but is also a key player in wound contraction, tissue maintenance and angiogenesis. To better understand the underlying physics of such phenomena, the paper describes a mathematical formulation of cell spreading and contraction that couples the processes of stress fiber formation, protrusion growth through actin polymerization at the cell edge and dynamics of cross-membrane protein (integrins) enabling cell-substrate attachment. The evolving cell’s cytoskeleton is modeled as a mixture of fluid, proteins and filaments that can exchange mass and generate contraction. In particular, besides self-assembling into stress fibers, actin monomers able to polymerize into an actin meshwork at the cell’s boundary in order to push the membrane forward and generate protrusion. These processes are possible via the development of cell-substrate attachment complexes that arise from the mechano-sensitive equilibrium of membrane proteins, known as integrins. After deriving the governing equation driving the dynamics of cell evolution and spreading, we introduce a numerical solution based on the extended finite element method, combined with a level set formulation. Numerical simulations show that the proposed model is able to capture the dependency of cell spreading and contraction on substrate stiffness and chemistry. The very good agreement between model predictions and experimental observations suggests that mechanics plays a strong role into the coupled mechanisms of contraction, adhesion and spreading of adherent cells.     

A new mathematical model of the dynamics of an age cohort population

A new generalized conception of an organism is given. Based on this conception, a new mathematical model of ontogenesis of an individual and the survival of the age cohort of population was proposed. By using real data on the dynamics of the survival of the age cohort of population, the model enables one to determine the parameters characterizing the relationship man-environment in the context of survival and calculate the dynamics (from birth to death) of the model variables of the state of the organism.     

A theory of contraction and a mathematical model of striated muscle

A theory of contraction and an associated model of striated muscle are presented, based on the assumption that chemical energy is being converted into electrical energy which, in turn, is being converted into mechanical energy and heat.The model, set up for the frog sartorius muscle, is able to predict the “rowing” motion of the cross-bridges, the force-velocity relation, the tension-length curve, the isometric force, all energy rates (heat and work rates), the metabolic rates and all known features of the stretched, stimulated muscle (no ATP-splitting, stretching tension higher than isometric tension, etc.). It also offers an alternative explanation for Hill's thermoelastic effect. The significance of Hill's force-velocity equation in the context of this theory is also discussed in detail.     

Calcium dynamics and homeostasis in a mathematical model of the principal cell of the cortical collecting tubule

Calcium (Ca) dynamics are incorporated into a mathematical model of the principal cell in the cortical collecting tubule developed earlier in Strieter et al. (1992a. Am. J Physiol. 263:F1063-1075). The Ca components are modeled after the Othmer-Tang model for IP(3)-sensitive calcium channels (1993, in Experimental and Theoretical Advances in Biological Pattern Formation, 295-319). There are IP(3)-sensitive Ca channels and ATP-driven pumps on the membrane of the endoplasmic reticulum. Calcium enters the cell passively down its electrochemical gradient. A Ca pump and Na/Ca exchange in the basolateral membrane are responsible for the extrusion of cytoplasmic calcium. Na/Ca exchange can also operate in reverse mode to transport Ca into the cell. Regulatory effects of cytoplasmic Ca on the apical Na channels are modeled after experimental data that indicate apical Na permeability varies inversely with cytoplasmic Ca concentration. Numerical results on changes in intracellular Ca caused by decreasing NaCl in the bath and the lumen are similar to those from experiments in Bourdeau and Lau (1990. Am. J Physiol. 258:F1497-1503). This match of simulation and experiment requires the synergistic action of the Na/Ca exchanger and the Ca regulated apical Na permeability. In a homogeneous medium, cytoplasmic Ca becomes oscillatory when extracellular Na is severely decreased, as observed in experiments of cultured principal cells (Koster, H., C. van Os and R. Bindels. 1993. Kidney Int.43:828-836). This essentially pathological situation arises because the hyperpolarization of membrane potential caused by Na-free medium increases Ca influx into the cell, while the Na/Ca exchanger is inactivated by the low extracellular Na and can no longer move Ca out of the cell effectively. The raising of the total amount of intracellular Ca induces oscillatory Ca movement between the cytoplasm and the endoplasmic reticulum. Ca homeostasis is investigated under the condition of severe extracellular Ca variations. As extracellular Ca is decreased, Ca regulation is greatly impaired if Ca does not regulate apical ionic transport. The simulations indicate that the Na/Ca exchanger alone has only limited regulatory capacity. The Ca regulated apical sodium or potassium permeability are essential for regulation of cytoplasmic Ca in the principal cell of the cortical collecting tubule.     

A mathematical model for coupling within-host and between-host dynamics in an environmentally-driven infectious disease

This work presents a new model for the linking of within- and between-host dynamics. We use this as a conceptual model for the dynamics of Toxoplasma gondii, in which the parasite’s life cycle includes interactions with the environment. We postulate the infection process to depend on the size of the infective inoculum that susceptible hosts may acquire by interacting with a contaminated environment. Because the dynamical processes associated with the within- and between-host occur on different time scales, the model behaviors can be analyzed by using a singular perturbation argument, which allows us to decouple the full model by separating the fast- and slow-systems. We define new reproductive numbers for the within-host and between host dynamics for both the isolated systems and the coupled system. Particularly, the reproduction number for the between-host (slow) system dependent on the parameters associated with the within-host (fast) system in a very natural way. We show that these reproduction numbers determine the stability of the infection-free and the endemic equilibrium points. Our model may present a so-called backward bifurcation.     

A mathematical model for dynamics of CD40 clustering

Ligand bound-receptors in a signalosome complex trigger signals to determine cellular functions. Upon ligand binding, the ligand–receptor complexes form clusters on cell membrane. Guided by the previous experimental reports on the cluster formation of CD40, a trans membrane receptor for CD40-ligand, we built a minimal model of the receptor cluster formation. In this model, we studied co-operative and non-co-operative clustering of a maximum of four CD40 molecules assuming a positive mediator of clustering such as cholesterol to be present in both cases. We observed that co-operative interactions between CD40 molecules resulted in more of the largest CD40 clusters than that observed with the non-co-operatively interacting CD40 molecules. We performed global sensitivity analysis on the model parameters and the analyses suggested that cholesterol influenced only the initial stage of the co-operatively clustering CD40 molecules but it affects both the initial and the final stages in case of the non-co-operatively clustering CD40 molecules. Robustness analyses revealed that in both co-operative and non-co-operative interactions, the higher order clusters beyond a critical size are more robust with respect to alterations in the environmental parameters including the cholesterol. Thus, the role of co-operative and non-co-operative interactions in environment-influenced receptor clustering is reported for the first time.     

A mathematical model of blood,cerebrospinal fluid and brain dynamics

Using first principles of fluid and solid mechanics a comprehensive model of human intracranial dynamics is proposed. Blood, cerebrospinal fluid (CSF) and brain parenchyma as well as the spinal canal are included. The compartmental model predicts intracranial pressure gradients, blood and CSF flows and displacements in normal and pathological conditions like communicating hydrocephalus. The system of differential equations of first principles conservation balances is discretized and solved numerically. Fluid–solid interactions of the brain parenchyma with cerebral blood and CSF are calculated. The model provides the transitions from normal dynamics to the diseased state during the onset of communicating hydrocephalus. Predicted results were compared with physiological data from Cine phase-contrast magnetic resonance imaging to verify the dynamic model. Bolus injections into the CSF are simulated in the model and found to agree with clinical measurements.

---

     

A mathematical model of the dynamics of granulocytopoiesis in mammals

A mathematical model has been developed for the dynamics of granulocytopoiesis in mammals subjected to chronic irradiation. The model involves a chalones mechanism of haemopoiesis regulation and comprises 12 nonlinear differential equations. The simulation results agree with the experimental data concerning the dynamics of granulocytopoiesis in rats affected by radiation within a wide range of dose rates.     

A mathematical model of the dynamics of odontogenic cyst growth

OBJECTIVE: To formulate a mathematical model of odontogenic cyst growth and establish the dynamics of cyst enlargement and role of osmotic pressure forces throughout its growth. STUDY DESIGN: The model assumed a spherical cyst with a semipermeable lining of living cells and a core consisting of degraded cellular material, including generic osmotic material, fed by the continuous death of epithelial cells in the lining. The lining cells were assumed to have both elastic and viscous properties, reflecting the action of physical stresses by the surrounding cyst capsule, composed of fibroblasts and collagen fibers. The model couples the cyst radius and osmotic pressure differences resulting in a system of 2 nonlinear ordinary differential equations. RESULTS: The model predicts that in all parameter regimens the long-time behavior of the cyst is the same and that linear radial expansion results. CONCLUSION: In the early and intermediate stages of cystic growth, osmotic pressure differences play an important role; however, in very large cysts, this role becomes negligible, and cell birth in the lining dominates growth.     

A mathematical model for the population dynamics of army ants

A stochastic cellular automata model for the population dynamics of the army antEciton burchelli on Barro Colorado Island in Panama is set up. It is simulated on the computer and shown to give good agreement with biological data. It is analysed using two approximations akin to the mean field approximation in statistical mechanics, and good agreement with the simulations is obtained. Finally, the role of distance between successive statary phase bivouacs is discussed with regard to the rate of colony growth. There are two aspects of the biological system studied here that make it of general importance. First, the population is structured, since the size of each colony of army ants is crucial. Second, the spatial behaviour of the population, as in many others, is not diffusion-like, although it is random. This has implications for the kind of model that is chosen.     

A mathematical model of population dynamics for Batesian mimicry system

We analyse a mathematical model of the population dynamics among a mimic, a corresponding model, and their common predator populations. Predator changes its search-and-attack probability by forming and losing its search image. It cannot distinguish the mimic from the model. Once a predator eats a model individual, it comes to omit both the model and the mimic species from its diet menu. If a predator eats a mimic individual, it comes to increase the search-and-attack probability for both model and mimic. The predator may lose the repulsive/attractive search image with a probability per day. By analysing our model, we can derive the mathematical condition for the persistence of model and mimic populations, and then get the result that the condition for the persistence of model population does not depend on the mimic population size, while the condition for the persistence of mimic population does depend the predator's memory of search image.     

A mathematical model of population dynamics for Batesian mimicry system

We analyse a mathematical model of the population dynamics among a mimic, a corresponding model, and their common predator populations. Predator changes its search-and-attack probability by forming and losing its search image. It cannot distinguish the mimic from the model. Once a predator eats a model individual, it comes to omit both the model and the mimic species from its diet menu. If a predator eats a mimic individual, it comes to increase the search-and-attack probability for both model and mimic. The predator may lose the repulsive/attractive search image with a probability per day. By analysing our model, we can derive the mathematical condition for the persistence of model and mimic populations, and then get the result that the condition for the persistence of model population does not depend on the mimic population size, while the condition for the persistence of mimic population does depend the predator's memory of search image.     

A 3D skeletal muscle model coupled with active contraction of muscle fibres and hyperelastic behaviour

This paper presents a three-dimensional finite element model of skeletal muscle which was developed to simulate active and passive non-linear mechanical behaviours of the muscle during lengthening or shortening under either quasi-static or dynamic condition. Constitutive relation of the muscle was determined by using a strain energy approach, while active contraction behaviour of the muscle fibre was simulated by establishing a numerical algorithm based on the concept of the Hill's three-element muscle model. The proposed numerical algorithm could be used to predict concentric, eccentric, isometric and isotonic contraction behaviours of the muscle. The proposed numerical algorithm and constitutive model for the muscle were derived and implemented into a non-linear large deformation finite element programme ABAQUS by using user-defined material subroutines. A number of scenarios have been used to demonstrate capability of the model for simulating both quasi-static and dynamic response of the muscle. Validation of the proposed model has been performed by comparing the simulated results with the experimental ones of frog gastrocenemius muscle deformation. The effects of the fusiform muscle geometry and fibre orientation on the stress and fibre stretch distributions of frog muscle during isotonic contraction have also been investigated by using the proposed model. The predictability of the present model for dynamic response of the muscle has been demonstrated by simulating the extension of a squid tentacle during a strike to catch prey.     

A mathematical model for progression and heterogeneity in colorectal cancer dynamics

This paper deals with the development of a mathematical model that describes cancer dynamics at the cellular scale.The selected case study concerns colon and rectum cancer, which originates in colorectal crypts. Cells inside the crypts are assumed to be organized according to a compartmental-like arrangement and to be homogeneously mixing. A mathematical model for cancer progression is proposed here. This model describes the generation of multiple clonal sub-populations of cells at different progression stages in a single crypt.Asymptotic analysis and simulations are developed with an exploratory aim. The obtained results offer some insights into the role played by mutation, proliferation and differentiation phenomena on cancer dynamics. In particular, the acquisition of an additional growing power and a reduction for cellular differentiation seem more likely to be the driving force behind carcinogenesis rather than an increase in the mutation rate. The mutation rate instead seems to affect progression dynamics and intra-tumor heterogeneity. The role played by cells, at different differentiation stages, in the onset and progression of colorectal cancer is highlighted. The results support the fact that stem cells play a key role in cancer development and the idea that transit-amplifying cells could also take on an active role in carcinogenesis.     

A mathematical model of the response of principal cells of collecting ducts to hypotonic shock

The regulatory decrease in the volume of principal cells of collecting ducts to hypoosmotic shock has been investigated experimentally and using the mathematical modeling. A mathematical model of the response of collecting duct principal cells to hypotonic shock has been constructed on the basis of the experimental time course of changes in cell volume measured by the fluorescent dye Calcein. It was shown that the regulatory decrease in volume under hypotonic conditions occurs via a marked release of osmolytes and is accompanied by a decrease in water permeability of the cell membrane. The mathematical modeling of transmembrane transport processes allowed us to quantitatively estimate the changes in membrane water permeability, which decreased tenfold, from 2 x 10(-1) cm/s to 2 x 10(-2) cm/s. It was also shown that the effective regulatory decrease in the volume of collecting duct principal cells in hypotonic medium results from a significant increase in membrane permeability for K+, Cl-, and organic anions.