PubMed Health⌕ Search

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

Fluorescence quenching of dyes by tryptophan: interactions at atomic detail from combination of experiment and computer simulation.

Fluorescence spectroscopy and molecular dynamics (MD) simulation are combined to characterize the interaction of two organic fluorescent dyes, rhodamine 6G (R6G) and an oxazine derivative (MR121), with the amino acid tryptophan in aqueous solution. Steady-state and time-resolved fluorescence quenching experiments reveal the formation of essentially nonfluorescent ground-state dye/Trp complexes. The MD simulations are used to elucidate the molecular interaction geometries involved. The MD-derived probability distribution of the distance r between the centers of geometry of the dye and quencher ring systems, P(r), extends to higher distances for R6G than for MR121 due to population in the R6G/Trp system of fluorescent interaction geometries between Trp and the phenyl ring and ester group of the dye. The consequence of this is the experimental finding that under the conditions used in the simulations about 25% of the R6G dye is fluorescent in comparison with 10% of the MR121. Combining the above findings allows determination of the "quenching distance", r, above which no quenching occurs. r is found to be very similar (approximately 5.5 A) for both dye/Trp systems, corresponding to close to van der Waals contact. Both experimental dynamic Stern-Volmer analysis and the MD trajectories demonstrate that the main determinant of the fluorescence intensity is static quenching. The approach presented is likely to be useful in the structural interpretation of data obtained from fluorescent conjugates commonly used for monitoring the binding and dynamics of biomolecular systems.

Computer Simulation↗

Spatial distribution of polyelectrolyte and counterions in nanocapsules: a computer simulation study.

We study a spatial distribution of polyelectrolyte chains and counterions inside nanometer-size capsules by means of molecular dynamics simulation on the level of a colloidal model in which polyelectrolyte coils are modeled as soft charged spheres. The capsule shell is treated as a semipermeable membrane, impermeable for the polyelectrolyte chains, but allowing free diffusion of solvent molecules and counterions. As a result, counterions leak out from the capsule immersed into a fluid of low ionic strength. This counterion leakage leads to a formation of characteristic polyelectrolyte density profiles with the central plateau and large peaks at the wall. We show that a nonuniform distribution of the inner polyelectrolyte depends on the capsule radius, surface charge, concentration of encapsulated polyelectrolyte, and the volume fraction of capsules.

Adsorption↗

Computer simulation of the accident with nine victims.

In the paper we wish to emphasise the significance of vehicle driving dynamics analysis in the collision phase and occupant load analysis by means of using a software environment. Thereby we also wish to present the results of the simulation of the course of a traffic accident with nine victims that arose from a collision between an Audi A6 passenger car and the VW Caravelle van. In treating the traffic accident the forensic expert was faced with the questions about what caused the injuries to the front passenger in the Audi A6 passenger car, about the way the two vehicles had collided, about their collision velocities, about the way the two vehicles were handled and about the causes that originated the traffic accident. The critical situation on the road was a consequence of the tiredness of the van driver, the inadequate use of the passive safety systems and overloading the van.

Acceleration↗

Computer simulation of final height after growth-promoting therapy.

This short communication discusses the possible influence of final height prediction on the estimation of the growth response to therapy. These influences are illustrated in terms of height at 8 and 18 years of age for 28 prepubertally short boys. On the basis of several treatment simulations and assumptions, we conclude that the mean residual final height; that is, the mean of individual differences between final height and predicted final height, is a robust outcome measure in growth-promoting studies. However, the treatment residual final height (the difference between final height and predicted final height following growth-promoting treatment) was only partially correlated (r2 = 0.27) with the gain due to treatment, but was highly correlated (r2 = 0.75) with the residual final height; thus, growth-response estimations made on the basis of treatment residual final height may be more dependent on the predicted interval, or confidence interval, of final height than on the effect of growth-promoting treatment.

Adolescent↗

