PubMed HealthSearch

SEARCH · PubMed Health

Results for “Computer Simulation”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original PubMed record.

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

At least 127 records · Page 7Linked to original sources

Computer simulation of optimum personnel assignment in hospital pharmacy using a work-sampling method.

A computer simulation of the dispensing work performed in a hospital pharmacy was undertaken, based on analysis data from a work-sampling method and utilizing a personal computer. As a result it was found to be possible to estimate the complete work volume of a day based on the number of prescriptions in the day, and hence to predict fairly precisely the number of pharmacists required to complete that work. The method also enables the users to assign the optimum number of staff for dispensing duties according to the estimated number of prescriptions and the optimum waiting time of patients. A correlation has been noted between the total residence time of prescriptions and the mean waiting time of patients. Using this simple and effective technique to efficiently assign personnel makes it possible to cope with an increasing workload and limited number of staff.

Computer Simulation

Programmatic computer simulation model for medical school planning.

A comprehensive programmatic computer simulation planning model for a school of medicine was generated by integrating several separate simulation models with on-site cost study information. An elementary validation of the model was achieved. The model generated program costs in terms of both faculty hours and dollars. Results indicated that the size of the medical class could be increased from 75 to 100 students within the present resource limitations by transferring faculty time to education from other programs. The maximum class size was limited by the availability of clinical material. The basic science departments could handle this class size easily without significant reduction in other programs, but the clinical departments could not do so unless inpatient levels increased significantly.

Computers

Transformation of a polygonal cellular pattern during sexual maturation of the avian oviduct epithelium: computer simulation.

A peculiar cellular pattern resembling a checkerboard has been observed on the luminal surface of the oviduct epithelium of an adult Japanese quail. The epithelium is a monolayer cell sheet and consists of two types of columnar cells, ciliated cells (C-cells) and gland cells (G-cells) assembled in alternating blocks. The pattern develops, during sexual maturation, from a kagome-like pattern (in which large C-cells are surrounded by small G-cells) characteristic of the immature oviduct. In the present paper, computer simulations of the pattern transformation from kagome to checkerboard were performed assuming a few properties of individual cells. The adult checkerboard-like pattern is not strictly rectangular, but is deformed toward a honeycomb pattern. In theoretical considerations the assumption that adhesion is stronger between unlike cells than between like cells formed an ideal checkerboard pattern, because all cell boundaries in it are edges along which unlike cells meet. On the other hand, a honeycomb pattern formed after assuming that the boundary length of the cellular pattern is minimized (caused by contraction of bundles of microfilaments running along lateral boundaries of the columnar epithelial cell while keeping contact between neighbouring cells). The actual checkerboard-like pattern was considered to be in a balanced state between the effects of (1) the strong adhesion between unlike cells, and (2) the boundary contraction. Using a computational analysis, this consideration enabled us to obtain a quantitative parameter value for the difference between cell adhesions of unlike cells and of like cells. C-cells divided once during the kagome-checkerboard transformation, while G-cells did not divide. We performed computer simulations starting with the kagome pattern in which all C-cells divided once. The computer program of the boundary shortening procedure we used involved the quantitative parameter value for differential cell adhesion obtained as described above. A checkerboard pattern was successfully generated in the simulation. It is concluded that the strong adhesion between unlike cells and the boundary shortening have important roles in formation and maintenance of the kagome and checkerboard patterns of the avian oviduct epithelium.

Animals

Simulated bipolar cells in fovea of human retina. I. Computer simulation.

This static bipolar cell (BC) model of the human fovea is based on a number of reasonable assumptions. The human fovea is directly responsible for visual acuity and color vision. The fovea can be considered as having two parts; a central fovea with only red- and green-sensitive cones and a parafovea with blue-sensitive cones added to the other two. A cone mosaic can be precisely organized spatially into unit hexagons that specify inputs to horizontal cells (HC) and BCs. The retina up to and including BCs is piece-wise linear, i.e. at a given steady-state adapting light intensity BC outputs are linear functions of the physical image. BC centers receive inputs directly from weighted cones, while antagonistic surrounds receive inverted inputs from HCs. Appropriate optical and chromatic filtering due to the eye that are taken from human data are incorporated into the model. Chromatic aberrations are simulated by three separate point spread functions that also are taken from human data. Automatic gain control of cones is a function of intensity and wavelength of the steady adapting light.

