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 757 records · Page 42Linked to original sources

Computational modeling of a binding conformation of the intermediate L-histidinal to histidinol dehydrogenase.

Histidinol dehydrogenase (HDH) is one of the enzymes involved in the L-histidine biosynthesis pathway. HDH is a dimer that contains one Zn2+ ion in each identical subunit. In this study, we predicted a possible binding conformation of the intermediate L-histidinal, which is experimentally not known, using a computational modeling method and three potent HDH inhibitors whose structures are similar to that of L-histidinal. At first, a set of the most probable active conformations of the potent inhibitors was determined using two different pharmacophore mapping techniques, the active analogue approach and the distance comparison method. From the most probable active conformations of the three potent inhibitors, the common parts of the L-histidinal structure were extracted and refined by energy minimization to obtain the binding conformation of L-histidinal. This predicted conformation of L-histidinal agrees with an experimentally determined conformation of L-histidine in a single crystal, suggesting that it is an experimentally acceptable conformation. The capability in this conformation to coordinate a Zn2+ ion was examined by comparing the spatial relative geometry of its functional groups with those of ligands that coordinate with a Zn2+ ion in Zn proteins of the Protein Data Bank. This comparison supported our predicted conformation.

Alcohol Oxidoreductases↗

Electric pacing of intact and transected canine small intestine and its computer model.

In 10 conscious, fasted dogs with electrodes chronically implanted on the intestine, current pulses (8 mA, 50 ms) at frequencies the same as, or faster than, that of the natural intestinal pacemaker always entrained pacesetter potentials (PP) along the proximal frequency plateau, but not along the distal frequency gradient. As duodenal PP's were paced faster, the proximal plateau shortened by orad extension of the distal gradient. Entrained PP's propagated orally and aborally. Their velocity slowed caudally while varying inversely with their frequency. After midduodenal transection had reduced the frequency of the natural PP's distal to the cut, pacing entrained PP's in all areas of the bowel and restored the proximal plateau and distal gradient. But no area could be paced faster than before transection. PP's were not propagated across sites of transection. A digital computer model of coupled relaxation oscillators gave similar results. We conclude that the frequency of the natural pacemaker and the declining gradient of maximal driven frequency determine the frequency pattern of the small intestinal PP.

Animals↗

Post-lesion lateralisation shifts in a computational model of single-word reading.

The mechanisms underlying lateralisation of language are incompletely understood. Existing data is inconclusive, for example, in determining which underlying asymmetries in hemispheric anatomy/physiology lead to lateralisation, the precise role of interhemispheric connections in this process, and exactly how and why lateralisation can shift following focal brain damage. Although these issues will ultimately be settled by experimentation, it is likely that computational modelling can be used to suggest, focus, and even interpret such empirical work. We have recently studied the emergence of lateralisation in an artificial neural network model having paired cerebral hemispheric regions, as the model learned to generate the correct pronunciation for simple words. In this paper we extend this previous work by examining the immediate and longer-term changes in lateralisation that occur following simulated acute hemispheric lesions. Among other things, the results demonstrate that the extent to which the non-lesioned model hemispheric region contributes to recovery is a function of lesion size, prelesion lateralisation, and assumptions about the excitatory/inhibitory influences of the corpus callosum. The relevance of these results to the currently controversial suggestion that language lateralisation shifts following focal damage to language areas, and that the unlesioned hemisphere contributes to recovery from stroke-induced aphasia in adults, is discussed.

Journal Article↗

Computer modelling: future directions.

Recent developments in computing and in the theory of simulation have extended greatly the successes of the modelling of ionic crystals pioneered by Mott and Littleton. This has changed the way in which computer experiments are brought to bear on an increasing range of solid-state phenomena. Yet applied science creates new demands, both in the form of new types of system and in terms of the complexity and subtlety of what is studied. The author's brief survey looks at some of the successes and gaps from interfaces and catalysts to neurotransmitters and from superconductors to slags.

Computer Simulation↗

Reflection after delayed excitation in a computer model of a single fiber.

