PubMed Health⌕ Search

Biomedical subjects

M Vauhkonen

Publications and source records attributed to M Vauhkonen.

At least 19 recordsLinked to original sources

A MATLAB package for the EIDORS project to reconstruct two-dimensional EIT images.

The EIDORS (electrical impedance and diffuse optical reconstruction software) project aims to produce a software system for reconstructing images from electrical or diffuse optical data. MATLAB is a software that is used in the EIDORS project for rapid prototyping, graphical user interface construction and image display. We have written a MATLAB package (http://venda.uku.fi/ vauhkon/) which can be used for two-dimensional mesh generation, solving the forward problem and reconstructing and displaying the reconstructed images (resistivity or admittivity). In this paper we briefly describe the mathematical theory on which the codes are based on and also give some examples of the capabilities of the package.

Algorithms↗

A new computational approach for cortical imaging.

Estimation of current or potential distribution on the cortex is used to obtain information about neural sources from the scalp recorded electroencephalogram. If the active sources in the brain are superficial, the estimated field distribution on the cortex also yields information about the active source configuration. In these cases, these methods can be used as source localization methods. In this study, we concentrate on finite-element-based cortex potential estimation. Usually these methods require surface interpolation of the recorded voltages at the electrodes onto the entire scalp surface. We propose a new computational approach which does not require the use of surface interpolation but does it implicitly and uses only the recorded data at the electrodes. We refer to this method as the systematic approach (SA). We compare the SA with the surface interpolation approach (IA) and show that the SA is able to produce somewhat better accuracy than the IA. However, the main asset is that the sensitivity of the cortical potential maps to the regularization parameter is significantly lower than with the IA.

Brain Mapping↗

Effects of electrode properties on EEG measurements and a related inverse problem.

A trend in EEG measurements is to increase the number of measurement electrodes in order to improve the spatial resolution of the recorded voltage distribution at the scalp. It is assumed that this would implicate better accuracy in the EEG inverse estimates. However, this does not necessarily hold. The reason for this is that the electrodes create a well conducting shunting "layer" on the scalp which affects the voltage distribution. This may decrease the information obtained and may therefore worsen the inverse estimates. Electrodes in EEG inverse problems are commonly modeled as point electrodes. This model cannot take into account the possible shunting effect of the electrodes. In this study the measurement electrodes are modeled using the so-called complete electrode model which takes into account the actual size of the electrode, the contact impedance between the skin and the electrode and also the shunting effect of the electrodes. In this paper the effects of the electrode size and the contact impedance on the voltage distribution are studied by simulations. It is shown that, depending on the size and the contact impedance of the electrodes, increasing the number of electrodes does not necessarily improve the accuracy of the inverse estimates. We also conclude that the use of the point electrode model is quite adequate in normal EEG studies. The use of a complete electrode model is necessary if electrodes cover more than 50% of the surface area.

Electric Impedance↗

Simultaneous reconstruction of internal tissue region boundaries and coefficients in optical diffusion tomography.

In this paper we propose a new numerical method to the inverse problem in optical diffusion tomography. We consider the reconstruction of the diffusion and absorption coefficients (kappa, mu(a)) within a domain omega which is known to consist of a set of disjoint regions of distinct tissue types. The assumption is that the regions of different tissues are bounded by smooth boundary curves and have constant absorption and diffusion coefficients. The goal in the proposed method is to reconstruct simultaneously the boundaries of the tissue regions together with the absorption and diffusion coefficients within these regions. The solution of the problem is based on the finite element method and subdivision of the elements. The performance of the proposed method is evaluated by simulations in which the optical parameters (kappa, mu(a)) are relevant in medical applications of optical tomography. It is shown that the proposed method is able to recover both the boundaries and the coefficients with good accuracy.

Image Processing, Computer-Assisted↗

Sensitivity matrix and reconstruction algorithm for EIT assuming axial uniformity.

In electrical impedance tomography (EIT) two-dimensional models continue to be applied despite their known inability to provide correct reconstruction. In this paper, a reconstruction algorithm that assumes a translationally invariant conductivity distribution is described. A more precise forward solver is obtained by taking off-slice currents into consideration. An appropriate sensitivity matrix is derived. Numerical evidence for the improvement in precision compared to two-dimensional reconstruction is given.

Algorithms↗

Errors due to the truncation of the computational domain in static three-dimensional electrical impedance tomography.

In electrical impedance tomography (EIT), an approximation for the internal resistivity distribution is computed based on the knowledge of the injected currents and measured voltages on the surface of the body. The currents spread out in three dimensions and therefore off-plane structures have a significant effect on the reconstructed images. A question arises: how far from the current carrying electrodes should the discretized model of the object be extended? If the model is truncated too near the electrodes, errors are produced in the reconstructed images. On the other hand if the model is extended very far from the electrodes the computational time may become too long in practice. In this paper the model truncation problem is studied with the extended finite element method. Forward solutions obtained using so-called infinite elements, long finite elements and separable long finite elements are compared to the correct solution. The effects of the truncation of the computational domain on the reconstructed images are also discussed and results from the three-dimensional (3D) sensitivity analysis are given. We show that if the finite element method with ordinary elements is used in static 3D EIT, the dimension of the problem can become fairly large if the errors associated with the domain truncation are to be avoided.

Biometry↗

Effects of local skull inhomogeneities on EEG source estimation.

The accuracy of the head model affects the solutions of the EEG inverse problems. If a simple three-sphere model and standard conductivity values for brain, skull and scalp regions are used, significant errors may occur in the dipole localisation. One of the most sensitive head model parameters is the conductivity of the skull. A realistic three-dimensional finite-element model provides a method to study the effect of inhomogeneities of the skull on the solutions of EEG inverse problems. In this paper the effect of a local skull conductivity inhomogeneity on source estimation accuracy is analyzed by computer simulations for different numbers of electrodes. It is shown that if the inhomogeneity of the skull conductivity is not taken into account, localisation errors of approximately 1 cm can be encountered in the equivalent current dipole estimation. This modelling error introduces a bias to the solution which cannot be compensated by increasing the number of electrodes.

Biophysical Phenomena↗

Subspace regularization method for the single-trial estimation of evoked potentials.

A method for the single-trial estimation of the evoked potentials is proposed. The method is based on the so-called subspace regularization approach in which the second-order statistics of the set of the measurements is used to form a prior information model for the evoked potentials. The method is closely related to the Bayesian estimation. The performance of the proposed method is evaluated using realistic simulations. As a specific application the method is applied to the estimation of the target responses in the P300 test.

Bayes Theorem↗

Three-dimensional electrical impedance tomography based on the complete electrode model.

In electrical impedance tomography an approximation for the internal resistivity distribution is computed based on the knowledge of the injected currents and measured voltages on the surface of the body. It is often assumed that the injected currents are confined to the two-dimensional (2-D) electrode plane and the reconstruction is based on 2-D assumptions. However, the currents spread out in three dimensions and, therefore, off-plane structures have significant effect on the reconstructed images. In this paper we propose a finite element-based method for the reconstruction of three-dimensional resistivity distributions. The proposed method is based on the so-called complete electrode model that takes into account the presence of the electrodes and the contact impedances. Both the forward and the inverse problems are discussed and results from static and dynamic (difference) reconstructions with real measurement data are given. It is shown that in phantom experiments with accurate finite element computations it is possible to obtain static images that are comparable with difference images that are reconstructed from the same object with the empty (saline filled) tank as a reference.

Electric Conductivity↗

Iodine-123 labeled nor-beta-CIT as a potential tracer for serotonin transporter imaging in the human brain with single-photon emission tomography.

Iodine-123 labelled 2beta-carbomethoxy-3beta-(4-iodophenyl) (nor-beta-CIT) is an analogue of beta-CIT, which has high affinity to the serotonin transporter. Initial single-photon emission tomography (SPET) studies with [123I]nor-beta-CIT were performed in five healthy volunteers. In addition, its metabolism in plasma was investigated with gradient high performance liquid chromatography. [123I]nor-beta-CIT was prepared by a method which gave a specific radioactivity of more than 180 GBq/micromol. Unchanged [123I]nor-beta-CIT in plasma accounted for 43% and 19% of total radioactivity after 30 and 180 min, respectively. The dynamic SPET studies demonstrated a high and rapid uptake of radioactivity in the brain (6%/ID at 30 min). Highest accumulation was observed in the striatum, the mid-brain and the thalamus. The specific binding in the mid-brain was 33% higher compared with that of [123I]beta-CIT. The high radioactivity in the mid-brain is assumed to represent the accumulation of [123I]nor-beta-CIT in the serotonin transporter-rich regions, which indicates that [123I]nor-beta-CIT might be a potential tracer for visualization of serotonin transporter sites in the human brain with SPET.

Adult↗

A Kalman filter approach to track fast impedance changes in electrical impedance tomography.

In electrical impedance tomography (EIT), an estimate for the cross-sectional impedance distribution is obtained from the body by using current and voltage measurements made from the boundary. All well-known reconstruction algorithms use a full set of independent current patterns for each reconstruction. In some applications, the impedance changes may be so fast that information on the time evolution of the impedance distribution is either lost or severely blurred. In this paper, we propose an algorithm for EIT reconstruction that is able to track fast changes in the impedance distribution. The method is based on the formulation of EIT as a state-estimation problem and the recursive estimation of the state with the aid of the Kalman filter. The performance of the proposed method is evaluated with a simulation of human thorax in a situation in which the impedances of the ventricles change rapidly. We show that with optimal current patterns and proper parameterization, the proposed approach yields significant enhancement of the temporal resolution over the conventional reconstruction strategy.

Algorithms↗

Tikhonov regularization and prior information in electrical impedance tomography.

The solution of impedance distribution in electrical impedance tomography is a nonlinear inverse problem that requires the use of a regularization method. The generalized Tikhonov regularization methods have been popular in the solution of many inverse problems. The regularization matrices that are usually used with the Tikhonov method are more or less ad hoc and the implicit prior assumptions are, thus, in many cases inappropriate. In this paper, we propose an approach to the construction of the regularization matrix that conforms to the prior assumptions on the impedance distribution. The approach is based on the construction of an approximating subspace for the expected impedance distributions. It is shown by simulations that the reconstructions obtained with the proposed method are better than with two other schemes of the same type when the prior is compatible with the true object. On the other hand, when the prior is incompatible with the true object, the method will still give reasonable estimates.

Algorithms↗

Assessment of errors in static electrical impedance tomography with adjacent and trigonometric current patterns.

In electrical impedance tomography (EIT), difference imaging is often preferred over static imaging. This is because of the many unknowns in the forward modelling which make it difficult to obtain reliable absolute resistivity estimates. However, static imaging and absolute resistivity values are needed in some potential applications of EIT. In this paper we demonstrate by simulation the effects of different error components that are included in the reconstruction of static EIT images. All simulations are carried out in two dimensions with the so-called complete electrode model. Errors that are considered are the modelling error in the boundary shape of an object, errors in the electrode sizes and localizations and errors in the contact impedances under the electrodes. Results using both adjacent and trigonometric current patterns are given.

Computer Simulation↗

Degradation of pyrene-labelled phospholipids by lysosomal phospholipases in vitro. Dependence of degradation on the length and position of the labelled and unlabelled acyl chains.

The hydrolysis of pyrenylacyl phosphatidylcholines (PyrnPCs)(n indicates the number of aliphatic carbons in the pyrene-chain) by crude lysosomal phospholipases in vitro was investigated. PyrnPCs consist of several sets in which the length of the pyrene-labelled or the unlabelled acyl chain, linked to the sn-1 or sn-2 position, was systematically varied. Lysophosphatidylcholine and fatty acid were the only fluorescent breakdown products detected, thus indicating that PyrnPCs were degraded by A-type phospholipases and lysophospholipases. Of these, mainly A1-type phospholipases appear to be involved, as determined from the relative amounts of labelled fatty acid and lysolipid released from the positional isomers. Based on the effects of the length and position of the pyrene-labelled and unlabelled chains it is suggested that (1) the lysosomal A-type phospholipases acting on PyrnPCs recognize the carboxy-terminal part of the lipid acyl chains and (2) the relevant part of the binding site is relatively narrow. Thus phospholipids with added bulk in the corresponding region, such as those that are peroxidized and polymerized, may not be good substrates for the lysosomal phospholipases mentioned. The impaired hydrolysis of the most hydrophobic PyrnPCs indicates that lysosomal phospholipases may not be able to penetrate significantly into the substrate interphase, but upward movement of the lipid may be required for efficient hydrolysis. Finally, the rate of hydrolysis of many pyrenyl derivatives was found to be comparable to that of a natural phosphatidylcholine species, both in micelles and in lipoprotein particles, indicating that these derivatives can be used as faithful reporters of lysosomal degradation of natural lipids in vivo and in vitro.

Acylation↗

Lateral organization of liquid-crystalline cholesterol-dimyristoylphosphatidylcholine bilayers. Evidence for domains with hexagonal and centered rectangular cholesterol superlattices.

The lateral organization of fluid cholesterol-dimyristoylphosphatidylcholine (DMPC) bilayers was studied by measuring the response of fluorescent membrane probes, dipyrenylphosphatidylcholines (diPyrxPCs) or merocyanine 540, to the variation of cholesterol concentration. Parallel absorbance and light-scattering measurements were also carried out. The excimer-to-monomer ratio of diPyrxPCs displayed abrupt deviations at particular cholesterol mole fractions (CMFs). The most notable of these occurred at CMFs of 0.15, 0.33, and 0.67. Deviations were also frequently observed at CMFs of 0.12, 0.20, 0.25, and 0.40. Merocyanine 540 reproducibly reported deviations at CMFs of 0.15 and 0.33 and frequently reported values close to 0.12, 0.20, and 0.25. In absorbance (turbidity) and light scattering versus CMF plots, well-defined kinks were observed at CMFs of 0.16, 0.33, 0.52, and 0.67. The occurrence of kinks or other deviations at those particular CMFs is most readily explained in terms of a superlattice model previously developed to explain the lateral distribution of pyrenylphospholipids in bilayers [Somerharju, et al. (1985) Biochemistry 24, 2773-2781; Virtanen, J. A., et al. (1988) J. Mol. Electron. 4, 233-236]. This model is based on the assumptions that (i) each cholesterol molecule replaces a single acyl chain in a hexagonal lattice, (ii) cholesterol molecules, because of their larger size, perturb the lattice, (iii) this perturbation is minimized when the cholesterol molecules are maximally separated from each other, and (iv) the maximal separation is achieved when the cholesterol molecules form a hexagonal or centered rectangular superlattice. All detected critical CMFs, except that at CMF 0.67, are predicted by the model, thus strongly supporting its validity. The critical CMF at 0.67 is a limiting case, which can be accounted for by assuming that cholesterol and phospholipid molecules form alternating rows, i.e., formation of a cholesterol superlattice with rectangular symmetry. As predicted by the superlattice model, composition-driven order-to-disorder transitions occur between the critical CMFs, as indicated by increased data scatter and sample fluctuations in those regions. Another important prediction of the superlattice model is that domains with different cholesterol superlattices should coexist at most cholesterol concentrations. Such domains do not have to be extensive to account for the critical events observed here; rather, they are expected to be dynamic entities of limited size. It is very likely that such microscopic domains with distinct cholesterol superlattices also coexist in biological membranes. This is expected to have remarkable effects on both the structure and functions of these membranes.

Binding Sites↗

Transversal distribution of acyl-linked pyrene moieties in liquid-crystalline phosphatidylcholine bilayers. A fluorescence quenching study.

Quenching of the fluorescence of pyrene-labeled phospholipids by dibromolipids was used to determine the chain length dependence of the bilayer depths of the pyrenyl moieties. Six 1-palmitoyl-2-(pyrenyl-n-acyl)-phosphatidylcholines (PyrnPC) were examined, with end-labeled pyrenyl chains varying in length, n, from 4 to 14 carbons. These lipids were incorporated, at a concentration of 0.3 mol%, into bilayers composed of various mixtures of 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) and of one of three 1-palmitoyl-2-(x,y-dibromostearoyl)phosphatidylcholine quencher lipids (Brx,yPC; x,y = 6,7; 9,10; or 11,12). Parallel experiments were carried out with bilayers containing 50 mol % cholesterol. Quenching in these systems is dynamic, as demonstrated by the identical dependence of steady-state fluorescence intensities and excited state lifetimes of Pyr8PC on the mole fraction of Br6,7PC. Stern--Volmer analysis of the Brx,PC mole fraction dependence of PyrnPC fluorescence yielded apparent quenching constants, KSV, which show a systematic relation with both the length of the pyrenyl acyl chain and the position of the bromine atoms. The quenching data were further analyzed by plotting KSV as a function of n (defined above), or b (the average of the two bromine positions for each PyrnPC), or n--b (the separation between pyrenes and bromines). In all cases, the data were fit by Gaussian functions yielding estimates of the centers and the apparent 1/e half-widths of the transversal distributions of the pyrenyl moieties in methylene units (mu). Both in the absence and in the presence of cholesterol, the position of each PyrnPC Gaussian center is equal to the sum of n plus a constant d approximately 2.5 mu, corresponding to the distance from the effective center of the pyrenyl moiety to its point of attachment to the acyl chain.(ABSTRACT TRUNCATED AT 250 WORDS)