Color Perception

Testing models of insulin binding in rat adipocytes using network thermodynamic computer simulations.

Many different models have been proposed to explain the complex binding behavior of insulin to its receptor, but a systematic comparison of models with experimental results is lacking. We have used network thermodynamic computer simulations to compare models of insulin binding against the results of several experimental tests designed to differentiate between the models. Six models of insulin binding were tested (simple, diffusion-reaction, conversion, dissociation, heterogeneous site, and two-step intramembrane) against results reported in the literature for isolated rat adipocytes. Although still a matter of experimental controversy, the criteria selected for modeling were curvilinear Scatchard plots, bi-or multi-exponential dissociation, insulin-accelerated dissociation, lack of dependence of the overall dissociation constant on receptor number, and receptor reserve. Using a given set of parameter values most appropriate for each model, none was able to account for all of the observed experimental results. This indicates both the complexity of the binding reaction and the need for further model development. The approach of using computer simulations to systematically test models against experimental results affords not only insight into the critical features of a model enabling it to pass a test, but also indicates potential experiments which might differentiate between models.

Adipose Tissue

Compliance in microcomputer-assisted conventional insulin therapy: computer simulation study results.

Compliance in diabetes self-management is a complex issue. It involves the interdependent daily actions of self-measurement of blood glucose and adherence to a prescribed schedule of daily activities. This impacts strongly on lifestyle because it necessitates precise meal timing as well as control of size and carbohydrate content. We sought to identify how strongly relaxing the lifestyle constraints per se would impact the ability to achieve improved metabolic control. To isolate these effects from those that result from poor measurement compliance, we used a computer simulator called OMNI et al. Furthermore, to standardize the "clinical" therapy, a second microprocessor device called an "Insulin Dosage Computer" was used to adjust insulin doses based on the usual clinical practice of four times a day precibal blood glucose measurements. Ten type 1 diabetic patients were stimulated and each followed for 120 simulated days. In each such subject, the simulation was repeated three times to include three different levels of lifestyle compliance ranging from excellent to poor. In all three protocols, starting from a level of poor control of diabetes, mean blood glucose values were significantly improved without significant differences after 120 days of computer-simulated treatment. Only the standard deviations, expressing the fluctuations of blood glucose and hence its stability, increased with decreasing lifestyle compliance. This computer simulation predicts that consistent self-monitoring of four blood glucose values per day is the cornerstone of diabetic self-control and that the use of these data according to a standardized therapeutic algorithm for insulin adjustment may successfully stabilize even patients with poor lifestyle compliance. Clinical studies must follow.

Blood Glucose

Reduction of community alcohol problems: computer simulation experiments in three counties.

A series of alcohol abuse prevention strategies was evaluated using computer simulation for three counties in the United States: Wake County, North Carolina, Washington County, Vermont and Alameda County, California. A system dynamics model composed of a network of interacting variables was developed for the pattern of alcoholic beverage consumption in a community. The relationship of community drinking patterns to various stimulus factors was specified in the model based on available empirical research. Stimulus factors included disposable income, alcoholic beverage prices, advertising exposure, minimum drinking age and changes in cultural norms. After a generic model was developed and validated on the national level, a computer-based system dynamics model was developed for each county, and a series of experiments was conducted to project the potential impact of specific prevention strategies. The project concluded that prevention efforts can both lower current levels of alcohol abuse and reduce projected increases in alcohol-related problems. Without such efforts, already high levels of alcohol-related family disruptions in the three counties could be expected to rise an additional 6% and drinking-related work problems 1-5%, over the next 10 years after controlling for population growth. Of the strategies tested, indexing the price of alcoholic beverages to the consumer price index in conjunction with the implementation of a community educational program with well-defined target audiences has the best potential for significant problem reduction in all three counties.

