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 433 records · Page 24Linked to original sources

Mathematical models for the cellular concentrations of cyclin and MPF.

Several mathematical models have been proposed for regulation of the cell cycle in early embryos by cyclin and maturation-promoting factor (MPF). In this paper the previously proposed models for cyclin and MPF activity are analyzed, and the validity of those models based on the mathematical behavior of their solutions and on physical considerations are discussed. In addition, three further models are proposed that exhibit the periodic behavior necessary for modeling the mitotic clock but that do not have certain of the limitations of the other models.

Animals↗

Intraperitoneal carboplatin infusion may be a pharmacologically more reasonable route than intravenous administration as a systemic chemotherapy. A comparative pharmacokinetic analysis of platinum using a new mathematical model after intraperitoneal vs. intravenous infusion of carboplatin--a Sankai Gynecology Study Group (SGSG) study.

OBJECTIVE: To clarify the pharmacological advantage of carboplatin-based intraperitoneal chemotherapy using the three-compartment mathematical model. METHODS: Eleven consecutive patients in one institution underwent intraperitoneal administration of carboplatin, and 11 consecutive patients in another institution received intravenous administration. Carboplatin (AUC=6 mg x min/ml) was diluted in 500 ml 5% glucose and administered either as an intraperitoneal bolus infusion or intravenous drip infusion during 1 h. Patients undergoing intravenous injection also received an infusion of 500 ml 5% glucose to obtain intraperitoneal samples. Intraperitoneal fluid and blood samples were obtained, immediately and 1, 2, 4, 8, 12, and 24 h after administration. The mathematical model consisting of a three-compartment model was applied to analyze the pharmacokinetics. The model was created with simultaneous differential equations and was solved by the Runge-Kutta method. RESULTS: The rate constants of platinum diffusion from the peritoneal cavity to serum, serum to peritoneal cavity, serum to peripheral space, peripheral space to serum, and elimination were 0.94+/-0.79 (mean+/-SD), 1.28+/-2.50, 16.50+/-9.26, 0.99+/-0.62, and 4.14+/-1.45 (h-1), respectively. When the theoretical pharmacological concentration of platinum was calculated using this mathematical model, 24-h platinum AUC in the serum was exactly the same regardless of intraperitoneal or intravenous administration of carboplatin. However, the 24-h platinum AUC in the peritoneal cavity was approximately 17 times higher when carboplatin was administered by the intraperitoneal route. CONCLUSION: The present pharmacological analysis suggests that intraperitoneal infusion of carboplatin is feasible not only as an intraperitoneal regional therapy but also as a more reasonable route for systemic chemotherapy.

Antineoplastic Combined Chemotherapy Protocols↗

A mathematical model of the intracerebral steal phenomenon in regional and focal ischaemia.

The objective of the present work was to mathematically estimate the extent and dynamics of intracerebral steal which may occur in response to cerebral vasodilation in regional and focal cerebral ischaemia. To this end, a spatially distributed mathematical model of regional cerebral blood flow (rCBF) was developed. The model contained a parallel system of intracerebral vascular resistances which were connected in series to a lumped extracerebral artery resistance and, for the focal ischaemia model, also a lumped pial collateral resistance. The rCBF was measured at 30 min of ischaemia in the following models: (1) bilateral carotid occlusion in spontaneously hypertensive rats (SHR), and (2) occlusion of the middle cerebral artery (MCA) in normotensive rats. The measured 3-dimensional rCBF data were used to set up the initial values of intracerebral resistance components. Cerebral vasodilation induced by inhalation of CO2 was simulated in the model by decreasing the values of both intracerebral and collateral resistance. Vascular responsiveness was specified to decrease with the ischaemic rCBF. In addition, a long term change in rCBF and resistance distribution was introduced to account for: (1) gradual rise in intracerebral resistance due to ischaemic oedema, and (2) adaptive decrease in collateral resistance. The following were predicted by the mathematical model. (1) At 60% maximum intracerebral dilatation a small intracerebral steal (5-10%) occurs at flow levels below 30-50 ml/100 g/min in both ischaemic models. (2) In focal ischaemia, the steal can be compensated by the 5% to 20% decrease in the collateral vascular resistance. (3) The rate of collateral adaptation overcomes the rate of intracerebral resistance rise and, therefore, eliminates the intracerebral steal after an adequately long period of time (on the order of a few hours). (4) An inverse steal effect can be demonstrated at the end of vasodilatation, provided that the time constant of collateral adaptation selected is longer (about 5:1) than the time constant of the intracerebral resistance rise. We conclude that the prediction of rCBF response to vasodilatation in cerebral ischaemia requires a knowledge of resting rCBF and of the response characteristics of both intracerebral and pial arterial segments.