A parsimonious description of motoneuron dendritic morphology using computer simulation.

Most quantitative descriptions of neuronal dendrite morphology involve tabulations of measurements and correlations among them. The present work is an attempt to extract from such data a parsimonious set of parameters that are sufficient to describe the quantitative features of individual and pooled dendrites, including their statistical variability. A relatively simple stochastic (Monte Carlo) model was devised to simulate branching dendritic trees. The necessary parameters were then derived directly from measurements of 64 completely reconstructed dendrites belonging to six gastrocnemius alpha-motoneurons, labeled by intracellular injection of HRP. Comparison of actual and simulated dendrites was used to guide the process of parameter extraction. The model included only two processes, one to generate individual branches given their starting diameters and the second to select starting diameters for the daughter branches produced at dichotomous branching points. The stochastic process for branch generation was controlled by probability functions for branching (Pbr) and for terminating (Ptrm), together with a constant rate of branch taper. All model parameters were fixed by motoneuron measurements except for branch taper rate, which was allowed to vary within limits consistent with observed taper rates in order to generate the appropriate total number of branches. The simplest model (model 1), in which Pbr and Ptrm depended only on local branch diameter, produced simulated dendrites that fit many, but not all, characteristics of actual motoneuron dendrites. Two additional properties produced significant improvements in the fit: (1) a small but significant dependence of daughter diameters on the normalized starting diameter of the parent branch, and (2) a dependence of Pbr and Ptrm on distance from the soma as well as on local branch diameter. The process of developing this model revealed unsuspected relations in the original data that suggest the existence of fundamental mechanisms for morphological control. The final model succinctly describes a large amount of data and will enable quantitative comparisons between the dendritic structures of different types of neurons, regardless of their relative sizes.

Animals↗

Computer simulation of re-entry sources in myocardium in two and three dimensions.

Re-entry in a simple model of cardiac muscle is examined. The properties of vortices in two dimensions are studied at different values of the ratio of the action potential duration to its rise time. The vortex changes from stationary to non-stationary rotation as the action potential duration is increased. A test for identifying non-stationary rotation from monophasic action potential recordings at a single site in cardiac tissue is proposed. In three dimensions re-entry of the scroll ring type is examined. Contraction of the scroll ring and its vertical drift occurs in this model. The drift velocities are shown to depend inversely on the scroll ring radius. The dependence of these drift velocities on re-entry cycle length is estimated for cardiac tissue.

Action Potentials↗

A visual study of surface potentials and Laplacians due to distributed neocortical sources: computer simulations and evoked potentials.

A "picture book" of surface potentials, Laplacians, and magnetic fields due to distributed, neocortical sources is presented. The mathematically simulated data is based on 4200 current sources at the macrocolumn scale. Estimated scalp surface maps are based on the three-concentic spheres model of the head. Emphasis is placed on the effects of sampling with a limited number of electrodes, the choice of reference electrode, and the use of the spline Laplacian to improve spatial resolution. The spline Laplacian is applied to median and ulnar nerve somatosensory evoked potentials and to auditory evoked potentials including P300. Substantial improvement in spatial resolution over conventional methods is obtained. The implementation of practical high resolution EEG systems based on the spline Laplacian is considered.

Algorithms↗

Prediction of infusion rates: validation of a computer simulation using vecuronium.

In Phase 1 of the study a loading dose of vecuronium was given to 24 patients, followed by a sequence of infusions and boluses. The relationship between the duration of a bolus dose and the infusion rate necessary to maintain the required degree of neuromuscular block was demonstrated. In Phase 2 the data acquired from Phase 1 were used in an attempt to produce a predetermined degree of neuromuscular block. The aim was a block of neuromuscular function so that the height of the first twitch of a train-of-four was 15% of control. In the 10 patients studied the average height was 14.7% (SD 4.95%). It was concluded that it is possible to predict accurately the infusion dose of vecuronium required to achieve a predetermined degree of neuromuscular block, according to the duration of action of one or two test doses.

