Ellipsoidal harmonic (Lame) MRI shims

Ellipsoidal harmonics are presented as a basis function set for the design of shim coils for magnetic resonance imaging (MRI) or spectroscopy. MR shim coils may be either superconductive or resistive. Ellipsoidal harmonics form an orthogonal set over an ellipsoid and hence are appropriate in circumstances where the imaging or spectroscopic region of a magnet more closely conforms to an ellipsoid rather than a sphere. This is often the case in practice. The Cartesian form of ellipsoidal harmonics is discussed. A method for the design of streamline coil designs is detailed and patterns for third-order ellipsoidal (Lame) shims wound on a cylindrical surface are presented.

Capillary network models for transport in packed beds: Considerations of pore aspect ratio

The conventional analysis for the estimation of the tortuosity factor for transport in porous media is modified here to account for the effect of pore aspect ratio. Structural models of the porous medium are also constructed for calculating the aspect ratio as a function of porosity. Comparison of the model predictions with the extensive data of Currie (1960) for the effective diffusivity of hydrogen in packed beds shows good agreement with a network model of randomly oriented intersecting pores for porosities upto about 50 percent, which is the region of practical interest. The predictions based on this network model are also found to be in better agreement with the data of Currie than earlier expressions developed for unconsolidated and grainy media.

A method for the design of MRI radiofrequency coils based on triangular and pulse basis functions

This paper presents a numerical technique for the design of an RF coil for asymmetric magnetic resonance imaging (MRI) systems. The formulation is based on an inverse approach where the cylindrical surface currents are expressed in terms of a combination of sub-domain basis functions: triangular and pulse functions. With the homogeneous transverse magnetic field specified in a spherical region, a functional method is applied to obtain the unknown current coefficients. The current distribution is then transformed to a conductor pattern by use of a stream function technique. Preliminary MR images acquired using a prototype RF coil are presented and validate the design method. (C) 2002 Elsevier Science B.V. All rights reserved.

Alterations in cardiac dihydropyridine receptors and calcium channel function in mdx mice

Duchenne muscular dystrophy (DMD) is a fatal neuromuscular condition affecting approximately one in 3500 live male births resulting from the lack of the myocyte protein dystrophin. The absence of dystrophin in cardiac myocytes is associated with calcium overload which in turn activates calcium-dependent proteolytic enzymes contributing to congestive heart failure, muscle necrosis and fibrosis. To date, the basis for the calcium overload has not been determined. Since L-type calcium channels are a major mediator of calcium influx we determined their potential contribution to the calcium overload. Male muscular dystrophy (mdx) mice and control C57BL10ScSn (C57) mice aged 12– 16 weeks were used in all experiments. In tissue bath studies, isolated contracting left atria from mdx revealed a reduced potency to the dihydropyridine (DHP) agonist BayK8644 and antagonist nifedipine (P < 0.05). Similarly, radioligand binding studies using the DHP antagonist [3H]-PN 200-110 showed a reduced potency (P < 0.05) in isolated membranes, associated with an increased receptor density (P < 0.05). The increased receptor density was supported by RT-PCR experiments revealing increased RNAfor the DHP receptor. Patch clamp studies revealed the presence of a diltiazem sensitive calcium current that showed delayed inactivation in isolated mdx myocytes (P < 0.01). In conclusion, the increased number of DHP binding sites and the delay in L-type current inactivation may both contribute to increased calcium influx and hence calcium overload in the dystrophin deficient mdx cardiac myocytes.

Mathematically gifted male adolescents activate a unique brain network during mental rotation

Mental rotation involves the creation and manipulation of internal images, with the later being particularly useful cognitive capacities when applied to high-level mathematical thinking and reasoning. Many neuroimaging studies have demonstrated mental rotation to be mediated primarily by the parietal lobes, particularly on the right side. Here, we use fMRI to show for the first time that when performing 3-dimensional mental rotations, mathematically gifted male adolescents engage a qualitatively different brain network than those of average math ability, one that involves bilateral activation of the parietal lobes and frontal cortex, along with heightened activation of the anterior cingulate. Reliance on the processing characteristics of this uniquely bilateral system and the interplay of these anterior/posterior regions may be contributors to their mathematical precocity.

Functional MR Imaging of visuo-spatial network in Autism Spectrum Disorder