Animals↗

THE SQUID GIANT AXON. MATHEMATICAL MODELS.

The voltage clamp results of Hodgkin and Huxley have been reanalyzed in terms of alternative mathematical models. The model used for the potassium conductance changes is similar to that of the HH model except that an empirical functional relationship replaces the fourth power Law used by HH and the twenty-fifth power law used by Cole and Moore. The model used for the sodium conductance changes involves the explicit use of one variable only rather than the two variables m and h of HH. The rise and fall of the sodium conductance during a depolarizing voltage clamp is obtained by specifying that this one variable satisfies a second order differential equation which results from the coupling of two first order equations. Not only can the adjustable parameters of these models be made to give good fit to the clamp conductance data but the models can also then be used to compute action potential curves. Theoretical interpretations can also be given to these mathematical models.

Action Potentials↗

A non-linear mathematical model for the in vivo evaluation of the RES phagocytic function.

A new non-linear mathematical model was constructed in order to perform in vivo quantification of the RES phagocytic function. This method is based on the same technical facilities as used for the routine liver-spleen scintigraphy with radiocolloids [1, 2]. But kinetic modeling of dynamic Tc-99m-sulfur colloid data produced estimations of the functional RE-parameters: the clearance rate of the colloidal particles, the rate of phagocytosis, and the RES functional volume, which can not be obtained by classical approaches. This non-linear model was designed on the basis of the principal characteristics of particulate material interaction with macrophages (attachment, phagocytosis, digestion) [3, 4, 5]. The theoretically examined behavior of this in vivo mathematical model corresponds with the experimental behavior of the RES. The mathematical expression of the dynamics is the system of non-linear differential equations with constant coefficients that have no analytical solution. Fitting of the normalized heart blood time-activity curve was obtained to identify the unknown model parameters via non-linear regression. For this purpose general interactive PASCAL procedure IDPAR for a PDP-11/34 computer was used (an IBM PC version is also available). Two to three iterations were needed to estimate the set of unknown parameters for any patient study (1-1.5 min). A very good fitting was obtained between experimental and model curves in every case of different pathologies (error of the approximation is about 2-3%). Studies were performed using an in vivo bolus injection of 3.6 mg/80 kg commercially available colloid KOREN labeled with 3m-Ci 99m-Tc (analog of TCK-1). Our method was used to determine the RES functional parameters for patient groups with different levels of the RES dysfunction. Obtained results illustrate the possibilities of our technique to quantitatively estimate not only great pathology (portal cirrhosis), but also small changes of the RE-function (case of hyperlipidemia and ulcer gaster). In all patient groups marked changes of Tc-99m-sulfur colloid turnover were observed. In general, tracer clearance from the circulation was decreased, and the rate of phagocytosis and the RES volume were diminished compared with controls. The effect of a reduction of phagocytosis increases when the RES dysfunction becomes stronger. It can be shown that a non-parametric Wilcoxon-Mann-Whitney test gives a significant difference (P95%) for these patient groups. Further, we represent the possibility of using the model for monitoring changes of the RES-function parameters during and after therapy. The quantitative test of the RES function can significantly enhance the diagnosis and management of different diseases. Serial colloidal studies may document changes in the RES-function for the tumors, cirrhosis, hyperlipidemia, reticulosis, hepatitis, thrombosis, infection, AIDS, burn injury, shock and trauma patients. The technique may be useful for the different RES investigations with laboratory animals. Created computer software can be used as a tool for kinetic models, simulation, and unknown parameters identification.

Humans↗