Adolescent↗

Direct computer simulation of water-mediated force between supported phospholipid membranes.

The grand canonical Monte Carlo technique is used to calculate the water-mediated force operating between two supported 1,2-dilauroyl-DL-phosphatidylethanolamine (DLPE) membranes in the short separation range. The intra- and intermolecular interactions in the system are described with a combination of an AMBER-based force field for DLPE and a TIP4P model for water. The long range contributions to the electrostatic interaction energy are treated in the dipole-dipole group-based approximation. The total water-mediated force is analyzed in terms of its hydration component and the component due to the direct interaction between the membranes. The latter is, in addition, partitioned into the electrostatic, van der Waals, and steric repulsion contributions to give an idea of their relative significance in the water-mediated interaction of the membranes.

Computer Simulation↗

Active site dynamics of the HhaI methyltransferase: insights from computer simulation.

A molecular dynamics study was performed on the DNA methyltransferase M. Hha I in a ternary complex with DNA and AdoMet in solution. Methylation involves addition of the Cys81 sulfhydryl anion to the 6-position of Cyt18, followed by a nucleophilic attack of the resultant carbanion at C5 on the AdoMet methyl group. It was found in this simulation that the distances between the sulfhydryl group (SG) of Cys81 to the C6 of Cyt18 (SG-C6) and methyl carbon (CH3) of AdoMet to the C5 of cytosine (CH3-C5) are dependent on the dihedral angle chi (O4'-C1'-N1-C2) of the nucleotide. When the chi angle of Cyt18 is low (< -80 degrees), the SG-C6 and CH3-C5 distances are large. A high chi angle (> -80 degrees) for the target cytosine residue reduces the distances for both SG-C6 and CH3-C5, and the angles formed between the cytosine ring and AdoMet correspond well to values for the transition state structures formed during methylation of cytosine from ab initio calculations. Two possible proton sources for protonation of N3 of the cytosine residue upon formation of the covalent intermediate were found in the simulation. The protonated amine group of AdoMet could provide a proton via a water bridge, or Arg163 could also be the source of the proton for N3 via a water bridge. The simulation provides insights into how the H5 of cytosine could go from the active site into solvent. Conserved residues Asn304 and Gln82 stabilize a water network within the active site of M. Hha I which provides a route for H5 to diffuse into bulk solvent. An initially distant water molecule was able to diffuse into the active site of the enzyme and replace a position of a crystallographic water molecule in close proximity to the C5 of cytosine. The movement of this water molecule showed that a channel exists between Gln82 and the AdoMet in M. Hha I which allows both water and protons to easily gain access to the active site of the enzyme.

Bacterial Proteins↗

A computer simulation of functional group contributions to free energy in water and a DPPC lipid bilayer.

A series of all-atom molecular dynamics simulations has been performed to evaluate the contributions of various functional groups to the free energy of solvation in water and a dipalmitoylphospatidylcholine lipid bilayer membrane and to the free energies of solute transfer (Delta(DeltaG(o))X) from water into the ordered-chain interior of the bilayer. Free energies for mutations of the alpha-H atom in p-toluic acid to six different substituents (-CH3, -Cl, -OCH3, -CN, -OH, -COOH) were calculated by a combined thermodynamic integration and perturbation method and compared to literature results from vapor pressure measurements, partition coefficients, and membrane transport experiments. Convergence of the calculated free energies was indicated by substantial declines in standard deviations for the calculated free energies with increased simulation length, by the independence of the ensemble-averaged Boltzmann factors to simulation length, and the weak dependence of hysteresis effects on simulation length over two different simulation lengths and starting from different initial configurations. Calculated values of Delta(DeltaG(o))X correlate linearly with corresponding values obtained from lipid bilayer transport experiments with a slope of 1.1 and from measurements of partition coefficients between water and hexadecane or decadiene, with slopes of 1.1 and 0.9, respectively. Van der Waals interactions between the functional group of interest and the acyl chains in the ordered chain region account for more than 95% of the overall potential energy of interaction. These results support the view that the ordered chain region within the bilayer interior is the barrier domain for transport and that solvation interactions within this region resemble those occurring in a nonpolar hydrocarbon.