Adolescent

Computer simulations of N-methyl-D-aspartate receptor-induced membrane properties in a neuron model.

1. To evaluate the role of N-methyl-D-aspartate (NMDA) receptors in simulations of the lamprey spinal locomotor network, we developed a computer-simulated electrical model of a neuron that contains NMDA channels in addition to voltage-gated Na+, K+, and Ca2+ channels and Ca(2+)-activated K+ channels [K(Ca) channels]. 2. The voltage dependence of the Mg2+ block of the Na(+)-K+ current flow through the NMDA channel was modeled according to a scheme of open-channel block. To account for the regulation of K(Ca) channels by NMDA and membrane voltage, we modeled two separate Ca2+ pools that had different voltage dependencies and dynamics. 3. Pacemaker-like membrane potential oscillations could be elicited in the model neuron, which resembled those observed experimentally in the presence of bath-applied NMDA and tetrodotoxin. The effect of changing different channel parameters were tested to determine under which conditions such membrane potential oscillations could occur. 4. The oscillation amplitude was determined by the potential levels at which the NMDA channels and voltage-dependent K+ channels, respectively, were activated. The oscillation frequency and the relative durations of the de- and hyperpolarized phases of the oscillations were determined by the balance between the depolarizing (NMDA channels) and hyperpolarizing [K(Ca) channels] currents. 5. Simulated alterations of the Mg2+ concentration and the K+ conductance as well as injection of constant current caused changes of the oscillations corresponding to those observed experimentally. The de- and hyperpolarizing phases could be reset by brief current pulses. 6. We conclude that the present model can account for the effects of bath-applied NMDA on spinal neurons. This permits an incorporation of NMDA-receptor-mediated properties in simulation models of the lamprey locomotor network.

Action Potentials

Blood glucose control by intermittent loop closure in the basal mode: computer simulation studies with a diabetic model.

A semiclosed loop, bedside insulin infusion system using a simple basal infusion algorithm consisting of a linear transition between two insulin delivery rates as blood glucose (BG) increases has been developed. A theoretical study using computer simulation has now been undertaken to examine the effect of BG sampling frequency and algorithm parameters on BG control. A model for BG control by exogenous insulin in the individual with diabetes was developed from a model for healthy subjects and from clinical data in the literature. Results of computer simulation using this model showed a decrease in BG stability as the sampling interval increased from 1 to 4 h. Simulations also showed a decrease in BG stability as the sensitivity of the control algorithm increased. Choice of an appropriate basal control algorithm involved a compromise between stability, sampling interval, and metabolic control. We conclude that satisfactory metabolic control can be obtained using intermittent BG sampling in the basal state; sampling at intervals of 3 h combined with a basal control algorithm whereby insulin delivery rate increases linearly from 0.5 to 2.5 U/h over the BG range 2-12 mmol/L appears suitable for most diabetic persons. Three-hour sampling offers a good compromise between degree of metabolic control and clinical effort involved.

Blood Glucose

Computer simulation of biological systems. Current trends.

The current status of mathematical models of biological systems is reviewed. Advances in supercomputer hardware allows more complex models to be constructed. The new generation of microcomputers are quite adequate for many computer simulations of biological systems. A theory of modeling is being developed to improve the relationship between the real biological system and the model. Deterministic models, stochastic models and applications of control theory and optimization methods are discussed. Examples given include models of molecular structure, of experimental techniques, and of biochemical reactions. It is recommended that experimental biologists consider the use of microcomputers to model the system under study as a part of their research program.

Computer Simulation

Cellular sodium transport in arterial wall--summary of studies by digital computer simulation.

