Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Mathematical Model”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54Linked to original sources

Mathematical model of a virus-neutralizing immunglobulin response.

We present a mathematical model to simulate the kinetics of B-cell activation and the virus-neutralizing immunoglobulin response in the spleen of mice after infection with vesicular stomatitis virus (VSV). Our model combines data from in vitro experiments and in vivo kinetic observations. A system of eight nonlinear differential equations was used in the computer experiments and numerically solved. The isotype switch from IgM to IgG in the presence of T-cell help was modelled by a time variable function, used as a parameter. The model solutions indicate fast kinetics of the generation of VSV-neutralizing IgM antibodies within 2-3 days post immunization peaking on day 5 at a serum concentration of approximately 80 microgram ml-1 IgM, which is the equivalent to about 10% of the total IgM serum concentration. The frequency of virus-specific B cells increases about 1000-fold within the first 4 days after immunization. Protective levels of VSV-neutralizing IgG antibodies (>/=10 microgram ml-1) are reached within 5 to 6 days post immunization. Fitting the model solution to the experimentally observed neutralizing serum titers suggests an increase in the neutralizing activity of IgGs occurring between days 5 and 8 post-infection. The model indicates that less than 10 VSV-specific B cells have to be triggered daily to maintain protective IgG serum titers during the memory phase.

Animals↗

A mathematical model of electron transfer within the mitochondrial respiratory cytochromes.

A simple mathematical model of electron flow along the mitochondrial respiratory cytochrome assembly and the transfer of electrons to molecular oxygen is presented. First, an expression for the current-voltage relationship for a biological oxygen electrode is derived, and from this the relationship between oxygen consumption rate and oxygen partial pressure is determined. An independent relationship between mitochondrial oxygen partial pressure and oxygen supply rate is then derived. By eliminating oxygen partial pressure from these two expressions, we may obtain a relationship between oxygen supply rate and oxygen consumption rate. This model is then used to investigate the effects of tissue dysoxia, uncoupling of oxidative phosphorylation, increased cellular diffusional resistance and inhomogeneities in oxygen supply on oxygen consumption. It is concluded that each of the above contribute in varying degrees to the phenomenon of "pathological oxygen supply dependency".

Critical Illness↗

A mathematical model for germinal centre kinetics and affinity maturation.

We present a mathematical model which reproduces experimental data on the germinal centre (GC) kinetics of the primed primary immune response and on affinity maturation observed during the reaction. We show that antigen masking by antibodies which are produced by emerging plasma cells can drive affinity maturation and provide a feedback mechanism by which the reaction is stable against variations in the initial antigen amount over several orders of magnitude. This provides a possible answer to the long-standing question of the role of antigen reduction in driving affinity maturation. By comparing model predictions with experimental results, we propose that the selection probability of centrocytes and the recycling probability of selected centrocytes are not constant but vary during the GC reaction with respect to time. It is shown that the efficiency of affinity maturation is highest if clones with an affinity for the antigen well above the average affinity in the GC leave the GC for either the memory or plasma cell pool. It is further shown that termination of somatic hypermutation several days before the end of the germinal centre reaction is beneficial for affinity maturation. The impact on affinity maturation of simultaneous initiation of memory cell formation and somatic hypermutation vs. delayed initiation of memory cell formation is discussed.

Antibody Affinity↗

A mathematical model for estimation of infarct size from serum creatine phosphokinase.

A mathematical model for kinetics of creatine phosphokinase (CK) after myocardial infarct is presented. The transport of CK from the infarct to the distribution space is modeled by a gamma variate. CK is released from the distribution space by a first-order process. The model permits simultaneous estimation of the disappearance rate parameter KD, the integrated appearance function, and the baseline from a set of observations of serum CK concentration.

Computers↗

Mathematical modeling of polymer-induced flocculation by charge neutralization.