Crystallization↗

Dipyrenylphosphatidylcholines as membrane fluidity probes. Pressure and temperature dependence of the intramolecular excimer formation rate.

We have measured the pressure dependence of the intramolecular excimer formation rate, K(p), for di-(1'-pyrenedecanoyl)-phosphatidylcholine (dipy10PC) probes in single-component lipid multilamellar vesicles (MLV) as a function of temperature. Apparent volumes of activation (V(a)) for intramolecular excimer formation are obtained from the slopes of plots of log K(p) versus P. For liquid-crystalline saturated lipid MLV (DMPC and DPPC), these plots are linear and yield a unique V(a) at each temperature, whereas for unsaturated lipids (POPC and DOPC) they are curvilinear and V(a) appears to decrease with pressure. The isothermal pressure induced phase transition is marked by an abrupt drop in the values of K(p). The pressure to temperature equivalence values, dPm/dT, estimated from the midpoint of the transitions, are 47.0, 43.5, and 52.5 bar degree C-1 for DMPC, DPPC, and POPC, respectively. In liquid-crystalline DMPC, V(a) decreases linearly as a function of temperature, with a coefficient -dVa/dT = 0.65 +/- 0.11 ml degree C-1 mol-1. Using a modified free volume model of diffusion, we show that this value corresponds to the thermal expansivity of DMPC. Both the apparent energy and entropy of activation, Ea and delta Sa, increase with pressure in DMPC, whereas both decrease in POPC and DOPC. This difference is attributed to the sensitivity of the dynamics and/or packing of the dipy10PC probes to the location of the cis-double bonds in the chains of the unsaturated host phospholipids. Finally, the atmospheric pressure values of Ea and delta Sa for the four host MLV examined are shown to be linearly related. The relevance of this finding with respect to the structure of the excimers formed by the dipy10PC probes is briefly discussed.

Biophysical Phenomena↗

Distraction bone healing.

Bone formation by distraction was studied using three different experimental models: (1) Physeal distraction of the sheep radius was performed in 20 animals. (2) Distraction after osteotomy of the radius was carried out in 39 sheep. (3) Mandibular distraction after osteotomy was performed in 17 sheep. Formation of the organic matrix and osteogenesis were studied by radiographic, histologic, and biochemical methods as well as by electron microscopy. The mode of osteogenesis was essentially similar in all of these distraction models. Bone formation was preceded by organization of the collagenous matrix in the distraction area. In the beginning of the distraction, the gap was composed of a heterogeneous cell population, with large polymorphic fibroblast-like cells. The cells in the central part differentiated into fibroblasts, which remained functionally active as long as distraction proceeded. During physeal distraction, bone formed from the epiphyseal and metaphyseal sides as well as from the surrounding perichondrium. Also, in osteotomy distraction of both tubular bone and mandible, bone formed centripetally from the osteotomized bone ends toward the center of the gap. The organic matrix was composed almost solely of Type I collagen in the earliest stages, suggesting that the mode of osteogenesis differs from bone repair by fracture callus. The structure of the distracted segment was mainly lamellar trabecular. Corticalization of the lengthened bone segment occurred gradually after several months.

Animals↗