Reflection (reflected reentry) is a case of reentry in a one-dimensional structure, divided into proximal and distal segments, in which tissue excited by a wave front propagating in a forward direction is reexcited by electrical activity coming backward from the original direction of propagation. Cases of reflection have been demonstrated in Purkinje fibers and in ventricular muscle preparations containing multiple fibers. Several mechanisms possibly responsible for reflected reentry have been proposed. However, the difficulty in the interpretation of the experimental results, as well as the limited number of different conditions in which reflection was obtained, has kept open the question about conditions and mechanisms for reflection. We have developed a computer model in which reflection occurs. The model involves a single fiber and uses the DiFrancesco-Noble equations for the Purkinje fiber to model the ionic currents. The results show that reflection is possible in a single fiber and that diastolic depolarization (automaticity) is not a requirement for reflection. Active membrane responses to a just-above-threshold stimulus were important for achieving the necessary time delay. Systematic simulations showed further that reflection occurred only when the right coupling conditions linked a short or long proximal fiber to a short distal segment.

Action Potentials↗

Computer model of clonazepam's effect in thalamic slice.

In the thalamus, paradoxical changes in response to augmentation of inhibition can occur as a result of either cellular or network effects. Clonazepam, a GABA(A) agonist, produces a paradoxical reduction in evoked thalamocortical neuron inhibitory postsynaptic potential (IPSP) in thalamic slice. This has been hypothesized to be a result of augmentation in inhibitory to inhibitory connections. In a computer model, orthodromic simulation produced an increase in initial IPSP, a result contrary to that found experimentally. This failure was traced to the inability of orthodromic activation to produce fast enough recurrent inhibition to alter initial reticularis neuron firing. Simulated antidromic stimulation was able to reduce this initial spike train and reproduced the experimental finding.

Animals↗

Computational models of predictive and memory-related functions of the hippocampus.

We discuss the role of the hippocampus in information processing in the brain and hypothesise that the hippocampus monitors the stability of sensory cues it receives from the external world, using the current context to predict the next sensory event in the episodic sequence by learning from experience, and memorising these sequences of sensory events. Two computational models are presented here. The predictive theory and model are closely related to experimental evidence and use dynamic synapses with an asymmetric learning rule to develop predictive neural activity of a leaky integrate-and-fire model of a pyramidal CA3 cell. The oscillatory model of the hippocampus for memorising sequences of sensory events is developed as a chain of interacting neural oscillators forced by oscillatory inputs from the entorhinal cortex and from the medial septum.

Animals↗

Computational models for predicting interactions with cytochrome p450 enzyme.

Cytochrome p450 (CYP) enzymes are predominantly involved in Phase 1 metabolism of xenobiotics. As only 6 isoenzymes are responsible for approximately 90 % of known oxidative drug metabolism, a number of frequently prescribed drugs share the CYP-mediated metabolic pathways. Competing for a single enzyme by the co-administered therapeutic agents can substantially alter the plasma concentration and clearance of the agents. Furthermore, many drugs are known to inhibit certain p450 enzymes which they are not substrates for. Because some drug-drug interactions could cause serious adverse events leading to a costly failure of drug development, early detection of potential drug-drug interactions is highly desirable. The ultimate goal is to be able to predict the CYP specificity and the interactions for a novel compound from its chemical structure. Current computational modeling approaches, such as two-dimensional and three-dimensional quantitative structure-activity relationship (QSAR), pharmacophore mapping and machine learning methods have resulted in statistically valid predictions. Homology models have been often combined with 3D-QSAR models to impose additional steric restrictions and/or to identify the interaction site on the proteins. This article summarizes the available models, methods, and key findings for CYP1A2, 2A6, 2C9, 2D6 and 3A4 isoenzymes.

Computational Biology↗

A comparison of radiation dose measured in CT dosimetry phantoms with calculations using EGS4 and voxel-based computational models.

