A novel phantom and method for comprehensive 3-dimensional measurement and correction of geometric distortion in magnetic resonance imaging

A phantom that can be used for mapping geometric distortion in magnetic resonance imaging (MRI) is described. This phantom provides an array of densely distributed control points in three-dimensional (3D) space. These points form the basis of a comprehensive measurement method to correct for geometric distortion in MR images arising principally from gradient field non-linearity and magnet field inhomogeneity. The phantom was designed based on the concept that a point in space can be defined using three orthogonal planes. This novel design approach allows for as many control points as desired. Employing this novel design, a highly accurate method has been developed that enables the positions of the control points to be measured to sub-voxel accuracy. The phantom described in this paper was constructed to fit into a body coil of a MRI scanner, (external dimensions of the phantom were: 310 mm x 310 mm x 310 mm), and it contained 10,830 control points. With this phantom, the mean errors in the measured coordinates of the control points were on the order of 0.1 mm or less, which were less than one tenth of the voxel's dimensions of the phantom image. The calculated three-dimensional distortion map, i.e., the differences between the image positions and true positions of the control points, can then be used to compensate for geometric distortion for a full image restoration. It is anticipated that this novel method will have an impact on the applicability of MRI in both clinical and research settings. especially in areas where geometric accuracy is highly required, such as in MR neuro-imaging. (C) 2004 Elsevier Inc. All rights reserved.

Analysis of cardiac ventricular wall motion based on a three-dimensional electromechanical biventricular model

This paper describes a biventricular model, which couples the electrical and mechanical properties of the heart, and computer simulations of ventricular wall motion and deformation by means of a biventricular model. In the constructed electromechanical model, the mechanical analysis was based on composite material theory and the finite-element method; the propagation of electrical excitation was simulated using an electrical heart model, and the resulting active forces were used to calculate ventricular wall motion. Regional deformation and Lagrangian strain tensors were calculated during the systole phase. Displacements, minimum principal strains and torsion angle were used to describe the motion of the two ventricles. The simulations showed that during the period of systole, (1) the right ventricular free wall moves towards the septum, and at the same time, the base and middle of the free wall move towards the apex, which reduces the volume of the right ventricle; the minimum principle strain (E3) is largest at the apex, then at the middle of the free wall and its direction is in the approximate direction of the epicardial muscle fibres; (2) the base and middle of the left ventricular free wall move towards the apex and the apex remains almost static; the torsion angle is largest at the apex; the minimum principle strain E3 is largest at the apex and its direction on the surface of the middle wall of the left ventricle is roughly in the fibre orientation. These results are in good accordance with results obtained from MR tagging images reported in the literature. This study suggests that such an electromechanical biventricular model has the potential to be used to assess the mechanical function of the two ventricles, and also could improve the accuracy ECG simulation when it is used in heart torso model-based body surface potential simulation studies.

CFD simulations of hydrocyclones with an air core - Comparison between large eddy simulations and a second moment closure

CFD simulations of the 75 mm, hydrocyclone of Hsieh (1988) have been conducted using Fluent TM. The simulations used 3-dimensional body fitted grids. The simulations were two phase simulations where the air core was resolved using the mixture (Manninen et al., 1996) and VOF (Hirt and Nichols, 1981) models. Velocity predictions from large eddy simulations (LES), using the Smagorinsky-Lilly sub grid scale model (Smagorinsky, 1963; Lilly, 1966) and RANS simulations using the differential Reynolds stress turbulence model (Launder et al., 1975) were compared with Hsieh's experimental velocity data. The LES simulations gave very good agreement with Hsieh's data but required very fine grids to predict the velocities correctly in the bottom of the apex. The DRSM/RANS simulations under predicted tangential velocities, and there was little difference between the velocity predictions using the linear (Launder, 1989) and quadratic (Speziale et al., 1991) pressure strain models. Velocity predictions using the DRSM turbulence model and the linear pressure strain model could be improved by adjusting the pressure strain model constants.