A detailed mathematical model for flocculation of colloidal suspensions in presence of salts and polymers is described and validated. In former case, the classical DLVO theory, which accounts for relevant variables such as pH and salt concentration, is incorporated into a geometrically sectioned discrete population balance model. For processes involving polymers, flocculation via simple charge neutralization is modeled using a modified DLVO theory in which the effect of adsorbed polymer layers on van der Waals attraction is included. The fractal dimension of aggregates is obtained by dynamic scaling of experimental data for time evolution of mean aggregate size. The particle surface potential is assumed to be approximately equal to the zeta potential. The model predictions are in close agreement with experimental results for flocculation of colloidal hematite suspensions in the presence of KCl and polyacrylic acid at different concentrations. In particular, given values of model parameters, e.g., Hamaker constant, fractal dimension, surface potential, and thickness of adsorbed polymer layer, the model can realistically describe the kinetics of flocculation by a simple charge neutralization mechanism and track the evolution of floc size distribution. Representative examples of sensitivity of the flocculation model to perturbations in surface potential and fractal dimension and to modification in the DLVO theory for polymer-coated particles are included.

Journal Article↗

A mathematical model of action potential heterogeneity in adult rat left ventricular myocytes.

Mathematical models were developed to reconstruct the action potentials (AP) recorded in epicardial and endocardial myocytes isolated from the adult rat left ventricle. The main goal was to obtain additional insight into the ionic mechanisms responsible for the transmural AP heterogeneity. The simulation results support the hypothesis that the smaller density and the slower reactivation kinetics of the Ca(2+)-independent transient outward K(+) current (I(t)) in the endocardial myocytes can account for the longer action potential duration (APD), and more prominent rate dependence in that cell type. The larger density of the Na(+) current (I(Na)) in the endocardial myocytes results in a faster upstroke (dV/dt(max)). This, in addition to the smaller magnitude of I(t), is responsible for the larger peak overshoot of the simulated endocardial AP. The prolonged APD in the endocardial cell also leads to an enhanced amplitude of the sustained K(+) current (I(ss)), and a larger influx of Ca(2+) ions via the L-type Ca(2+) current (I(CaL)). The latter results in an increased sarcoplasmic reticulum (SR) load, which is mainly responsible for the higher peak systolic value of the Ca(2+) transient [Ca(2+)](i), and the resultant increase in the Na(+)-Ca(2+) exchanger (I(NaCa)) activity, associated with the simulated endocardial AP. In combination, these calculations provide novel, quantitative insights into the repolarization process and its naturally occurring transmural variations in the rat left ventricle.

Action Potentials↗

A mathematical model and computer simulation study of insulin receptor regulation.

A homeomorphic mathematical model of cell surface insulin receptor regulation is developed. The overall structure of the model is based on molecular mechanisms suggested by in vivo and in vitro experimental evidence from many different cell types. Model parameters correspond to cellular processes which are constrained by known boundry value conditions. As an example, computer simulation results are compared with published data from BC3H-1 myocytes in culture. With appropriate parameter choice, this model is able to simulate data from other cell types. Cellular processes which are explicitly represented in the model include: bound and unbound receptor endocytosis, receptor recycling, intracellular receptor degradation, and state-dependent receptor synthesis. Most of these processes are represented as first-order events. Using more complex representations of the model structure with higher order rate constants or saturable pathways does not qualitatively improve simulation results. Simulations are able to reproduce ligand-induced down and up regulation of receptors as well as the initial spontaneous display of surface insulin receptors. To demonstrate the behavior of our model and illustrate its utility for explaining insulin receptor regulation for a variety of conditions, simulations for which experimental data is unavailable for direct comparison are also shown. We believe the structure of our model is sufficient to explain insulin receptor regulation in a wide variety of cell types. In addition our model may aid in understanding the receptor component of insulin resistance (decreased sensitivity or responsiveness to insulin) seen in pathological states such as obesity and diabetes mellitus. Finally, this model may be applicable to the study of the regulation of other polypeptide hormone receptors.

Cell Membrane↗

A mathematical model of the dynamics of scrapie in a sheep flock.

