Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “computational modeling”

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 1,765 records · Page 98Linked to original sources

A new computational model for expiratory flow from nonhomogeneous human lungs.

A model has been developed for expiration from human lungs in which the mechanical properties of the airways and parenchyma can be varied between regions. The model is based on an existing homogeneous model. The fluid mechanical problem of the merging of dissimilar flows from adjacent regions is underspecified by the conservation laws of mass and energy. An existing, empirically derived result, provides the required extra equation. Model simulation of a nonhomogeneously distributed mild constriction of the peripheral airways gives results for maximal flows and alveolar pressure differences which are in good agreement with recent experimental findings.

Forced Expiratory Flow Rates↗

Comparisons of point and average organ dose within an anthropomorphic physical phantom and a computational model of the newborn patient.

Pediatric radiographic examinations yield medical benefits and/or diagnostic information that must be balanced against potential risk from patient radiation exposure. Consequently, clinical tools for measuring internal organ dose are needed for medical risk assessment. In this study, a physical phantom and Monte Carlo simulation model of the newborn patient were developed based upon their stylized mathematical expressions (ORNL and MIRD model series). The physical phantom was constructed using tissue equivalent substitutes for soft tissue, lung, and skeleton. Twenty metal-oxide-semiconductor field effect transistor (MOSFET) dosimeters were then inserted at three-dimensional positions representing the centroids of organs assigned in the ICRP's definition of the effective dose. MOSFET-derived point estimates of organ dose were shown to be in reasonable agreement with Monte Carlo estimates for representative newborn head, chest, and abdomen radiographic exams. Ratios of average organ dose assessed via MCNP simulations to the MOSFET-derived point doses (point-to-organ dose scaling factors, SF(POD)) are tabulated for subsequent use in clinical irradiations of the newborn phantom/MOSFET system. Values of SF(POD) indicate that MOSFET measurements of point dose for in-field exposures need to be adjusted only to within 10% to report volume-averaged organ dose. Larger adjustments to point doses are noted for organs out-of-field. For walled organs, point estimates of organ dose at the content centroid are shown to underestimate the average wall dose when the organ is within the primary field: SF(POD) of 1.19 for the stomach (AP chest exam), and SF(POD) of 1.15 for the urinary bladder (AP abdomen exam).

Abdomen↗

A computational model for electron backscattering in electron dosimetry.

A one-dimensional electron transport algorithm based on the ideas of invariant imbedding has been developed and applied to the problem of electron backscattering in electron dosimetry. Results have been compared with experimental and Monte Carlo results. Agreement appears to be reasonable at low incident energies (0.5 to 10 MeV) and low Z(Z < or = 65) materials.

Electrons↗

Regional septal dysfunction in a three-dimensional computational model of focal myofiber disarray.

MLC2v/ras transgenic mice display a phenotype characteristic of hypertrophic cardiomyopathy, with septal hypertrophy and focal myocyte disarray. Experimental measurements of septal wall mechanics in ras transgenic mice have previously shown that regions of myocyte disarray have reduced principal systolic shortening, torsional systolic shear, and sarcomere length. To investigate the mechanisms of this regional dysfunction, a three-dimensional prolate spheroidal finite-element model was used to simulate filling and ejection in the hypertrophied mouse left ventricle with septal disarray. Focally disarrayed septal myocardium was modeled by randomly distributed three-dimensional regions of altered material properties based on measured statistical distributions of muscle fiber angular dispersion. Material properties in disarrayed regions were modeled by decreased systolic anisotropy derived from increased fiber angle dispersion and decreased systolic tension development associated with reduced sarcomere lengths. Compared with measurements in ras transgenic mice, the model showed similar heterogeneity of septal systolic strain with the largest reductions in principal shortening and torsional shear in regions of greatest disarray. Average systolic principal shortening on the right ventricular septal surface of the model was -0.114 for normal regions and -0.065 for disarrayed regions; for torsional shear, these values were 0.047 and 0.019, respectively. These model results suggest that regional dysfunction in ras transgenic mice may be explained in part by the observed structural defects, including myofiber dispersion and reduced sarcomere length, which contributed about equally to predicted dysfunction in the disarrayed myocardium.

Animals↗

A computational model for expiratory flow.