A proposed scheme for comprehensive characterization of the measured geometric distortion in magnetic resonance imaging using a three-dimensional phantom

Recently, a 3-dimensional phantom that can provide a comprehensive, accurate and complete measurement of the geometric distortion in MRI has been developed. In this paper, a scheme for characterizing the measured geometric distortion using the 3-D phantom is described. In the proposed scheme, a number of quantitative measures are developed and used to characterize the geometric distortion. These measures encompass the overall and spatial aspects of the geometric distortion. Two specific types of volume of interest, rectangular parallelepipeds (including cubes) and spheres are considered in the proposed scheme. As an illustration, characterization of the geometric distortion in a Siemens 1.5T Sonata MRI system using the proposed scheme is presented. As shown, the proposed scheme provides a comprehensive assessment of the geometric distortion. The scheme can be potentially used as a standard procedure for the assessment of geometric distortion in MRI. (C) 2004 American Association of Physicists in Medicine.

Use of real-time three-dimensional echocardiography to measure left atrial volume: Comparison with other echocardiographic techniques

Objectives: Left atrial (LA) volume (LAV) is a prognostically important biomarker for diastolic dysfunction, but its reproducibility on repeated testing is not well defined. LA assessment with 3-dimensional. (3D) echocardiography (3DE) has been validated against magnetic resonance imaging, and we sought to assess whether this was superior to existing measurements for sequential echocardiographic follow-up. Methods: Patients (n = 100; 81 men; age 56 +/- 14 years) presenting for LA evaluation were studied with M-mode (MM) echocardiography, 2-dimensional (2D) echocardiography, and 3DE. Test-retest variation was performed by a complete restudy by a separate sonographer within 1 hour without alteration of hemodynamics or therapy. In all, 20 patients were studied for interobserver and intraobserver variation. LAVs were calculated by using M-mode diameter and planimetered atrial area in the apical. 4-chamber view to calculate an assumed sphere, as were prolate ellipsoid, Simpson's biplane, and biplane area-length methods. All were compared with 3DE. Results: The average LAV was 72 +/- 27 mL by 3DE. There was significant underestimation of LAV by M-mode (35 +/- 20 mL, r = 0.66, P < .01). The 3DE and various 2D echocardiographic techniques were well correlated: LA planimetry (85 +/- 38 mL, r = 0.77, P < .01), prolate ellipsoid (73 +/- 36 mL, r = 0.73, P = .04), area-length (64 +/- 30 mL, r = 0.74, P < .01), and Simpson's biplane (69 +/- 31 mL, r = 0.78, P = .06). Test-retest variation for 3DE was most favorable (r = 0.98, P < .01), with the prolate ellipsoid method showing most variation. Interobserver agreement between measurements was best for 3DE (r = 0.99, P < .01), with M-mode the worst (r = 0.89, P < .01). Intraobserver results were similar to interobserver, the best correlation for 3DE (r = 0.99, P < .01), with LA planimetry the worst (r = 0.91, P < .01). Conclusions. The 2D measurements correlate closely with 3DE. Follow-up assessment in daily practice appears feasible and reliable with both 2D and 3D approaches.

Anisotropic viscous models of large-deformation Mohr-Coulomb failure

We have developed a way to represent Mohr-Coulomb failure within a mantle-convection fluid dynamics code. We use a viscous model of deformation with an orthotropic viscoplasticity (a different viscosity is used for pure shear to that used for simple shear) to define a prefered plane for slip to occur given the local stress field. The simple-shear viscosity and the deformation can then be iterated to ensure that the yield criterion is always satisfied. We again assume the Boussinesq approximation, neglecting any effect of dilatancy on the stress field. An additional criterion is required to ensure that deformation occurs along the plane aligned with maximum shear strain-rate rather than the perpendicular plane, which is formally equivalent in any symmetric formulation. We also allow for strain-weakening of the material. The material can remember both the accumulated failure history and the direction of failure. We have included this capacity in a Lagrangian-integration-point finite element code and show a number of examples of extension and compression of a crustal block with a Mohr-Coulomb failure criterion. The formulation itself is general and applies to 2- and 3-dimensional problems.

