Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “Models, Cardiovascular”

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 487 records · Page 27Linked to original sources

Structural effects of an LQT-3 mutation on heart Na+ channel gating.

Computational methods that predict three-dimensional structures from amino acid sequences have become increasingly accurate and have provided insights into structure-function relationships for proteins in the absence of structural data. However, the accuracy of computational structural models requires experimental approaches for validation. Here we report direct testing of the predictions of a previously reported structural model of the C-terminus of the human heart Na(+) channel. We focused on understanding the structural basis for the unique effects of an inherited C-terminal mutation (Y1795C), associated with long QT syndrome variant 3 (LQT-3), that has pronounced effects on Na(+) channel inactivation. Here we provide evidence that this mutation, in which a cysteine replaces a tyrosine at position 1795 (Y1795C), enables the formation of disulfide bonds with a partner cysteine in the channel. Using the predictions of the model, we identify the cysteine and show that three-dimensional information contained in the sequence for the channel protein is necessary to understand the structural basis for some of the effects of the mutation. The experimental evidence supports the accuracy of the predicted structural model of the human heart Na(+) channel C-terminal domain and provides insight into a structural basis for some of the mutation-induced altered channel function underlying the disease phenotype.

Amino Acid Substitution↗

Explaining glomerular pores with fiber matrices. A visualization study based on computer modeling.

The extracellular space of the glomerular capillary wall is occupied by a complex meshwork of fibrous molecules. Little is understood about how the size, shape, and charge recognition properties of glomerular ultrafiltration arise from this space-filling fiber matrix. We studied the problem of size recognition by visualizing the void volume accessible to hard spheres in computer-generated three-dimensional homogeneous random fiber matrices. The spatial organization of the void volume followed a complex "blob-and-throat" pattern in which circumscribed cavities of free space within the matrix ("blobs") were joined to adjacent cavities by narrower throats of void space. For sufficiently small solutes, chains of blobs and throats traversed the matrix, providing pathways for trans-matrix permeation. The matrices showed threshold or gating properties with respect to permeation: solutes whose radius exceeded a critical value, at which a throat on the last connected trans-matrix pathway pinched off, could not cross, whereas smaller solutes had nonzero permeability. The thresholds may give the glomerular fiber matrix porelike response properties and explain why pore models have been such a useful means of treating permselectivity.

Capillaries↗

Deformation and flow of red blood cells in a synthetic lattice: evidence for an active cytoskeleton.

We introduce the use of microfabrication techniques to construct on a silicon wafer a synthetic capillary bed with 2.5- to 4-micron (mu)-wide channels. Establishment of a fluid pressure gradient allowed us to observe simultaneously using optical microscopy hundreds of cells flowing through the bed at physiological speeds. We find a large distribution of mobilities among red cells flowing through the structure; smaller channels provide a greater impedance to flow than larger ones, indicating that kinetic drag variations provide the origin of the distribution. The mobility of a particular cell is not correlated with the cell diameter but appears to be inversely correlated with intracellular calcium concentration of the cell, as determined by fluorescence of the calcium-binding dye fluo-3 AM. Also, we are able to use the parallel processing nature of our arrays to observe isolated events where the rigidity of the red cell seems to change suddenly over several orders of magnitude as it blocks a channel in the array.

Capillaries↗

Factors shaping the confocal image of the calcium spark in cardiac muscle cells.

The interpretation of confocal line-scan images of local [Ca2+]i transients (such as Ca2+ sparks in cardiac muscle) is complicated by uncertainties in the position of the origin of the Ca2+ spark (relative to the scan line) and by the dynamics of Ca(2+)-dye interactions. An investigation of the effects of these complications modeled the release, diffusion, binding, and uptake of Ca2+ in cardiac cells (producing a theoretical Ca2+ spark) and image formation in a confocal microscope (after measurement of its point-spread function) and simulated line-scan images of a theoretical Ca2+ spark (when it was viewed from all possible positions relative to the scan line). In line-scan images, Ca2+ sparks that arose in a different optical section or with the site of origin displaced laterally from the scan line appeared attenuated, whereas their rise times slowed down only slightly. These results indicate that even if all Ca2+ sparks are perfectly identical events, except for their site of origin, there will be an apparent variation in the amplitude and other characteristics of Ca2+ sparks as measured from confocal line-scan images. The frequency distributions of the kinetic parameters (i.e., peak amplitude, rise time, fall time) of Ca2+ sparks were calculated for repetitive registration of stereotyped Ca2+ sparks in two experimental situations: 1) random position of the scan line relative to possible SR Ca(2+)-release sites and 2) fixed position of the scan line going through a set of possible SR Ca(2+)-release sites. The effects of noise were incorporated into the model, and a visibility function was proposed to account for the subjective factors that may be involved in the evaluation of Ca(2+)-spark image parameters from noisy experimental recordings. The mean value of the resulting amplitude distributions underestimates the brightness of in-focus Ca2+ sparks because large numbers of out-of-focus Ca2+ sparks are detected (as small Ca2+ sparks). The distribution of peak amplitudes may split into more than one subpopulation even when one is viewing stereotyped Ca2+ sparks because of the discrete locations of possible SR Ca(2+)-release sites in mammalian ventricular heart cells.