CT is a high-dose examination and possibly the dominant contributor to dose from diagnostic radiology. Estimates of organ doses are obtained from Monte Carlo calculations and used to quantify radiation risk. To ensure the validity of using Monte Carlo calculations to estimate actual dose, measurements must be compared with calculations. We have measured doses to CT head and chest dosimetry phantoms and compared them with Monte Carlo (EGS4) calculated doses in voxel-based computational models of the phantoms. The simulation used an x-ray spectrum calculated from the specified values of the scanner's x-ray tube parameters. The scanner's beam-shaping filter was included in the modelling. Measured and calculated doses to both the head and chest phantoms agreed to within 7%. The inclusion of Rayleigh scattering in the calculations has a significant effect if only one slice is scanned but not if multiple slices are scanned.

Head↗

Computer modeling of laser damage to the eye.

The effects of laser damage in the eye have been studied with an electrophysiological measure, the Visual Evoked Potentials (VEP). The VEP is generated by placing surface electrodes at the back of a subject's head, and having the subject watch a pattern of alternating light and dark bars. The spatial frequency of the bar pattern, measured in cycles per degree, increases as the bars become finer. The bar pattern changes in time, as the light bars become dark, and visa versa. The rate of this change is the temporal frequency, measured in Hertz. The VEP amplitudes referred to in this paper are a sum of the first three even harmonics of the temporal frequency (since the pattern change two times for every temporal cycle). As spatial frequency increases, a point will be reached where the subject will no longer be able to resolve the bars, and the VEP amplitude will become indistinguishable from noise. This spatial frequency is a measure of visual acuity. A computer model has been developed to simulate vision loss from lesions on the retina. It has been checked against experimental data from Rhesus monkeys with burns made under laboratory conditions. Vision was checked before and after the burns with VEP at a temporal frequency of six Hertz. This model divides the visual field into concentric rings centered about the visual axis. Each ring contributes to the VEP by an amount that depends on its eccentricity, or radius in degrees, as well as the spatial frequency.(ABSTRACT TRUNCATED AT 250 WORDS)

Animals↗

Computer model for study of cardiac repolarization.

INTRODUCTION: We propose a new and simple method to model repolarization in the left ventricle and the corresponding T waves on the surface ECG. METHODS AND RESULTS: We modeled the cardiac cell action potentials (APs) in the left ventricle (LV) with differences in only the duration of the plateau phase. Using published experimental data on the epicardial and endocardial repolarization sequences, for each point on the left ventricular surface we set a different AP repolarization starting time, determined by the duration of the plateau phase. The surface source model was used to compute potentials on the surface of the torso, generated by repolarization of the LV. Both the torso and the LV had homogeneous and isotropic conductivity. We simulated T waves on the 12-lead ECG and compared our results with measured T waves from five normal subjects. The orientation and shape in each lead were reproduced. In each lead we computed the root mean square error between simulated and measured T waves. The average error across the 12 leads was small, with a mean value of 0.11 mV across all the subjects. CONCLUSION: Repolarization of the LV can be modeled independently of the depolarization sequence and AP duration gradients. This method is an easy and powerful tool to describe the ECG features of repolarization.

Action Potentials↗

Kinetic analysis of the human erythrocyte glyoxalase system using 1H NMR and a computer model.

1H NMR was used with methylglyoxal, purified by an HPLC technique, to study the kinetics of the human erythrocyte glyoxalase system. 1H NMR enabled the direct measurement of the time-dependent changes in concentrations of the two hydrates of methylglyoxal, which have not previously been directly measurable, as well as measurement of substrates and products of the glyoxalase enzyme system in the human red blood cell. A computer model of the reaction scheme was developed and NMR data numerically analyzed, thus allowing a complete kinetic description of the reactions. The rate constants describing the chemical equilibria between the hydrated species of methylglyoxal were determined by this numerical analysis or by a saturation-transfer technique, and found to be much slower (by several orders of magnitude) than previously determined by other methods. The kinetic parameters describing the enzyme-catalyzed reactions were also determined from experiments using a dilute haemolysate that was added to solutions of methylglyoxal and reduced glutathione (GSH). The maximal velocity of glyoxalase 1 is threefold greater (Vmax = 70.4 +/- 4.7 mmol.min-1.1 packed cells-1) than glyoxalase 2(Vmax = 24 +/- 5 mmol.min-1.1 packed cells-1) and it exhibits threefold-greater affinity for its substrate (Km = 0.46 +/- 0.04 mM) than the second enzyme (Km = 1.5 +/- 0.4 mM). Both enzymes are subject to competitive inhibition; glyoxalase 1 by reduced glutathione (KiGSH = 7.88 +/- 0.16 mM) and glyoxalase 2 by the hemithioacetal (HTA) of methylglyoxal and GSH (KiHTA = 0.29 +/- 0.04 mM).