Individuals with Autism Spectrum Disorder (ASD) are generally thought to have impaired attentional and executive function upon which all their cognitive and behaviour functions are based. Mental Rotation is a recognized visuo-spatial task, involving spatial working memory, known to involve activation in the fronto-parietal networks. To elucidate the functioning of fronto-parietal networks in ASD, the aim of this study was to use fMRI techniques with a mental rotation task, to characterize the underlying functional neural system. Sixteen male participants (seven highfunctioning autism or Asperger's syndrome; nine ageand performance IQ-matched controls) underwent fMRI. Participants were presented with 18 baseline and 18 rotation trials, with stimuli rotated 3- dimensionaUy (45°-180°). Data were acquired on a 3- Tesla scanner. The most widely accepted area reported to be involved in processing of visuo-spatial information. Posterior Parietal Cortex, was found to be activated in both groups, however, the ASD group showed decreased activation in cortical and subcortical frontal structures that are highly interconnected, including lateral and medial Brodmann area 6, frontal eye fields, caudate, dorsolateral prefrontal cortex and anterior cingulate. The suggested connectivity between these regions indicates that one or more circuits are impaired as a result of the disorder. In future it is hoped that we are able to identify the possible point of origin of this dysfunction, or indeed if the entire network is dysfunctional.

The dorsomedial visual areas in New World and Old World monkeys: homology and function

The extrastriate cortex near the dorsal midline has been described as part of an 'express' pathway that provides visual input to the premotor cortex. This pathway is considered important for the integration of sensory information about the visual field periphery and the skeletomotor system, especially in relation to the control of arm movements. However, a better understanding of the functional contributions of different parts of this complex has been hampered by the lack of data on the extent and boundaries of its constituent visual areas. Recent studies in macaques have provided the first detailed view of the topographical organization of this region in Old World monkeys. Despite differences in nomenclature, a comparison of the visuotopic organization, myeloarchitecture and connections of the relevant visual areas with those previously studied in New World monkeys reveals a remarkable degree of similarity and helps to clarify the subdivision of function between different areas of the dorsomedial complex. A caudal visual area, named DM or V6, appears to be important for the detection of coherent patterns of movement across wide regions of the visual field, such as those induced during self-motion. A rostral area, named M or V6A, is more directly involved with visuomotor integration. This area receives projections both from DM/V6 and from a separate motion analysis channel, centred on the middle temporal visual area (or V5), which detects the movement of objects in extrapersonal space. These results support the suggestion, made earlier on the basis of more fragmentary evidence, that the areas rostral to the second visual area in dorsal cortex are homologous in all simian primates. Moreover, they emphasize the importance of determining the anatomical organization of the cortex as a prerequisite for elucidating the function of different cortical areas.

Effect of the toxic milk mutation (tx) on the function and intracellular localization of Wnd, the murine homologue of the Wilson copperATPase

Wilson disease is an autosomal recessive copper transport disorder resulting from defective biliary excretion of copper and subsequent hepatic copper accumulation and liver failure if not treated. The disease is caused by mutations in the ATP7B (WND) gene, which is expressed predominantly in the liver and encodes a copper-transporting P-type ATPase that is structurally and functionally similar to the Menkes protein (MNK), which is defective in the X-linked copper transport disorder Menkes disease. The toxic milk (tx) mouse has a clinical phenotype similar to Wilson disease patients and, recently, the tx mutation within the murine WND homologue (Wnd) of this mouse was identified, establishing it as an animal model for Wilson disease. In this study, cDNA constructs encoding the wild-type (Wnd-wt) and mutant (Wnd-tx) Wilson proteins (Wnd) were generated and expressed in Chinese hamster ovary (CHO) cells. The fx mutation disrupted the copper-induced relocalization of Wnd in CHO cells and abrogated Wnd-mediated copper resistance of transfected CHO cells. In addition, co-localization experiments demonstrated that while Wnd and MNK are located in the trans-Golgi network in basal copper conditions, with elevated copper, these proteins are sorted to different destinations within the same cell, Ultrastructural studies showed that with elevated copper levels, Wnd accumulated in large multivesicular structures resembling late endosomes that may represent a novel compartment for copper transport. The data presented provide further support for a relationship between copper transport activity and the copper-induced relocalization response of mammalian copper ATPases, and an explanation at a molecular level for the observed phenotype of fx mice.

Efficient simulation of a Tandem Jackson Network

The two-node tandem Jackson network serves as a convenient reference model for the analysis and testing of different methodologies and techniques in rare event simulation. In this paper we consider a new approach to efficiently estimate the probability that the content of the second buffer exceeds some high level L before it becomes empty, starting from a given state. The approach is based on a Markov additive process representation of the buffer processes, leading to an exponential change of measure to be used in an importance sampling procedure. Unlike changes of measures proposed and studied in recent literature, the one derived here is a function of the content of the first buffer. We prove that when the first buffer is finite, this method yields asymptotically efficient simulation for any set of arrival and service rates. In fact, the relative error is bounded independent of the level L; a new result which is not established for any other known method. When the first buffer is infinite, we propose a natural extension of the exponential change of measure for the finite buffer case. In this case, the relative error is shown to be bounded (independent of L) only when the second server is the bottleneck; a result which is known to hold for some other methods derived through large deviations analysis. When the first server is the bottleneck, experimental results using our method seem to suggest that the relative error is bounded linearly in L.

A monolithic smoothing-gap algorithm for contact-impact based on the signed distance function

A new algorithm has been developed for smoothing the surfaces in finite element formulations of contact-impact. A key feature of this method is that the smoothing is done implicitly by constructing smooth signed distance functions for the bodies. These functions are then employed for the computation of the gap and other variables needed for implementation of contact-impact. The smoothed signed distance functions are constructed by a moving least-squares approximation with a polynomial basis. Results show that when nodes are placed on a surface, the surface can be reproduced with an error of about one per cent or less with either a quadratic or a linear basis. With a quadratic basis, the method exactly reproduces a circle or a sphere even for coarse meshes. Results are presented for contact problems involving the contact of circular bodies. Copyright (C) 2002 John Wiley Sons, Ltd.

Characterization of the retinoid orphan-related receptor-alpha coactivator binding interface: A structural basis for ligand-independent transcription

The retinoid orphan-related receptor-alpha (RORalpha) is a member of the ROR subfamily of orphan receptors and acts as a constitutive activator of transcription in the absence of exogenous ligands. To understand the basis of this activity, we constructed a homology model of Rill using the closely related TRbeta as a template. Molecular modeling suggested that bulky hydrophobic side chains occupy the RORa ligand cavity leaving a small but distinct cavity that may be involved in receptor stabilization. This model was subject to docking simulation with a receptor-interacting peptide from the steroid receptor coactivator, GR-interacting protein-1, which delineated a coactivator binding surface consisting of the signature motif spanning helices 3-5 and helix 12 [activation function 2 (AF2)]. Probing this surface with scanning alanine mutagenesis showed structural and functional equivalence between homologous residues of RORalpha and TRbeta. This was surprising (given that Rill is a ligand-independent activator, whereas TRbeta has an absolute requirement for ligand) and prompted us to use molecular modeling to identify differences between Rill and TRbeta in the way that the All helix interacts with the rest of the receptor. Modeling highlighted a nonconserved amino acid in helix 11 of RORa (Phe491) and a short-length of 3.10 helix at the N terminus of AF2 which we suggest i) ensures that AF2 is locked permanently in the holoconformation described for other liganded receptors and thus 2) enables ligand-independent recruitment of coactivators. Consistent with this, mutation of RORa Phe491 to either methionine or alanine (methionine is the homologous residue in TRbeta), reduced and ablated transcriptional activation and recruitment of coactivators, respectively. Furthermore, we were able to reconstitute transcriptional activity for both a deletion mutant of Ill lacking All and Phe491 Met, by overexpression of a GAL-AF2 fusion protein, demonstrating ligand-independent recruitment of AF2 and a role for Phe491 in recruiting AF2.

Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