A mathematical model for prediction of drug molecule diffusion across the blood-brain barrier.

BACKGROUND: Predicting the ability of drugs to enter the brain is a longstanding problem in neuropharmacology. The first step in creating a much-needed computational algorithm for predicting whether a drug will enter brain is to devise a rigorous mathematical model. METHODS: Employing two experimental measures of blood-brain barrier (BBB) penetrability (brain/plasma ratio and the brain-uptake index) and 14 theoretically derived biophysical predictors, a mathematical model was developed to quantitatively correlate molecular structure with ability to traverse the BBB. RESULTS: This mathematical model employs Stein's hydrogen bonding number and Randic's topological descriptors to correlate structure with ability to cross the BBB. The final model accurately predicts the ability of test molecules to cross the BBB. CONCLUSIONS: A mathematical method to predict blood-brain barrier penetrability of drug molecules has been successfully devised. As a result of bioinformatics, chemoinformatics and other informatics-based technologies, the number of small molecules being developed as potential therapeutics is increasing exponentially. A biophysically rigorous method to predict BBB penetrability will be a much-needed tool for the evaluation of these molecules.

Algorithms↗

A mathematical model for ocular tear and solute balance.

PURPOSE: In this paper we develop a mathematical model that can predict the steady-state tear film thickness and the dynamic tear film thickness and the solute concentration after instillation of a solute-laden fluid in the eye. METHODS: The mathematical model developed in this paper is based on a balance between the inflow and outflow of tears into the tear film. It incorporates a tear drainage model and a model that relates the tear film thickness to the meniscus radius of curvature. To predict the solute concentrations, the tear balance is coupled with the solute balance. The differential equations for the unsteady balances are solved numerically. RESULTS: The model predicts that the tear film thickness depends on a number of physiological factors, such as rates of tear production and evaporation, geometry and modulus of the canaliculi, and surface tension and viscosity of tears, and varies from about 3 to 15 microm. The model also predicts that the drainage time for an instilled volume of 15 microl is 1283 s. Additionally, the time required for the tracer concentration to decay to 1% of the value immediately after instillation of a drug-laden 40 microl drop is about 2480 s. Similarly, the time for intensity decay for a radioactive tracer after 25 microl instillation is about 1566 s. Also, the model predicts that the fraction of the instilled drug that reaches the cornea is about 1.3% for topical application of timolol. CONCLUSIONS: The predicted results agree reasonably with the reported experimental results, at least qualitatively. The model developed here can serve as a useful tool to develop a more precise understanding of various issues related to tears and also evaluate the effect of various parameters on the tear volume.

Eye↗

Mathematical modeling of wire-duct single-stage electrostatic precipitators.

A two-dimensional mathematical model was developed to simulate the performance of wire-duct single-stage electrostatic precipitators (ESP). The model presented by Talaie et al. [M.R. Talaie, M. Taheri, J. Fathikaljahi, A new method to evaluate the voltage-current characteristics applicable for a single-stage electrostatic precipitator, J. Electrostat., 53 (3) (2001) 221-233] was used for prediction of electric field strength distribution and V-I characteristic for high-voltage wire-plate configuration. Simple Lagrangian approach was used to predict particle movement. Normal k-epsilon turbulent flow model with considering electrical body force due to ion and charged particle flow was used to evaluate gas velocity distribution. Ignoring the effect of particle movement and fluid flow, the results of electrical part of mathematical model are in good agreement with experimental data of Penny and Matick [G.W. Penny, R.E. Matrick, Potential in DC corona field, Trans. AIEE Part 1, 79 (1960) 91-99]. The prediction of corona sheath radius and its variation with particle loading and applied voltage is the main distinguishing feature of the present model. This fact was not included in the earlier models.

Filtration↗

A new mathematical model of dynamic cerebral autoregulation based on a flow dependent feedback mechanism.