Animals↗

Numerical simulation of local calcium movements during L-type calcium channel gating in the cardiac diad.

Computer simulation was used to investigate the calcium levels after sarcolemmal calcium influx through L-type calcium channels (DHPRs) into the narrow diadic space of cardiac muscle. The effect of various cytosolic and membranebound buffers, diad geometry, DHPR properties (open time and current), and surface charge were examined. The simulations showed that phospholipid binding sites on the sarcolemmal membrane are the major buffer affecting free calcium ([Ca2+]) levels in the diad. The inclusion of surface charge effects calculated from Gouy-Chapman theory resulted in a marked decrease in [Ca2+] levels at all times and a faster decay of [Ca2+] after termination of DHPR influx. For a DHPR current of 200 fA, [Ca2+] at the center of the diad reached peak levels of approximately 73 microM. In larger diads (> or = 400 nm diameter), [Ca2+] decayed more slowly than in smaller diads (100-200 nm diameter), although peak [Ca2+] levels reached during typical DHPR open times were similar. For a wide range of DHPR single-channel current magnitudes (Ica = 25-200 fA), [Ca2+] levels in the diad were approximately proportional to ICa. The decrease in calculated [Ca2+] levels due to the effects of surface charge can be interpreted as resulting from an effective "volume expansion" of the diad space. Furthermore, the layer of increased [Ca2+] close to the sarcolemmal membrane can act as a fast buffer.

Binding Sites↗

Calcium concentration and movement in the ventricular cardiac cell during an excitation-contraction cycle.

This paper extends the model for Ca movement in the cardiac ventricular cell from the diadic cleft space to the entire sarcomere. The model predicts the following: 1) Shortly after SR release there is a [Ca] gradient >3 orders of magnitude from cleft center to M-line which, 50 ms after release, is still >30. Outside the cleft, 40 ms after cessation of release, the axial gradient from Z to M-line is >3. 2) At the end of SR release, >50% of the total Ca released is bound to low-affinity inner sarcolemmal phospholipid binding sites within the cleft. 3) Halving the SR release almost doubles the fraction of release removed from the cell via Na/Ca exchange and reduces average sarcomeric free [Ca] by 70%. 4) Adding 100 microM fluo-3, which doubles the buffering capacity of the cytoplasm, reduces peak average sarcomeric [Ca] by >50% and increases the initial half-time for [Ca] decrease by approximately twofold. 5) A typical Ca "spark" can be generated by an SR release 20% of maximum (4 x 10(-20) moles) over 2 ms. Fluo-3 (100 microM) significantly "shrinks" the spark. 6) The "spark" is a consequence of elementary events within the diadic cleft space. For example, removal of cleft binding sites would cause average sarcomeric Ca to increase by >10 fold, fall 10 times more rapidly, decrease latency for appearance of the spark by >20 times, and reduce spark duration by 85%. 7) Dividing SR Ca release between cleft and corbular SR produces a secondary [Ca] peak and a "flattening" of the sarcomeric [Ca] transient. These changes probably could not be resolved with current confocal microscopic techniques.

Aniline Compounds↗

Cardiac Ca2+ dynamics: the roles of ryanodine receptor adaptation and sarcoplasmic reticulum load.