The distribution of electrolytes between extracellular and intracellular compartments is of greater functional significance in biological tissues than the total amount. A method commonly employed for the study of a substance distribution among different spaces is the continuous outflow and radioisotope counting of a fragment of tissue previously incubated with a radioactive tracer (here 22Na+). The experimental data can be analyzed within the frame of a compartmental analysis system which is defined by transport rate constants (kij) and solved when all kij are numerically determined. In practical applications the direct solution requires measurements in different compartments or knowledge of the initial condition for each compartment, but this is not possible in most biological tissues where the experimentally observable data represent the summation of tracer in all compartments at each point in time. The availability of the SAAM (Simulation, Analysis And Modelling) computer program allows for a powerful simulation procedure to arrive at values for kij. Experimental data and specifications for a chosen model are entered in the SAAM format. A simulated 'inflow' experiment is made by the computer by assigning to kij some estimate values. All initial conditions are now zero, since at the beginning there is no radioactive tracer in any compartment. After a time equivalent in the experimental realm to achieving constant specific activity, a 'time change' programmed into the computer takes place so that the outflow part of the experiment is developed with the same kij as for the inflow part, the final conditions for the inflow before the time change being the initial conditions for the outflow. An iterative procedure modifies the initial estimates of kij to seek the best fit. By coupling a powerful simulation procedure to an advanced radionuclide tracer biological methodology, we have been able to find effects of the hormone aldosterone and of experimental hypertension on Na+ distribution in arterial wall years before such effects were demonstrated by purely conventional experimental techniques. This computer simulation procedure may be suited to other experimental situations in which outflow data are the only available.

Aldosterone

Action potential refractory period in axonal demyelination: a computer simulation.

Axonal demyelination leads to an increase in the refractory period for propagation of the action potential. Computer simulations were used to investigate the mechanism by which changes in the passive properties of the internodal membrane increase the refractory period. The properties of the voltage dependent ion channels can be altered to restore conduction in demyelinated nerve fibers. The ability of these alterations to decrease the refractory period of demyelinated model nerve fibers was compared. The model nerve fiber contained six nodes. The action potential was stimulated at node one and propagated to node six. The internode between nodes three and four was demyelinated in a graded manner. The absolute refractory period for propagation of the action potential through the demyelinated internode increased as the number of myelin wraps was reduced to less than 25% of the normal value. The increase in refractory period was found to be due to a reduction in the rate or repolarization of the action potential at node three. The delay in repolarization reduced the rate of recovery of inactivated Na channels and slowed the closing of K channels. The rate of repolarization of node three was reduced by the conduction delay for the depolarization of node four caused by demyelination of the preceding internode. In these simulations the increase in refractory period due to demyelination was eliminated by slowing the onset of Na channel inactivation. A small reduction of the K conductance also decreased the refractory period. However, larger reductions eliminated this effect.

Action Potentials

Computer simulation as a tool for clinical trial design.

Use of a computer simulation technique as a tool for optimising the design of a clinical trial is described. One of several aspects of trial design facilitated by use of the technique is the estimation of the number of patients required in the trial. The simulation tool provides the basis for structured and disciplined discussion of all the issues to be considered in designing a trial. Copies of the simulation program are available from the author.

Clinical Trials as Topic

Computer simulations of cardiac action potentials in two dimensions.

A novel method of calculation is presented which allows computer simulations of action potentials for a two-dimensional, planar reconstruction of the depolarization phase of a cardiac action potential to be followed through nonuniform tissue. The calculation is explicit in type and assumes an infinite, grounded extracellular region. A 70 by 70 point region is calculated swiftly on a PC/AT. Results show elliptical isochrones when slow and fast directions have different resistivities. The time derivative of Vm is, however, similar in both directions. Insertion of a test region of tissue with high resistance results in wave-like propagation of the A.P. around the 'bad' region.

Action Potentials

Investigation of parameter estimator and adaptive controller for assist pump by computer simulation.

The multi-output adaptive controller of a left ventricular assist device (LVAD) was studied by computer simulation. The controller regulated two outputs--mean aortic pressure (mAoP) and mean atrial pressure (mLAP)--by regulating vacuum pressure (input). The autoregressive models were used to describe the circulatory system. The parameters of the models were estimated by the recursive least squares method. Based on the autoregressive models, the vacuum pressure minimizing a performance index was searched. The index used was the weighted summation of the square errors. Responses of the adaptive controller were simulated when the contractility of the left ventricle was decreased at various rates and the peripheral resistance was changed. Both the mAoP and mLAP were controlled to their predicted values in the steady state. The steady-state errors of the mAoP were less than a few mm Hg, and those of the mLAP were lower than 1 mm Hg. Consequently, the estimated parameters can be regarded as true parameters, and the adaptive controller has the potential to control more than two outputs. The multioutput adaptive controller studied is useful in controlling the LVAD according to the change in circulatory condition.