1,2-Dipalmitoylphosphatidylcholine↗

Molecular basis for microbial adhesion to geochemical surfaces: computer simulation of Pseudomonas aeruginosa adhesion to goethite.

The adhesion of Pseudomonas aeruginosa to the goethite mineral is investigated using classical molecular simulation. A fragment model for goethite has been integrated into a fully atomistic membrane model. Properties for the resulting system are evaluated for a 1.5-ns simulation in the isothermal-isobaric ensemble. The response of the membrane to the presence of the mineral is investigated. Radial distribution functions are used to present an average picture of the hydrogen bonding. Orientational vectors, assigned to the saccharide groups, reveal the extent of the mineral's perturbations on the membrane. Significant structural changes were observed for the outermost saccharide groups, several of which rotate to form hydrogen bonds with the mineral surface. The structure of the inner core, and the corresponding integrity of the membrane, is maintained. The mineral surface dehydrates slightly in the presence of the membrane as saccharide hydroxyl groups compete with water molecules for hydrogen-bonding sites on its surface.

Bacterial Adhesion↗

Computer simulation for the prediction of separation as a function of pH for reversed-phase high-performance liquid chromatography. II. Resolution as a function of simultaneous change in pH and solvent strength.

The optimization of reversed-phase high-performance liquid chromatographic separation by the simultaneous variation of pH and solvent strength (%B) was studied for acidic (substituted benzoic acids) and basic samples (substituted anilines). The combination of these two variables was expected to be more useful than either variable alone. This proved to be the case for the benzoic acid sample, but not for the aniline sample. Column plate numbers were also studied for each sample and as a function of pH. With the exception of one compound (3,5-dimethylaniline) in one particular pH range (3.0-4.5), plate numbers of 12,000-20,000 were observed for each sample.

Aniline Compounds↗

A computer simulation of the oxygen balance in a cold climate winter storage WSP during the critical spring warm-up period.

The paper considers factors that determine the oxygen balance in extreme climate waste stabilisation ponds during the critical spring warm-up period. At this time BOD load on the pond is a maximum, due to accumulation of wastewater under the ice during the winter. The paper describes the operation of a typical cold climate WSP and the events leading to a balanced steady state system as spring develops into summer. A mathematical model to simulate conditions within a batch fed experimental pond over the transient period is described. To model temperature changes in the water body experimental data were fitted to a generalised equation based on diurnal fluctuations in air temperature. The results are plotted in a normalised form and show the diurnal fluctuation and time lapse as the depth of the pond increases. Maximum daily water temperature lags behind maximum light intensity. Bacterial growth is simulated by a Monod kinetic model in which growth rate depends on initial substrate concentration; temperature compensation is applied using a temperature activity coefficient. Oxygen utilisation is calculated from substrate removal. Algal growth rate is more complicated as it is affected by temperature and light availability. Algal oxygen production potential is considered in terms of its primary metabolite yield, which is then used in a Monod equation to estimate the growth rate. The model uses a mass balance approach to determine dissolved oxygen concentration in the pond. The model is still in a simple form but shows reasonable agreement, in terms of events and time lapses, to measured parameters in experimental ponds recovering from ice cover.

Bacteria↗

A four-dimensional computer simulation model of the in vivo response to radiotherapy of glioblastoma multiforme: studies on the effect of clonogenic cell density.