We construct a detailed mathematical model for Ca2+ regulation in the ventricular myocyte that includes novel descriptions of subcellular mechanisms based on recent experimental findings: 1) the Keizer-Levine model for the ryanodine receptor (RyR), which displays adaptation at elevated Ca2+; 2) a model for the L-type Ca2+ channel that inactivates by mode switching; and 3) a restricted subspace into which the RyRs and L-type Ca2+ channels empty and interact via Ca2+. We add membrane currents from the Luo-Rudy Phase II ventricular cell model to our description of Ca2+ handling to formulate a new model for ventricular action potentials and Ca2+ regulation. The model can simulate Ca2+ transients during an action potential similar to those seen experimentally. The subspace [Ca2+] rises more rapidly and reaches a higher level (10-30 microM) than the bulk myoplasmic Ca2+ (peak [Ca2+]i approximately 1 microM). Termination of sarcoplasmic reticulum (SR) Ca2+ release is predominately due to emptying of the SR, but is influenced by RyR adaptation. Because force generation is roughly proportional to peak myoplasmic Ca2+, we use [Ca2+]i in the model to explore the effects of pacing rate on force generation. The model reproduces transitions seen in force generation due to changes in pacing that cannot be simulated by previous models. Simulation of such complex phenomena requires an interplay of both RyR adaptation and the degree of SR Ca2+ loading. This model, therefore, shows improved behavior over existing models that lack detailed descriptions of subcellular Ca2+ regulatory mechanisms.

Action Potentials↗

Analysis of the mechanisms of mitochondrial NADH regulation in cardiac trabeculae.

We have previously shown that increased cardiac work initially caused a rapid Ca(2+)-independent fall of mitochondrial [NADH] ([NADH](m)) to a minimum level, and this was followed by a slow Ca(2+)-dependent recovery toward control level (Brandes and Bers, Biophys. J. 71:1024-1035, 1996; Brandes and Bers, Circ. Res. 80:82-87, 1997). The purpose of this study is to improve our understanding of the factors that control [NADH](m) during increased work. [NADH](m) was monitored using fluorescence spectroscopy in intact rat trabeculae isolated from the right ventricular wall. Work was increased by increasing sarcomere length, pacing frequency, external [Ca(2+)], or by decreased temperature. The results were: 1) The initial fall of [NADH](m) during increased pacing frequency depends independently on increased myofilament work and on increased Ca(2+)-transport ATPase activity. 2) The [NADH](m) recovery process depends on average cytosolic [Ca(2+)] (Av[Ca(2+)](c)), but not on absolute work level. 3) The initial fall of [NADH](m) and the [NADH](m) recovery are similar whether increased work is associated with low frequency and high Ca(2+)-transient amplitude or vice versa (at the same myofilament work level and Av[Ca(2+)](c)). 4) The mechanisms associated with the smaller fall and recovery of [NADH](m) at 37 degrees C versus 27 degrees C, may be explained by lowered Av[Ca(2+)](c) and myofilament work. The NADH control mechanisms that operate at lower temperature are thus qualitatively similar at more physiological temperatures.

Adenosine Triphosphate↗

6-[18F]fluoro-L-m-tyrosine: metabolism, positron emission tomography kinetics, and 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine lesions in primates.

The tracer 6-[18F]fluoro-L-m-tyrosine (FMT) was studied with regard to its biochemistry and kinetics, as well as its utility in evaluating brain dopaminergic function in primates before and after 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) treatment using positron emission tomography (PET). Plasma analysis of FMT and its F18-labeled metabolites 6-fluoro-3-hydroxyphenylacetic acid (FPAC) and 6-fluoro-3-hydroxyphenylethylamine (FMA) during PET scanning enabled kinetic analysis of FMT uptake. A separate study examined brain FMT metabolism in MPTP-naive monkeys euthanized 60 or 120 min after FMT injection. Almost 60% of total plasma F-18 activity was associated with FPAC and FMA 120 min after FMT injection. The FMT signal accumulated preferentially in dopaminergic areas such as caudate and putamen. This bilateral FMT signal was disrupted after unilateral intracarotid artery (ICA) MPTP infusion which reduced ipsilateral striatal activity. A three compartment three kinetic rate constant model for FMT uptake revealed reduced FMT decarboxylation (k3) in ipsilateral caudate and putamen after unilateral MPTP although a further decrease was not evident after intravenous MPTP. FPAC was the major F-18 species in all brain regions except in cerebellum where FMT was predominant 60 min post-mortem. FPAC was most concentrated in dopaminergic areas whereas lower levels occurred in areas containing few dopamine terminals. These data demonstrate preferential FMT metabolism and F-18 retention in dopaminergic tissue and support the use of FMT to evaluate normal and abnormal dopaminergic function.