A mathematical model is developed for the dynamics of an outbreak of scrapie in a single sheep flock with the aim of assisting the interpretation of field data. The model incorporates age structure of the sheep population, both horizontal and vertical transmission, genetic predisposition to infection, variable initial load of the infectious agent, and increasing infection load during an incubation period of the same order as sheep life expectancy. This leads to system of partial differential equations with respect to time, age and infection load. Numerical analyses using this model demonstrate that a scrapie outbreak is likely to be of long duration (several decades), will lead to reduction of scrapie susceptible allele frequency (but not to zero), and has different dynamics in homozygous and heterozygous susceptible sheep, even if these genotypes are equally susceptible, due to the different contributions of vertical infection to transmission to genotypes.

Age Factors↗

A mathematical model for crystal growth by aggregation of precursor metastable nanoparticles.

A mathematical model is developed to describe aggregative crystal growth, including oriented aggregation, from evolving pre-existing primary nanoparticles with composition and structure that are different from that of the final crystalline aggregate. The basic assumptions of the model are based on the ideas introduced in an earlier published report [Buyanov and Krivoruchko, Kinet. Katal. 1976, 17, 666-675] to describe the growth of low-solubility metal hydroxides (e.g., iron oxides) by oriented aggregation. It is assumed that primary particles can be described as pseudo-species A, B, and C, which have the following properties: (1) fresh primary particles (colloidally stable inert nanoparticles, denoted as A), (2) mature primary particles (partially transformed nanoparticles at an optimum stage of development for attachment to a growing crystal, denoted as B), and (3) nucleated primary particles (denoted as C1). The evolution of primary particles, A --> B --> C1, is treated as two first-order consecutive reactions. Crystal growth via crystal-crystal aggregation (Ci +) is described using the Smoluchowski equation. The new element of this model is the inclusion of an additional crystal growth mechanism via the addition of primary particles (B) to crystals (Ci): (B + ). Two distinct, but constant, kernels (K not equal K') are used. It is shown that, when K' = 0, a steplike crystal size distribution (CSD) is obtained. Within a range of K'/K values (e.g., K'/K = 10(3)), CSD with multiple peaks are obtained. Comparison with predictions of models that do not include the intermediate stage of primary particles (B) indicates pronounced differences. Despite its simplicity, the model is able to capture the qualitative features of CSD evolution that have been obtained from crystal growth experiments in hematite, which is a system that is believed to undergo oriented aggregation.

Journal Article↗

A mathematical model of the folate cycle: new insights into folate homeostasis.

A mathematical model is developed for the folate cycle based on standard biochemical kinetics. We use the model to provide new insights into several different mechanisms of folate homeostasis. The model reproduces the known pool sizes of folate substrates and the fluxes through each of the loops of the folate cycle and has the qualitative behavior observed in a variety of experimental studies. Vitamin B(12) deficiency, modeled as a reduction in the V(max) of the methionine synthase reaction, results in a secondary folate deficiency via the accumulation of folate as 5-methyltetrahydrofolate (the "methyl trap"). One form of homeostasis is revealed by the fact that a 100-fold up-regulation of thymidylate synthase and dihydrofolate reductase (known to occur at the G(1)/S transition) dramatically increases pyrimidine production without affecting the other reactions of the folate cycle. The model also predicts that an almost total inhibition of dihydrofolate reductase is required to significantly inhibit the thymidylate synthase reaction, consistent with experimental and clinical studies on the effects of methotrexate. Sensitivity to variation in enzymatic parameters tends to be local in the cycle and inversely proportional to the number of reactions that interconvert two folate substrates. Another form of homeostasis is a consequence of the nonenzymatic binding of folate substrates to folate enzymes. Without folate binding, the velocities of the reactions decrease approximately linearly as total folate is decreased. In the presence of folate binding and allosteric inhibition, the velocities show a remarkable constancy as total folate is decreased.

Allosteric Site↗

A mathematical model and a computerized simulation of PCR using complex templates.