Tumours behave as complex, self-organizing, opportunistic dynamic systems. In an attempt to better understand and describe the highly complicated tumour behaviour, a novel four-dimensional simulation model of in vivo tumour growth and response to radiotherapy has been developed. This paper presents the latest improvements to the model as well as a parametric validation of it. Improvements include an advanced algorithm leading to conformal tumour shrinkage, a quantitative consideration of the influence of oxygenation on radiosensitivity and a more realistic, imaging based description of the neovasculature distribution. The tumours selected for the validation of the model are a wild type and a mutated p53 gene glioblastomas multiforme. According to the model predictions, a whole tumour with larger cell cycle duration tends to repopulate more slowly. A lower oxygen enhancement ratio value leads to a more radiosensitive whole tumour. Higher clonogenic cell density (CCD) produces a higher number of proliferating tumour cells and, therefore, a more difficult tumour to treat. Simulation predictions agree at least semi-quantitatively with clinical experience, and particularly with the outcome of the Radiation Therapy Oncology Group (RTOG) Study 83-02. It is stressed that the model allows a quantitative study of the interrelationship between the competing influences in a complex, dynamic tumour environment. Therefore, the model can already be useful as an educational tool with which to study, understand and demonstrate the role of various parameters in tumour growth and response to irradiation. A long term quantitative clinical adaptation and validation of the model aiming at its integration into the treatment planning procedure is in progress.

Algorithms↗

Comparison of computer simulated and phantom measured phase variance in the study of trabecular bone.

A new method using magnetic resonance phase images for the assessment of trabecular bone structure has recently been proposed. To test this method, a mathematical model is developed which calculates the phase distribution in gradient echo acquired phase images of a structure of Pyrex glass rods immersed in a copper sulfate solution. Several experiments were performed using a phantom built in the same way as the structure used in the mathematical model. The results from the model are compared with those from the phantom tests, and the influence of resolution and bone area fraction on the phase dispersion is studied. The good correlation between theoretical and experimental results shows that phase variance increases with increasing resolution and bone density. However, the dependence of variance on bone density is less prominent for large pixel sizes.

Bone Density↗

A computer simulation of the static magnetic field distribution in the human head.

Distortion of the static magnetic field inside the human head is dependent on regional tissue susceptibility variations and geometrical shape. These effects result in resonance line broadening and frequency shifts and consequently, intensity and spatial errors in both magnetic resonance imaging (MRI) and magnetic resonance (MR) spectroscopy. To calculate the field distortion due to the susceptibility's geometry, two dimensional (2D) finite element analysis was applied to simulate the field distribution in a 2D model of the human head, placed in a uniform magnetic field. The model contains air-filled cavities and sinuses, and the remainder is treated as water. The magnetic field deviation was evaluated using gray scale plots and histograms of the magnetic field. The shifts in parts/million and broadening of the histograms correspond to the NMR of the sampled region. The field distribution of the human head was also experimentally mapped using the DANTE tagging sequence. The calculated and experimental field maps are in good agreement. Thus, geometric considerations with uniform susceptibilities are sufficient to explain most of the static magnetic field distribution in the human head.

Air↗

Computer simulation of strand break yields in plasmid pBR322: DNA damage following 125I decay.

This paper presents results of (125)I effects on plasmid pBR322 in aqueous solution, simulating the complete transport of Auger and X rays up to the chemical phase. In addition to new sampling algorithms, new electronic cross sections are included. Simulations were carried out both with (125)I, bound to plasmid, or free, in its vicinity. The influence of the hydroxyl radical scavenger dimethyl sulfoxyde (DMSO) has also been tested, underlying that, in naked DNA, double strand breaks (caused by the decay of bound (125)I) are mainly due to direct hits. The calculated yields of relaxation events (RE) and linearization events (LE) show good agreement with experimental ones: when (125)I is bound to the plasmid pBR322, 0.16 RE and 0.83 LE per decay (without DMSO) are then observed. Then, when 2 mol DMSO is added, RE and LE probabilities become 0.22 and 0.76. The very light differences with those from literature could arise from experimental conditions.

Base Pair Mismatch↗