Animals↗

Model-based assessment of pressure and flow-dependent coronary responses following abrupt pressure drops.

The response of the coronary vasculature to an experimental manoeuvre of step-like decrease of the perfusion pressure has been investigated with a model. The coronary vasculature was simulated using a 'windkessel' scheme. Proximal resistance and compliance were assumed to be pressure-independent. The distal resistance, on the contrary, was controlled by a feed-back loop which accounts for the smooth muscle activation induced by the pressure variation. Three more parameters were introduced, and namely the smooth muscle activation time constant and the pressure-induced and flow-induced gains. The parameter values were assessed by comparing the model predicted coronary flow with the one actually measured in animals.

Animals↗

Numerical modeling of the flow in stenosed coronary artery. The relationship between main hemodynamic parameters.

The severity of coronary arterial stenosis is usually measured by either simple geometrical parameters, such as percent diameter stenosis, or hemodynamically based parameters, such as the fractional flow reserve (FFR) or coronary flow reserve (CFR). The present study aimed to establish a relationship between actual hemodynamic conditions and the parameters that define stenosis severity in the clinical setting. We used a computational model of the blood flow in a vessel with a blunt stenosis and an autoregulated vascular bed to simulate a stenosed blood vessel. A key point in creating realistic simulations is to properly model arterial autoregulation. A constant flow regulation mechanism resulted in CFR and FFR values that were within the physiological range, while a constant wall-shear stress model yielded unrealistic values. The simulation tools developed in the present study may be useful in the clinical assessment of single and multiple stenoses by means of minimally invasive methods.

Blood Flow Velocity↗

Simulation of acute subdural hematoma and diffuse axonal injury in coronal head impact.

Coronal head impacts were simulated in a physical model, based on the hypothesis that acute subdural hematoma (ASDH) is related to cerebral vertex displacement and diffuse axonal injury (DAI) to local Green-Lagrange strain. The geometry of the 2D model was based on anatomical measurements taken from the MRI scans of 10 adult males. Silicone gel modelled the cerebrum, paraffin the CSF and elastic membranes the trabeculae of the sulci. Pendulum impacts gave peak angular acceleration of 7800 rad s(-2) in models with and without sulci. The motion of the gel and Green-Lagrange strain were calculated from tracked coordinates of Patrick markers. Worst-case bridging vein strains are produced on the contrecoup side and are approximately doubled by adding sulci. Given that axons in the corpus callosum are highly oriented, Green-Lagrange strain was resolved in the fibre direction. It is found to be close to the minimum principal strain, indicating a degree of natural, teleological protection for the axons. The data support the use of delta0peak as a suitable descriptor for the risk of DAI but not for ASDH.

Adult↗

A numerical simulation of mechanical heart valve closure fluid dynamics.

A computational fluid dynamics model for the analysis of the bileaflet mechanical heart valve closure process is presented. The objective of the study is to demonstrate the ability of the numerical model to simulate the leaflet motion during the closing phase in order to investigate the closure fluid dynamics and to evaluate the effect of alterations in the leaflet tip geometry. The model has been applied to six different combinations of the leaflet tip geometry and the gap width between the leaflet tip and the housing. The results show that the negative pressure quickly develops on the atrial side of the leaflet tip. The pressure becomes more negative as the leaflet closure progresses and the lowest pressure is reached before the leaflet comes to a stop in the closed position. The flow dynamics at the instant of valve closure is strongly dependent on the leaflet velocity during the closing phase. Decrease of the tip velocity by a factor of three in the last four degrees of leaflet motion leads to a 50% reduction in the negative pressure magnitude.

Biomechanical Phenomena↗

Role of tapering in aortic wave reflection: hydraulic and mathematical model study.