Erythrocytes↗

Distinct modes of collagen type I proteolysis by matrix metalloproteinase (MMP) 2 and membrane type I MMP during the migration of a tip endothelial cell: insights from a computational model.

Matrix metalloproteinases (MMPs) are a family of enzymes responsible for the proteolytic processing of extracellular matrix (ECM) structural proteins under physiological and pathological conditions. During sprouting angiogenesis, the MMPs expressed by a single "tip" endothelial cell exhibit proteolytic activity that allows the cells of the sprouting vessel bud to migrate into the ECM. Membrane type I matrix metalloproteinase (MT1-MMP) and the diffusible matrix metalloproteinase MMP2, in the presence of the tissue inhibitor of metalloproteinases TIMP2, constitute a system of proteins that play an important role during the proteolysis of collagen type I matrices. Here, we have formulated a computational model to investigate the proteolytic potential of such a tip endothelial cell. The cell expresses MMP2 in its proenzyme form, pro-MMP2, as well as MT1-MMP and TIMP2. The interactions of the proteins are described by a biochemically detailed reaction network. Assuming that the rate-limiting step of the migration is the ability of the tip cell to carry out proteolysis, we have estimated cell velocities for matrices of different collagen content. The estimated velocities of a few microns per hour are in agreement with experimental data. At high collagen content, proteolysis was carried out primarily by MT1-MMP and localized to the cell leading edge, whereas at lower concentrations, MT1-MMP and MMP2 were found to act in parallel, causing proteolysis in the vicinity of the leading edge. TIMP2 is a regulator of the proteolysis localization because it can shift the activity of MT1-MMP from its enzymatic toward its activatory mode, suggesting a tight mechanosensitive regulation of the enzymes and inhibitor expression. The model described here provides a foundation for quantitative studies of angiogenesis in extracellular matrices of different compositions, both in vitro and in vivo. It also identifies critical parameters whose values are not presently available and which should be determined in future experiments.

Cell Movement↗

Study of spinal cord evoked injury potential by use of computer modeling and in dogs with naturally acquired thoracolumbar spinal cord compression.

OBJECTIVE: To add objective measurements of the characteristics of evoked injury potentials (EIP) and their relations to clinical severity in dogs with thoracolumbar spinal cord damage. ANIMALS: 25 dogs with naturally acquired spinal cord compression attributable to disk extrusion or vertebral fracture at the level of the thoracolumbar junction and with various degrees of paresis/paralysis. PROCEDURE: Spinal cord potentials evoked by tibial nerve stimulation were recorded every 5 to 10 mm at the lamina level in the vicinity of the cord compression. This allowed an EIP to be recorded even in the least handicapped dogs. A computer model yielded information about the waveform changes of the EIP in the vicinity of conduction blocks. RESULTS: The EIP waveform changed from biphasic to monophasic a short distance caudad to the location of spinal cord compression. Location of a maximal conduction block was measured in relation to position of the electrodes recording this waveform change. The distance between the assumed conduction block and the actual spinal cord compression was larger in the most affected dogs. The amplitude of the EIP was not related to severity of the clinical picture; however, the proximity of the recording electrode to the spine influenced the amplitude and the waveform of the EIP. CONCLUSION AND CLINICAL RELEVANCE: Change in the EIP waveform from biphasic to monophasic makes it possible to estimate the conduction block location along the spinal cord. A large distance between the assumed conduction block and site of actual cord compression could be an objective argument to confirm severity of a lesion.