Aorta, Thoracic

Computer simulation of the segmental neural network generating locomotion in lamprey by using populations of network interneurons.

Realistic computer simulations of the experimentally established local spinal cord neural network generating swimming in the lamprey have been performed. Populations of network interneurons were used in which cellular properties, like cell size and membrane conductance including voltage dependent ion channels were randomly distributed around experimentally obtained mean values, as were synaptic conductances (kainate/AMPA, NMDA, glycine) and delays. This population model displayed more robust burst activity over a wider frequency range than the more simple subsample model used previously, and the pattern of interneuronal activity was appropriate. The strength of the reciprocal inhibition played a very important role in the regulation of burst frequency, and just by changing the inhibitory bias the entire physiological range could be covered. At the lower frequency range of bursting the segmental excitatory interneurons provide stability as does the activation of voltage dependent NMDA receptors. Spike frequency adaptation by means of summation of afterhyperpolarization (AHP) serves as a major burst terminating factor, and at lower rates the membrane properties conferred by the NMDA receptor activation. The lateral interneurons were not of critical importance for the burst termination. They may, however, be of particular importance for inducing a rapid burst termination during for instance steering and righting reactions. Several cellular factors combine to provide a secure and stable motor pattern in the entire frequency range.

Animals

Supercoiled DNA energetics and dynamics by computer simulation.

A new formulation is presented for investigating supercoiled DNA configurations by deterministic techniques. Thus far, the computational difficulties involved in applying deterministic methods to supercoiled DNA studies have generally limited computer simulations to stochastic approaches. While stochastic methods, such as simulated annealing and Metropolis-Monte Carlo sampling, are successful at generating a large number of configurations and estimating thermodynamic properties of topoisomer ensembles, deterministic methods offer an accurate characterization of the minima and a systematic following of their dynamics. To make this feasible, we model circular duplex DNA compactly by a B-spline ribbon-like model in terms of a small number of control vertices. We associate an elastic deformation energy composed of bending and twisting integrals and represent intrachain contact by a 6-12 Lennard Jones potential. The latter is parameterized to yield an energy minimum at the observed DNA-helix diameter inclusive of a hydration shell. A penalty term to ensure fixed contour length is also included. First and second partial derivatives of the energy function have been derived by using various mathematical simplifications. First derivatives are essential for Newton-type minimization as well as molecular dynamics, and partial second-derivative information can significantly accelerate minimization convergence through preconditioning. Here we apply a new large-scale truncated-Newton algorithm for minimization and a Langevin/implicit-Euler scheme for molecular dynamics. Our truncated-Newton method exploits the separability of potential energy functions into terms of differing complexity. It relies on a preconditioned conjugate gradient method that is efficient for large-scale problems to solve approximately for the search direction at every step. Our dynamics algorithm is numerically stable over large time steps. It also introduces a frequency-discriminating mechanism so that vibrational modes with frequencies greater than a chosen cutoff frequency are essentially frozen by the method. With these tools, we rapidly identify corresponding circular and interwound energy minima for small DNA rings for a series of imposed linking-number differences. These structures are consistent with available electron microscopy data. The energetic exchange of stability between the circle and the figure-8, in very good agreement with analytical results, is also detailed. Molecular dynamics trajectories at 100 femtosecond time steps then reveal the rapid folding of the unstable circular state into supercoiled forms. Significant bending and twisting motions of the interwound structures are also observed. Such information may be useful for understanding transition states along the folding pathway and the role of enzymes that regulate supercoiling.(ABSTRACT TRUNCATED AT 400 WORDS)

Algorithms