A mathematical model and a computer simulation were used to study PCR specificity. The model describes the occurrences of non-targeted PCR products formed through random primer-template interactions. The PCR simulation scans DNA sequence databases with primers pairs. According to the model prediction, PCR with complex templates should rarely yield non-targeted products under typical reaction conditions. This is surprising as such products are often amplified in real PCR under conditions optimized for stringency. The causes for this 'PCR paradox' were investigated by comparing the model predictions with simulation results. We found that deviations from randomness in sequences from real genomes could not explain the frequent occurrence of non-targeted products in real PCR. The most likely explanation to the 'PCR paradox' is a relatively high tolerance of PCR to mismatches. The model also predicts that mismatch tolerance has the strongest effect on the number of non-targeted products, followed by primer length, template size and product size limit. The model and the simulation can be utilized for PCR studies, primer design and probing DNA uniqueness and randomness.

Animals↗

Analysis of growth of multicellular tumour spheroids by mathematical models.

We wished to determine the applicability of previously proposed deterministic mathematical models to description of growth of multicellular tumour spheroids. The models were placed into three general classes: empirical, functional and structural. From these classes, 17 models were applied systematically to growth curves of multicellular tumour spheroids used as paradigms of prevascular and microregional tumour growth. The spheroid growth curves were determined with uniquely high density of measurements and high precision. The theoretical growth curves obtained from the models were fitted by the weighted least-squares method to the 15 measured growth curves, each corresponding to a different cell line. The classical growth models such as von Bertalanffy, logistic and Gompertz were considered as nested within more general models. Our results demonstrate that most models fitted the data fairly well and that criteria other than statistical had to be used for final selection. The Gompertz, the autostimulation and the simple spheroid models were the most appropriate for spheroid growth in the empirical, functional and structural classes of models, respectively. We also showed that some models (e.g. logistic, von Bertalanffy) were clearly inadequate. Thus, contrary to the widely held belief, the sigmoid character of a three or more parameter growth function is not sufficient for adequate fits.

Cell Division↗

A mathematical model of the outer medullary collecting duct of the rat.

A mathematical model of the outer medullary collecting duct (OMCD) has been developed, consisting of alpha-intercalated cells and a paracellular pathway, and which includes Na(+), K(+), Cl(-), HCO(3)(-), CO(2), H(2)CO(3), phosphate, ammonia, and urea. Proton secretion across the luminal cell membrane is mediated by both H(+)-ATPase and H-K-ATPase, with fluxes through the H-K-ATPase given by a previously developed kinetic model (Weinstein AM. Am J Physiol Renal Physiol 274: F856-F867, 1998). The flux across each ATPase is substantial, and variation in abundance of either pump can be used to control OMCD proton secretion. In comparison with the H(+)-ATPase, flux through the H-K-ATPase is relatively insensitive to changes in lumen pH, so as luminal acidification proceeds, proton secretion shifts toward this pathway. Peritubular HCO(3)(-) exit is via a conductive pathway and via the Cl(-)/HCO(3)(-) exchanger, AE1. To represent AE1, a kinetic model has been developed based on transport studies obtained at 38 degrees C in red blood cells. (Gasbjerg PK, Knauf PA, and Brahm J. J Gen Physiol 108: 565-575, 1996; Knauf PA, Gasbjerg PK, and Brahm J. J Gen Physiol 108: 577-589, 1996). Model calculations indicate that if all of the chloride entry via AE1 recycles across a peritubular chloride channel and if this channel is anything other than highly selective for chloride, then it should conduct a substantial fraction of the bicarbonate exit. Since both luminal membrane proton pumps are sensitive to small changes in cytosolic pH, variation in density of either AE1 or peritubular anion conductance can modulate OMCD proton secretory rate. With respect to the OMCD in situ, available buffer is predicted to be abundant, including delivered HCO(3)(-) and HPO(4)(2-), as well as peritubular NH(3). Thus, buffer availability is unlikely to exert a regulatory role in total proton secretion by this tubule segment.

Acid-Base Equilibrium↗

A mathematical model for patient skin dose assessment in cardiac catheterization procedures.