Method for generation of homogeneous multicellular tumor spheroids applicable to a wide variety of cell types

Multicellular tumor spheroids (MCTS) are used as organotypic models of normal and solid tumor tissue. Traditional techniques for generating MCTS, such as growth on nonadherent surfaces, in suspension, or on scaffolds, have a number of drawbacks, including the need for manual selection to achieve a homogeneous population and the use of nonphysiological matrix compounds. In this study we describe a mild method for the generation of MCTS, in which individual spheroids form in hanging drops suspended from a microtiter plate. The method has been successfully applied to a broad range of cell lines and shows nearly 100% efficiency (i.e., one spheroid per drop). Using the hepatoma cell line, HepG2, the hanging drop method generated well-rounded MCTS with a narrow size distribution (coefficient of variation [CV] 10% to 15%, compared with 40% to 60% for growth on nonadherent surfaces). Structural analysis of HepG2 and a mammary gland adenocarcinoma cell line, MCF-7, composed spheroids, revealed highly organized, three-dimensional, tissue-like structures with an extensive extracellular matrix. The hanging drop method represents an attractive alternative for MCTS production, because it is mild, can be applied to a wide variety of cell lines, and can produce spheroids of a homogeneous size without the need for sieving or manual selection. The method has applications for basic studies of physiology and metabolism, tumor biology, toxicology, cellular organization, and the development of bioartificial tissue. (C) 2003 Wiley Periodicals, Inc.

Structure of Petunia hybrida Defensin 1, a Novel Plant Defensin with Five Disulfide Bonds

The structure of a novel plant defensin isolated from the flowers of Petunia hybrida has been determined by H-1 NMR spectroscopy. P. hybrida defensin 1 (PhD1) is a basic, cysteine-rich, antifungal protein of 47 residues and is the first example of a new subclass of plant defensins with five disulfide bonds whose structure has been determined. PhD1 has the fold of the cysteine-stabilized alphabeta motif, consisting of an alpha-helix and a triple-stranded antiparallel beta-sheet, except that it contains a fifth disulfide bond from the first loop to the alpha-helix. The additional disulfide bond is accommodated in PhD1 without any alteration of its tertiary structure with respect to other plant defensins. Comparison of its structure with those of classic, four-disulfide defensins has allowed us to identify a previously unrecognized hydrogen bond network that is integral to structure stabilization in the family.

Microcin J25 has a threaded sidechain-to-backbone ring structure and not a head-to-tail cyclized backbone

Microcin J25 is a 21 amino acid bacterial peptide that has potent antibacterial activity against Gram-negative bacteria, resulting from its interaction with RNA polymerase. The peptide was previously proposed to have a head-to-tail cyclized peptide backbone and a tight globular structure (Blond, A., Peduzzi, J., Goulard, C., Chiuchiolo, M. J., Barthelemy, M., Prigent, Y., Salomon, R. A., Farias, R. N., Moreno, F. & Rebuffat, S. Eur. J. Biochem. 1999, 259, 747-755). It exhibits remarkable thermal stability for a peptide of its size lacking disulfide bonds and in part this was previously proposed to derive from its macrocyclic structure. We show here that in fact the peptide does not have a head-to-tail cyclic structure but rather a side chain to backbone cyclization between Glu8 and the N-terminus. This creates an embedded ring that is threaded by the C-terminal tail of the molecule, forming a noose-like feature. The three-dimensional structure deduced from NMR data suggests that slippage of the noose is prevented by two aromatic residues flanking the embedded ring. Unthreading does not occur even when the molecule is enzymatically digested with thermolysin. The new structural interpretation fully accounts for previously reported NMR and biophysical data and is consistent with the remarkable stability of this potent antimicrobial peptide.