An augmented Lanczos algorithm for the efficient computation of a dot-product of a function of a large sparse symmetric matrix

The Lanczos algorithm is appreciated in many situations due to its speed. and economy of storage. However, the advantage that the Lanczos basis vectors need not be kept is lost when the algorithm is used to compute the action of a matrix function on a vector. Either the basis vectors need to be kept, or the Lanczos process needs to be applied twice. In this study we describe an augmented Lanczos algorithm to compute a dot product relative to a function of a large sparse symmetric matrix, without keeping the basis vectors.

Renal function and cardiovascular risk markers in a remote Australian Aboriginal community

Background. Australian Aborigines living in remote areas have exceedingly high rates of renal failure together with increased cardiovascular morbidity and mortality. To examine the basis of this association, we studied markers of renal function and cardiovascular (CV) risk in a coastal Aboriginal community in a remote area of the Northern Territory of Australia. End-stage renal disease (ESRD) incidence rates in that community are 15 times the national non-Aboriginal rate and CV mortality rates in the region are increased 5-fold. Methods. A cross-sectional community survey was conducted. Markers of early renal disease examined included urine albumin/creatinine ratio (ACR), serum creatinine concentration and calculated glomerular filtration rate (GFR). CV risk markers included blood pressure as well as measures of glycaemia, diabetes and serum lipids. Results. The study group included 237 people, 58% of the adult population of the community. The crude prevalence of microalbuminuria (urine ACR: 3.4-33.9 g/mol, 30-299 mg/g) was 31% and of overt albuminuria (urine ACR: greater than or equal to34 g/mol, greater than or equal to300 mg/g), 13%. The prevalence of overt albuminuria increased with age, but the prevalence of microalbuminuria was greatest in the 45-54 year age group. Microalbuminuria was associated with increasing body mass index, whereas overt albuminuria was associated with increasing glycated haemoglobin (HbA1c) and systolic blood pressure and a history of diabetes. The prevalence of elevated serum creatinine concentration (greater than or equal to120 mumol/l) was 10%. GFR (calculated using the MDRD equation) was <60 ml/min/1.73m(2) in 12% and 60-79 ml/min/1.73 m(2) in a further 36% of the study population. Although many people with albuminuria had well preserved GFRs, mean GFR was lower in people with higher levels of albuminuria. Conclusions. The high prevalence of markers of renal disease in this community was consistent with their high rates of ESRD. The distribution of microalbuminuria suggested a 'cohort effect', representing a group who will progress to overt albuminuria. The powerful association of renal disease markers with CV risk factors confirms a strong link between renal and CV disease in the early, asymptomatic stages of each. Thus, pathologic albuminuria, in part, might be a manifestation of the metabolic/haemodynamic syndrome and both conditions might arise out of a common menu of risk factors. Hence, a single agenda of primary and secondary intervention may benefit both.