Pressure and flow have been measured simultaneously at six locations along the aorta of an anatomically correct 1:1 scale hydraulic elastic tube model of the arterial tree. Our results suggest a discrete reflection point at the level of the renal arteries based on (i) the quarter-wavelength formula and (ii) the comparison of foot-to-foot (c(ff)) and apparent phase velocity (c(app)). However, separation of the pressure wave into an incident and reflected wave at all six locations indicates continuous reflection: a reflected wave is generated at each location as the forward wave passes by. We did a further analysis using a mathematical transmission line model with a simple tapering geometry (length 50 cm, 31 and 11 mm proximal and distal diameter, respectively) for a low (0.32 ml/mmHg), normal (1.6 ml mmHg) and high (8 ml/mmHg) value of total arterial compliance. Using the quarter-wavelength formula, a discrete reflection point is found at x = 33 cm, the level of the renal arteries, independent of the value of total compliance. However, local analysis comparing c(ff) and c(app) does not reveal a marked reflection site, and the analysis of incident and reflected waves merely suggests a continuous reflection. We therefore conclude that the measured in vivo aortic wave reflection indices are the result of at least two interacting phenomena: a continuous wave reflection due to tapering, and local reflections arising from branches at the level of the diaphragm. The continuous reflection is hidden in the input impedance pattern. Using the quarter-wavelength formula or the classical wave separation theory, it appears as a reflection coming from a single discrete site, confusingly also located at the level of the diaphragm. Therefore, the quarter-wavelength formula and the linear wave separation theory should be used with caution to identify wave reflection zones in the presence of tapering, i.e., in most mammalian arteries.

Animals↗

ST-deviation reconstruction in missing leads on the 12-lead ECG: applicability in studies on ST-segment resolution during thrombolysis.

Quantitative analysis of ST-segment deviations (STdev) and their resolution by treatment (STR; calculated from a combined sum of STdev in multiple leads) are used in trials on reperfusion for myocardial infarction (MI). Unreadable or unavailable electrocardiogram (ECG) leads are a common reason for exclusion, decreasing the statistical power of the trials. We developed mathematical formulas for reconstruction of immeasurable STdev based on STdev from other available leads on the 12-lead ECG. Formulas were deducted from a database of computer-assisted STdev measurements in 2 ECGs (baseline and 180 min after thrombolysis) of 1121 pts. Their accuracy was later evaluated on a second dataset of 377 pts. Acceptable fits could be derived for absent single leads, or for groups of absent limb leads (I-II-III or aVL-aVF). The intraclass correlation coefficient between real and calculated STdev was >or= 0.80 for each (0.77 for V1 in inferior MI). The correlations between STR calculated from original data and from reconstructed STdev were very strong (all intraclass correlation >or=0.97), and discordance in STR subgroup categorization occurred in <or=10% of pts in all but one of the scenarios (I-II-III substituted in 180 min ECG in inferior MI). Scenarios with multiple missing precordial leads however are not substitutable, nor are calculated STdev reliable for STR evaluation in only the lead with highest elevation in baseline. STdev reconstruction formulas can reliably be used in trials where analysis of aggregate STR is an endpoint. Reliable substitution can significantly increase the number of evaluable patients and therefore strengthen the statistical power of these trials.

Electrocardiography↗

Modelling the periodicity of cardiac muscle.

The intractable problem of modelling cardiac muscle of arbitrary extent while preserving cellular structure has been solved using an analytical rather than numerical approach with a method called two-scale asymptotic analysis. In this method, the myocardium was modelled as a collection of bundles arranged periodically in space and connected by junctions, and the distribution of the steady-state potential and current density was determined. The potential both along and across fibers was found to contain a distinct periodic component that determines the transmembrane potential. The magnitude of the transmembrane potential depends on the gradient of applied potential, the dimensions of the bundles, and their internal conductivity. Current flows primarily in the extracellular space, and the extracellular pathway also determines the apparent conductivity of cardiac muscle.

Electric Conductivity↗

Numerical analysis of electrical defibrillation. The parallel approach.

Numerical modeling offers a viable tool for studying electrical defibrillation, allowing the behavior of field quantities to be observed easily as the different system parameters are varied. One numerical technique, namely the finite-element method, has been found particularly effective for modeling complex thoracic anatomies. However, an accurate finite-element model of the thorax often requires a large number of elements and nodes, leading to a large set of equations that cannot be solved effectively with the computational power of conventional computers. This is especially true if many finite-element solutions need to be achieved within a reasonable time period (eg, electrode configuration optimization). In this study, the use of massively parallel computers to provide the memory and reduction in solution time for solving these large finite-element problems is discussed. Both the uniform and unstructured grid approaches are considered. Algorithms that allow efficient mapping of uniform and unstructured grids to data-parallel and message-passing parallel computers are discussed. An automatic iterative procedure for electrode configuration optimization is presented. The procedure is based on the minimization of an objective function using the parallel direct search technique. Computational performance results are presented together with simulation results.

Algorithms↗