A mathematical model has been developed for the assessment of patient skin doses from cardiac catheterization procedures. This uses exposure and projection data stored in the DICOM image files. Since these contain only information about the acquisition runs, a correction is needed to estimate and include the contribution from fluoroscopy. Maximum skin doses calculated by the model were found to correlate well with those measured on Kodak EDR2 film. Three methods for including the contribution from fluoroscopy were investigated, and all successfully identified patients receiving skin doses in excess of 1 Gy. It is hoped to automate this tool for routine assessment of skin doses in our cardiac catheterization laboratories.

Cardiac Catheterization↗

Mathematical modelling of the anaerobic hybrid reactor.

This paper presents a new mathematical model for the anaerobic hybrid reactor (AHR) (a UASB reactor and an anaerobic filter in series) and its experimental calibration and verification. The model includes a biochemical part and a mass transport one, which considers the AHR as two contact reactors in series. The anaerobic process transformations are described by the model developed by Siegrist et al. The fraction (F) of solids in the clarification zone of the UASB reactor that leaves this first reactor is the key physical parameter to be estimated. The main parameters of the model were calibrated using experimental results from a bench-scale AHR fed with real slaughterhouse wastewater. The fraction of inert particulate COD in the influent and the factor F were estimated by a trial and error procedure comparing experimental and simulated results of the mass of solids in the lower tank and the VSS concentration in the AHR effluent. A good fit was obtained. The final verification was carried out by comparing a set of experiments with simulated data. The model's capability to predict the process performance was thus proved.

Anaerobiosis↗

Effects of semantic cues on mathematical modeling: evidence from word-problem solving and equation construction tasks.

Mathematical solutions to textbook word problems are correlated with semantic relations between the objects described in the problem texts. In particular, division problems usually involve functionally related objects (e.g., tulips-vases) and rarely involve categorically related objects (e.g., tulips-daisies). We examined whether middle school, high school, and college students use object relations when they solve division word problems (WP) or perform the less familiar task of representing verbal statements with algebraic equations (EQ). Both tasks involved multiplicative comparison statements with either categorically or functionally related objects (e.g., "four times as many cupcakes [commuters] as brownies [automobiles]"). Object relations affected the frequency of correct solutions in the WP task but not in the EQ task. In the latter task, object relations did affect the structure of nonalgebraic equation errors. We argue that students use object relations as "semantic cues" when they engage in the sense-making activity of mathematical modeling.

Adolescent↗

A mathematical model of phytochrome - the receptor of photomorphogenetic processes in plants.

A mathematical model of the phytochrome pigment system is considered, the parameters of which are determined by solving the problem of identification on the basis of the dynamic characteristics of the photoreceptor obtained experimentally. The model is based on two photochemical and four dark reactions occurring in the phytochrome system. A study of the model permitted establishing the dynamics of processes in the photoreceptor and determining the relation between the behavior of the active form of phytochrome and the intensity and spectral distribution of light. The results of the investigation allows explaining of the high-energy reaction of photomorphogenesis as a consequence of the dynamics of phytochrome processes.

Darkness↗

Mathematical models for predicting straining G-level tolerances in reclined subjects.

Mathematical models that predict straining G-level tolerances of subjects reclined > 55 degrees are described. Straining G-level tolerances are defined as those G-levels of exposures that require the use of the anti-G straining maneuver (AGSM). A subject reclined > 55 degrees has a reduced inspiratory volume that limits the maximum intrathoracic pressure that can be developed that is the basis for the AGSM. Therefore the inability to perform a maximum AGSM reduces the maximum G-level that can be tolerated. G-level tolerances of subjects at seat-back angles between 75 degrees and 55 degrees are modeled. The use of positive pressure breathing (PBG) to increase G-level tolerance is considered. Maximum G-level tolerances are dependent upon anti-G suit pressures (psi) used and are predicted as follows: 55 degrees = 11.7 and 12.7 G with 7 and 10 psi; 65 dergrees = 12.1 and 12.8 G with 6 and 9 psi; and 75 degrees = 13.4 and 14.2 G with 5 and 8 psi anti-G suit pressure.

Acceleration↗