Typing mineral deposits using their grades and tonnages in an artificial neural network

A test of the ability of a probabilistic neural network to classify deposits into types on the basis of deposit tonnage and average Cu, Mo, Ag, Au, Zn, and Pb grades is conducted. The purpose is to examine whether this type of system might serve as a basis for integrating geoscience information available in large mineral databases to classify sites by deposit type. Benefits of proper classification of many sites in large regions are relatively rapid identification of terranes permissive for deposit types and recognition of specific sites perhaps worthy of exploring further. Total tonnages and average grades of 1,137 well-explored deposits identified in published grade and tonnage models representing 13 deposit types were used to train and test the network. Tonnages were transformed by logarithms and grades by square roots to reduce effects of skewness. All values were scaled by subtracting the variable's mean and dividing by its standard deviation. Half of the deposits were selected randomly to be used in training the probabilistic neural network and the other half were used for independent testing. Tests were performed with a probabilistic neural network employing a Gaussian kernel and separate sigma weights for each class (type) and each variable (grade or tonnage). Deposit types were selected to challenge the neural network. For many types, tonnages or average grades are significantly different from other types, but individual deposits may plot in the grade and tonnage space of more than one type. Porphyry Cu, porphyry Cu-Au, and porphyry Cu-Mo types have similar tonnages and relatively small differences in grades. Redbed Cu deposits typically have tonnages that could be confused with porphyry Cu deposits, also contain Cu and, in some situations, Ag. Cyprus and kuroko massive sulfide types have about the same tonnages. Cu, Zn, Ag, and Au grades. Polymetallic vein, sedimentary exhalative Zn-Pb, and Zn-Pb skarn types contain many of the same metals. Sediment-hosted Au, Comstock Au-Ag, and low-sulfide Au-quartz vein types are principally Au deposits with differing amounts of Ag. Given the intent to test the neural network under the most difficult conditions, an overall 75% agreement between the experts and the neural network is considered excellent. Among the largestclassification errors are skarn Zn-Pb and Cyprus massive sulfide deposits classed by the neuralnetwork as kuroko massive sulfides—24 and 63% error respectively. Other large errors are the classification of 92% of porphyry Cu-Mo as porphyry Cu deposits. Most of the larger classification errors involve 25 or fewer training deposits, suggesting that some errors might be the result of small sample size. About 91% of the gold deposit types were classed properly and 98% of porphyry Cu deposits were classes as some type of porphyry Cu deposit. An experienced economic geologist would not make many of the classification errors that were made by the neural network because the geologic settings of deposits would be used to reduce errors. In a separate test, the probabilistic neural network correctly classed 93% of 336 deposits in eight deposit types when trained with presence or absence of 58 minerals and six generalized rock types. The overall success rate of the probabilistic neural network when trained on tonnage and average grades would probably be more than 90% with additional information on the presence of a few rock types.