A new mathematical model representing dynamic cerebral autoregulation as a flow dependent feedback mechanism is presented. Two modelling parameters are introduced, lambda, the rate of restoration, and tau, a time delay. Velocity profiles are found for a general arterial blood pressure, allowing the model to be applied to any experiment that uses changes in arterial blood pressure to assess dynamic cerebral autoregulation. Two such techniques, thigh cuffs and a lower body negative pressure box, which produce step changes and oscillatory variations in arterial blood pressure respectively, are investigated. Results derived using the mathematical model are compared with data from the two experiments. The comparisons yield similar estimates for lambda and tau, suggesting these parameters are independent of the pressure change stimulus and depend only on the main features of the dynamic cerebral autoregulation process. The modelling also indicates that for imposed oscillatory variations in arterial blood pressure a small phase difference between pressure and velocity waveforms does not necessarily imply impaired autoregulation. It is shown that the ratio between the variation in maximum velocity and pressure variation can be used, along with the phase difference, to indicate the nature of the autoregulatory response.

Blood Flow Velocity↗

A mathematical model of pattern formation.

This paper presents an explicit mathematical model describing pattern formation in monolayer epithelia. The approach is a generalization of the equations describing soap bubble configurations (Plateau, 1873; Thompson, 1917; Almgren & Taylor, 1976) that allows adjacent cells to adhere with differing intensities (Steinberg, 1962, 1978). The model is a system of simultaneous non-linear equations that considers cell-cell interactions in a two-dimensional sheet. The implementation involves using the equations of the model to predict explicitly the energy-minimizing configuration of a system of cells, based on the adhesivity of their membranes. The model can thus be used to explore the effects of varying adhesions on the dynamics of pattern formation. Following Chichilnisky (1985), such a descriptive system is introduced in this paper, and its predictive properties explored.

Cell Aggregation↗

Virtual and real brain tumors: using mathematical modeling to quantify glioma growth and invasion.

Over the last 10 years increasingly complex mathematical models of cancerous growths have been developed, especially on solid tumors, in which growth primarily comes from cellular proliferation. The invasiveness of gliomas, however, requires a change in the concept to include cellular motility in addition to proliferative growth. In this article we review some of the recent developments in mathematical modeling of gliomas. We begin with a model of untreated gliomas and continue with models of polyclonal gliomas following chemotherapy or surgical resection. From relatively simple assumptions involving homogeneous brain tissue bounded by a few gross anatomical landmarks (ventricles and skull) the models have recently been expanded to include heterogeneous brain tissue with different motilities of glioma cells in grey and white matter on a geometrically complex brain domain, including sulcal boundaries, with a resolution of 1 mm(3) voxels. We conclude that the velocity of expansion is linear with time and varies about 10-fold, from about 4 mm/year for low-grade gliomas to about 3 mm/month for high-grade ones.

Brain Neoplasms↗

Estimation of the mtDNA mutation rate in aging mice by proteome analysis and mathematical modeling.

The accumulation of mitochondria containing mutated genomes was proposed to be an important factor involved in aging. Although the level of mutated mtDNA has shown to increase over time, it is currently not possible to directly measure the mtDNA mutation rate within living cells. The combination of mathematical modeling and controlled experiments is an alternative approach to obtain an estimate for the mutation rate in a well-defined system. In order to judge the relevance of mitochondrial mutations for the aging process, we used a mouse model to study age-related alterations of the mitochondrial proteins. Based on these experimental data we constructed a mathematical model of the mitochondrial population dynamics to estimate mtDNA mutation rates. Mitochondria were isolated from mouse brain and liver at six different ages (newborn to 24-months). A large-gel 2D-electrophoresis-based proteomics approach was used to analyze the mitochondrial proteins. The expression of two respiratory chain complex I subunits and one complex IV subunit decreased significantly with age. One subunit of complex III and one subunit of complex V increased in expression during aging. Together, these data indicate that complex I and IV deficiency in aged tissues might be accompanied by feedback regulation of other protein complexes in the respiratory chain. When we fitted our experimental data to the mathematical model, mtDNA mutation rate was estimated to be 2.7x10(-8) per mtDNA per day for brain and 3.2x10(-9) per mtDNA per day for liver. According to our model and in agreement with the mitochondrial theory of aging, mtDNA mutations could cause the detrimental changes seen in mitochondrial populations during the normal lifespan of mice, while at the same time ensure that the mitochondrial population remains functional during the developmental and reproductive period of mice.

Aging↗