Twists, knots, and rings in proteins: Structural definition of the cyclotide framework

In recent years an increasing number of miniproteins containing an amide-cyclized backbone have been discovered. The cyclotide family is the largest group of such proteins and is characterized by a circular protein backbone and six conserved cysteine residues linked by disulfide bonds in a tight core of the molecule. These form a cystine knot in which an embedded ring formed by two of the disulfide bonds and the connecting backbone segment is threaded by a third disulfide bond. In the current study we have undertaken high resolution structural analysis of two prototypic cyclotides, kalata B1 and cycloviolacin O1, to define the role of the conserved residues in the sequence. We provide the first comprehensive analysis of the topological features in this unique family of proteins, namely rings (a circular backbone), twists (a cis-peptide bond in the Mobius cyclotides) and knots (a knotted arrangement of the disulfide bonds).

Does interhemispheric competition mediate motion-induced blindness? A transcranial magnetic stimulation study

Motion-induced blindness (MIB) is a phenomenon, perhaps related to perceptual rivalry, where stationary targets disappear and reappear in a cyclic mode when viewed against a background (mask) of coherent, apparent 3-D motion. Since MIB has recently been shown to share similar temporal properties with binocular rivalry, we probed the appearance-disappearance cycle of MIB using unilateral, single-pulse transcranial magnetic stimulation (TMS)-a manipulation that has previously been shown to influence binocular rivalry. Effects were seen for both hemispheres when the timing of TMS was determined prospectively on the basis of a given subject's appearance-disappearance cycle, so that it occurred on average around 300 ms before the time of perceptual switch. Magnetic stimulation of either hemisphere shortened the time to switch from appearance to disappearance and vice versa. However, TMS of left posterior parietal cortex more selectively shortened the disappearance time of the targets if delivered in phase with the disappearance cycle, but lengthened it if TMS was delivered in the appearance phase after the perceptual switch. Opposite effects were seen in the right hemisphere, although less marked than the left-hemisphere effects. As well as sharing temporal characteristics with binocular rivalry, MIB therefore seems to share a similar underlying mechanism of interhemispheric modulation. Interhemispheric switching may thus provide a common temporal framework for uniting the diverse, multilevel phenomena of perceptual rivalry.

Comparative habitat use in Juniper woodland bird community

We compared vegetation structure used by 14 bird species during the 1998 and 1999 breeding seasons to determine what habitat features best accounted for habitat division and community organization in Utah juniper (Juniperus osteosperma) woodlands of southwestern Wyoming. Habitat use was quantified by measuring 24 habitat variables in 461 bird-centered quadrats, each 0.04 ha in size. Using discriminant function analysis, we differentiated between habitat used by 14 bird species along 3 habitat dimensions: (1) variation in shrub cover, overstory juniper cover, mature tree density, understory height, and decadent tree density; (2) a gradient composed of elevation and forb cover; and (3) variation in grass cover, tree height, seedling/sapling cover, and bare ground/rock cover. Of 14 species considered, 9 exhibited substantial habitat partitioning: Mourning Dove (Zenaida macroura), Bewick's Wren (Thryomanes bewickii), Blue-gray Gnatcatcher (Polioptila caerulea), Mountain Bluebird (Sialia currucoides), Plumbeous Vireo (Vireo plumbeus), Green-tailed Towhee (Pipilo chlorurus), Brewer's Sparrow (Spizella breweri), Dark-eyed Junco (Junco hyemalis), and Cassin's Finch (Carpodacus cassinii). Our results indicate juniper bird communities of southwestern Wyoming are organized along a 3-dimensional habitat gradient composed of woodland maturity, elevation, and juniper recruitment. Because juniper birds partition habitat along successional and altitudinal gradients, indiscriminate woodland clearing as well as continued fire suppression will alter species composition. Restoration efforts should ensure that all successional stages of juniper woodland are present on the landscape.