A mathematical model of maximal expiratory flow was developed. Coupled equations describing the pressure losses in the flow and the pressure-area behavior of the airway were integrated along the airway from the periphery to the flow-limiting site. Equations describing pressure losses in the flow were adapted from studied of bronchial casts. The bronchial anatomy utilized was that described by Weibel. Bronchial mechanical properties were obtained from measurements in excised human lungs for the central airways and by extrapolations of these data for the peripheral airways. The maximal flow for air and helium predicted by the model agrees with that of five lungs from which mechanical properties were obtained. The model predictions agree with published values of density and viscosity dependence of maximal flow. At high and midlung volumes, maximal flow is determined primarily by the wave-speed mechanism. At low lung volumes, maximal flow is primarily determined by the coupling of viscous pressure losses and airway mechanical properties.

Biomechanical Phenomena↗

Mechanical compression-induced pressure sores in rat hindlimb: muscle stiffness, histology, and computational models.

Pressure sores affecting muscles are severe injuries associated with ischemia, impaired metabolic activity, excessive tissue deformation, and insufficient lymph drainage caused by prolonged and intensive mechanical loads. We hypothesize that mechanical properties of muscle tissue change as a result of exposure to prolonged and intensive loads. Such changes may affect the distribution of stresses in soft tissues under bony prominences and potentially expose additional uninjured regions of muscle tissue to intensified stresses. In this study, we characterized changes in tangent elastic moduli and strain energy densities of rat gracilis muscles exposed to pressure in vivo (11.5, 35, or 70 kPa for 2, 4, or 6 h) and incorporated the abnormal properties that were measured in finite element models of the head, shoulders, pelvis, and heels of a recumbent patient. Using in vitro uniaxial tension testing, we found that tangent elastic moduli of muscles exposed to 35 and 70 kPa were 1.6-fold those of controls (P < 0.05, for strains </=5%) and strain energy densities were 1.4-fold those of controls (P < 0.05, for strains >/=5%). Histological (phosphotungstic acid hematoxylin) evaluation showed that this stiffening accompanied extensive necrotic damage. Incorporating these effects into the finite element models, we were able to show that the increased muscle stiffness in widening regions results in elevated tissue stresses that exacerbate the potential for tissue necrosis. Interfacial pressures could not predict deep muscle (e.g., longissimus or gluteus) stresses and injuring conditions. We conclude that information on internal muscle stresses is required to establish new criteria for pressure sore prevention.

Animals↗

A computational model of muscle recruitment for wrist movements.

To execute a movement, the CNS must appropriately select and activate the set of muscles that will produce the desired movement. This problem is particularly difficult because a variety of muscle subsets can usually be used to produce the same joint motion. The motor system is therefore faced with a motor redundancy problem that must be resolved to produce the movement. In this paper, we present a model of muscle recruitment in the wrist step-tracking task. Muscle activation levels for five muscles are selected so as to satisfy task constraints (moving to the designated target) while also minimizing a measure of the total effort in producing the movement. Imposing these constraints yields muscle activation patterns qualitatively similar to those observed experimentally. In particular, the model reproduces the observed cosine-like recruitment of muscles as a function of movement direction and also appropriately predicts that certain muscles will be recruited most strongly in movement directions that differ significantly from their direction of action. These results suggest that the observed recruitment behavior may not be an explicit strategy employed by the nervous system, but instead may result from a process of movement optimization.

Humans↗

Long-range connections synchronize rather than spread intrathalamic oscillations: computational modeling and in vitro electrophysiology.

A thalamic network model was developed based on recent data regarding heterogeneous thalamic reticular (RE) cell axonal arborizations that indicate at least two projection patterns, short-range cluster projections and long-range diffuse projections. The model was constrained based on expected convergence and the biophysical properties of RE and thalamocortical (TC) cells and their synapses. The model reproduced in vitro synchronous slow (3-Hz) oscillatory activity and the known effects of T-channel blockade and cholecystokinin (CCK) application on this activity. Whereas previous models used the speed at which approximately 3-Hz oscillations propagate in vitro to infer the spatial extent of intrathalamic projections, we found that, so long as the gamma-aminobutyric acid-B synaptic conductance was adjusted appropriately, a network with only short-range projections and another network with both short- and long-range projections could both produce physiologically realistic propagation speeds. Although the approximately 3-Hz oscillations propagated at similar speeds in both networks, phase differences between oscillatory activity at different locations in the network were much smaller in the network containing both short- and long-range projections. We measured phase differences in vitro and found that they were similar to those that arise in the network containing both short- and long-range projections but are inconsistent with the much larger phase differences that occur in the network containing only short-range projections. These results suggest that, although they extend much further than do short-range cluster projections, long-range diffuse projections do not spread activity over greater distances or increase the speed at which intrathalamic oscillations propagate. Instead, diffuse projections may function to synchronize activity and minimize phase shifts across thalamic networks. One prediction of this hypothesis is that, immediately after a collision between propagating oscillations, phase gradients should vary smoothly across the thalamic slice. The model also predicts that phase shifts between oscillatory activity at different points along a thalamic slice should be unaffected by T-channel blockers and decreased by suppression of synaptic transmission or application of CCK.