[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.

Animals↗

Mathematical model of remission duration in acute myelogenous leukemia.

A mathematical model of disease-free survival in acute myelogenous leukemia is formulated in terms of proliferation of the leukemic stem cell compartment. A survivorship function is described which allows for the possibility that some patients may be cured of their disease and which captures the main features of experimental remission-duration curves. The model is used to investigate the effect of the amount of residual disease and the growth rate of the stem cell population on the length of remissions and the chances for cure. A method of comparing clinical trials, with potential for explanatory inference, is exemplified.

Cell Division↗

[Mathematical model of electromechanical coupling in the rat myocardium].

A mathematical modelling approach was used to study the negative staircase in the rat papillary muscle. This phenomenon was found to be associated with an excess of the steady-state Ca++-influx in the myocardium cells. The model with a single chamber intracellular Ca-pool simulates satisfactory the experimental data obtained with a standard set of parameter values. It is concluded that the rat myocardium sarcoplasmic reticulum is loaded by Ca++-ions which enter the cell as a potential-independent Ca-influx presumably.

Action Potentials↗

Mathematical models of the population biology of Ostertagia ostertagi and Teladorsagia circumcincta, and the economic evaluation of disease control strategies.

The construction and use of mathematical models of the population biology of Ostertagia ostertagi and Teladorsagia circumcincta is discussed. Simulated field trials implemented by deterministic mathematical models currently share with actual field trials the disadvantage that they convey no information concerning the risk associated with the net return demonstrated by the trial. This has important implications when it is necessary to rank disease control strategies in order of usefulness.

Animals↗

Mathematical modelling of 3-(3',4'-dichlorophenyl)-1,1-dimenthylurea action in plant leaves

A mathematical model of the action of a photosystem II herbicide 3-(3',4'-dichlorophenyl)-1, 1-dimenthylurea, DCMU, in plant leaves upon an external application is presented. The diffusion of DCMU in a plant tissue is described with the help of Fick's laws and the following reaction of the herbicide with the QB-binding site of photosystem II by the mass action theory. The model is used for a description of the effect of the herbicide on chlorophyll fluorescence induction (the O-J-I-P curve) measured with spring barley primary leaves submerged in the herbicide solution. The increase of the J step during the herbicide action is ascribed to an increase of the number of photosystem II centres with bound herbicide molecules and malfunctioning in the electron transport to the plastoquinone pool. The experimental data were fitted with the help of the mathematical model. Values of the diffusion coefficient and the second order rate constant of the reaction of the herbicide with photosystem II, obtained by the fitting procedure, are discussed.Copyright 1998 Academic Press Limited

Journal Article↗

Ratcheting in post-translational protein translocation: a mathematical model.

We have developed a non-steady-state mathematical model describing post-translational protein translocation across the endoplasmic reticulum membrane. Movement of the polypeptide chain through the channel in the endoplasmic reticulum membrane is considered to be a stochastic process which is biased at the lumenal side of the channel by the binding of BiP (Kar2p), a member of the Hsp70 family of ATPases (ratcheting model). Assuming that movement of the chain through the channel is caused by passive diffusion (Brownian ratchet), the model describes all available experimental data. The optimum set of model parameters indicates that the ratcheting mechanism functions at near-maximum rate, being relatively insensitive to variations of the association or dissociation rate constants of BiP or its concentration. The estimated rate constant for diffusion of a polypeptide inside the channel indicates that the chain makes contact with the walls of the channel. Since fitting of the model to the data required that the backward rate constant be larger than the forward constant during early diffusion steps, translocation must occur against a force. The latter may arise, for example, from the unfolding of the polypeptide chain in the cytosol. Our results indicate that the ratchet can transport polypeptides against a free energy of about 25 kJ/mol without significant retardation of translocation. The modeling also suggests that the BiP ratchet is optimized, allowing fast translocation to be coupled with minimum consumption of ATP and rapid dissociation of BiP in the lumen of the ER. Finally, we have estimated the maximum hydrophobicity of a polypeptide segment up to which lateral partitioning from the channel into the lipid phase does not result in significant retardation of translocation.

Adenosine Triphosphate↗