Animals↗

Application of discrete computer modeling to the dynamics of cell populations.

Computer simulation of the dynamics of cell populations is discussed in the paper. Fundamental features of CELLSIM simulation language are described. Acute leukemia cell population model was implemented on a digital computer using CELLSIM. Some of the computer outputs are shown. The results obtained so far are interesting and promissing for further research.

Acute Disease↗

[Stereoscopic simulation of ear surgery intervention with a novel 3D computer models].

INTRODUCTION: The presentation of the surgical anatomy of the temporal bone by standard anatomical figures is not suitable for otosurgical training. For the comprehension of its complex morphology temporal bone drilling is inalienable. Aim of the present cooperation was to gain an interactive real-3D program for the simulation of specific laterobasal surgical approaches. METHODS: The program was derived from a standard horizontal section of a human temporal bone using a Siemens Somatom Plus 4 Tomograph. The slice thickness was 1 mm, the image matrix was 512. The Voxel-Man-system was used to built up the application for unix workstations. RESULTS: Each step of a surgical approach to the temporal bone can be performed by the present computer model. Calculation in a stereo mode even allows spatial 3D-perception when using red/green glasses. CONCLUSIONS: This program is a novel tool to simulate critical aspects of otosurgical procedures on a computer. Up to now the lack of tactile and kinesthetic information does not allow to renounce individual temporal bone drilling.

Computer Simulation↗

Computational modelling of ErbB family phosphorylation dynamics in response to transforming growth factor alpha and heregulin indicates spatial compartmentation of phosphatase activity.

Members of the ErbB receptor family are associated with several cancers and appear to be providing useful targets for pharmacological therapeutics for tumours of the lung and breast. Further improvements of these therapies may be guided by a quantitative, dynamic integrative systems understanding of the complexities of ErbB dimerisation, trafficking and activation, for it is these complexities that render difficult intuiting how perturbations such as drug intervention will affect ErbB signalling activities. Towards this goal, we have developed a computational model implementing commonly accepted principles governing ErbB receptor interaction, trafficking, phosphorylation and dephosphorylation. Using this model, we are able to investigate several hypotheses regarding the compartmental localisation of dephosphorylation. Model results applied to experimental data on ErbB 1, ErbB2 and ErbB3 phosphorylation in H292 human lung carcinoma cells support a hypothesis that key dephosphorylation activity for these receptors occurs largely in an intracellular, endosomal compartment rather than at the cell surface plasma membrane. Thus, the endocytic trafficking-related compartmentalisation of dephosphorylation may define a critical aspect of the ErbB signalling response to ligand.

Cell Line↗

Comparison of retrospective luminescence dosimetry with computational modeling in two highly contaminated settlements downwind of the Chernobyl NPP.

The cumulative absorbed dose in bricks collected from six buildings in two heavily contaminated settlements (137Cs > 2,000 kBq m(-2)) located downwind of the Chernobyl Nuclear Power Plant was determined using luminescence techniques by six laboratories. The settlements, Vesnianoje in Ukraine and Zaborie in Russia, are located in, respectively, proximal and distal locations relative to the Chernobyl Nuclear Power Plant. The luminescence determinations of cumulative dose in brick, after subtraction of the natural background dose, were translated to absorbed dose in air at a Reference Location using conversion factors derived from Monte Carlo simulations of photon transport. The simulations employed source distributions inferred from contemporary soil contamination data and also took into account heterogeneity of fallout deposition. This translation enables the luminescence determinations to be compared directly with values of cumulative absorbed dose obtained by computational modeling and also other dose reconstruction methods. For each sampled location the cumulative dose was calculated using three deterministic models, two of which are based on the attenuation of dose-rate with migration of radionuclides in soil and the third on historic instrumental gamma dose-rate data. The results of the comparison of the two methods indicate overall agreement within margins of +/-25%. The methodology developed is generally applicable and adaptable to areas contaminated by much lower levels of radioactive fallout in which brick buildings are found.

Computer Simulation↗