Animals↗

Kinetics of triiodothyronine dissociation from bovine serum albumin: modification of the resin capture method with subsequent computer modeling.

There has been much debate about the mechanism by which thyroid hormone leaves the circulation and enters cells. Critical to the development of theoretical models of thyroid hormone uptake are precise estimates of the kinetic constants describing the reversible binding of the hormone to its binding proteins. A resin capture method was modified to permit the collection of dissociation data at intervals of 0.12 sec. Columns (30 X 0.6 cm) of Amberlite IRA 400 anion exchange resin were perfused at 30 ml/min with 3.5% dextran in modified Krebs Henseleit bicarbonate buffer at pH 7.4 and 37 C. Injection of [125I]triiodothyronine (T3) tracer in BSA-free medium established the rate constant for resin capture of free T3 in this system: kRes = 0.89 +/- 0.04 sec-1 (mean +/- SD, n = 4). [125I]T3 preequilibrated with BSA was infused to achieve final BSA concentrations of 2.3, 48, 120, 290, 507, and 725 microM, yielding apparent uptake rate constants (kapp) of 0.77 +/- 0.12 (n = 6), 0.23 +/- 0.01 (n = 4), 0.17 +/- 0.04 (n = 4), 0.057 +/- 0.005 (n = 4), 0.039 +/- 0.003 (n = 4), 0.020 +/- 0.002 (n = 4) sec-1 (SSE less than 9.4 X 10(-2), respectively, when fitted by nonlinear methods. Solving the equation kapp = kRes/(1 + [BSA]/KD') gave KD' = 1.80 +/- 0.07 X 10(-5) M (r2 = 0.96, P less than 0.001, n = 26) where KD' is the apparent or flow equilibrium constant. In a further experiment T3:BSA reassociation was eliminated by BSA dilution. This resulted in determination of the dissociation rate constant kdissoc = 0.63 +/- 0.16 sec-1 (n = 8) and therefore, the association rate constant kassoc = 3.5 X 10(4) M-1 sec-1. Mathematical modeling predicted that 14.4 microM BSA would decrease the rate of T3 uptake by the resin system by 65%. The free fraction of T3 during a single pass extraction was predicted to be greater than that measured by equilibrium dialysis, suggesting that dissociation of T3 from BSA is favored in the flux model relative to that measured at equilibrium in vitro.

Animals↗

Static forces and moments generated in the insect leg: comparison of a three-dimensional musculo-skeletal computer model with experimental measurements

As a first step towards the integration of information on neural control, biomechanics and isolated muscle function, we constructed a three-dimensional musculo-skeletal model of the hind leg of the death-head cockroach Blaberus discoidalis. We tested the model by measuring the maximum force generated in vivo by the hind leg of the cockroach, the coxa&shy;femur joint angle and the position of this leg during a behavior, wedging, that was likely to require maximum torque or moment production. The product of the maximum force of the leg and its moment arm yielded a measured coxa&shy;femur joint moment for wedging behavior. The maximum musculo-apodeme moment predicted by summing all extensor muscle moments in the model was adequate to explain the magnitude of the coxa&shy;femur joint moment produced in vivo by the cockroach and occurred at the same joint angle measured during wedging. Active isometric muscle forces predicted from our model varied by 3.5-fold among muscles and by as much as 70 % with joint angle. Sums of active and passive forces varied by less than 3.5 % over the entire range of possible joint angles (0&shy;125 &deg;). Maximum musculo-apodeme moment arms varied nearly twofold among muscles. Moment arm lengths decreased to zero and switched to the opposite side of the center of rotation at joint angles within the normal range of motion. At large joint angles (>100 &deg;), extensors acted as flexors. The effective mechanical advantage (musculo-apodeme moment arm/leg moment arm = 0.10) resulted in the six femoral extensor muscles of the model developing a summed force (1.4 N) equal to over 50 times the body weight. The model's three major force-producing extensor muscles attained 95 % of their maximum force, moment arm and moment at the joint angle used by the animal during wedging.

Journal Article↗

Wave-induced forces on the giant kelp Macrocystis pyrifera (Agardh): field test of a computational model

We propose a hydro-mechanical numerical model that predicts the maximal tension to which stipes of the giant kelp Macrocystis pyrifera will be subjected when exposed to ocean waves. Predicted maximal tensions are close to those measured in the field. The strength of Macrocystis pyrifera stipes was measured, allowing our prediction of forces to be translated into a prediction of the fraction of stipes broken. Predicted breakage is low even for extreme storm waves, a testament to the mechanical design of individual kelp fronds. However, empirically measured rates of kelp mortality can be high, considerably higher than those predicted on the basis of hydrodynamic forces acting alone. This indicates that factors not taken into account in our model (such as holdfast dislodgment, entanglement of stipes, damage from herbivory and wave breaking) contribute substantially to mortality in Macrocystis pyrifera.

Journal Article↗

Computational models to predict aqueous drug solubility, permeability and intestinal absorption.

In the last decade, poor intestinal absorption of candidate drugs intended for oral administration has been identified as a major bottleneck in drug development. Poor intestinal absorption can often be related to poor aqueous solubility and/or poor permeability across the intestinal wall. Other factors, such as poor stability and the metabolism of the compounds, can also decrease the amount of compound absorbed. In an effort to design compounds with enhanced absorption profile, theoretical predictions of solubility and permeability, among other factors, have gained increased interest, and a large number of papers have been published. In this review, the databases and techniques used for the development of in silico absorption models will be discussed. The focus is on aqueous drug solubility, which has become a major problem in drug development.

Intestinal Absorption↗

Development of computational models for the purpose of conducting individual livestock and premises traceback investigations utilizing National Animal Identification System-compliant data.

Many of the efforts surrounding the development of the National Animal Identification System have encompassed the identification of livestock production and handling premises as well as individuals or herds of animals, whereas little effort has been directed toward the ultimate goal of animal traceback within 48 h. A mock data set representative of the Colorado cattle population was created for modeling of cattle traceability. Using this data set, algorithms were developed to complete rapid and accurate traceback and traceforward of animals or premises or both. On July 19, 2005, the Colorado Department of Public Health and Environment, in conjunction with the Colorado Department of Agriculture, conducted a test exercise pertaining to homeland security. The exercise team randomly identified animal number 926,583 (of the 2 million total animals) as a potentially infected animal of interest and requested a traceback of this animal. Traceback was accomplished in 215 s, and 540 primary coresident animals were identified. However, due to animal movements, the number of coresidents (animals exposed, directly or indirectly, to the animal of interest) expanded with coresidency level (level 1 = direct contact; level 2 = direct contact with an animal that had direct contact with the animal of interest; level 3 = direct contact with an animal that had contact with an animal that had direct contact with the animal of interest, etc.) to more than 1.2 million coresidents at level 4, and more than 90% of all animals identified as a coresident at some level. In addition to the coresidency results, the premises containing the coresidents were identified and sorted by the number of coresidents. Because of animal movement, all 19,391 premises included in the data set had coresidents at some level. This exercise demonstrated the capability of the developed algorithms to complete rapid traceback and the complexity of the resulting animal traceback output.

Agriculture↗

Evaluating a computational model of perceptual grouping by proximity.

A formal approach to the phenomenon of perceptual grouping by proximity was investigated. Grouping judgments of random dot patterns were made by the CODE algorithm (van Oeffelen & Vos, 1982) and several related algorithms, and these judgments were compared with subjects' grouping judgments for the same stimuli. Each algorithm predicted significantly more subject judgments than would be expected by chance. The more subjects agreed on how a given dot pattern should be grouped, the more successful was the algorithms' ability to match the judgments for that pattern. CODE predicted significantly fewer subject judgments than did some of the other algorithms, largely because of its overemphasis on the extent of interactivity among dots as they are being grouped.

Adult↗