Numerical modeling of 11.1 T MRI of a human head using MoM/FDTD methods

In a recent study, severe distortions in the proton images of an excised, fixed, human brain in an 11.1 Tesla/40 cm MR instrument have been observed, and the effect modeled on phantom images using a finite difference time domain (FDTD) model. in the present study, we extend these simulations to that of a complete human head, employing a hybrid FDTD and method of moments (MoM) approach, which provides a validated method for simulating biological samples in coil structures. The effect of fixative on the image distortions is explored. importantly, temperature distributions within the head are also simulated using a bioheat method based on parameters derived from the electromagnetic simulations. The MoM/FDTD simulations confirm that the transverse magnetic field (B,) from a ReCav resonator exhibits good homogeneity in air but strong inhomogeneity when loaded with the head with or without fixative. The fixative serves to increase the distortions, but they are still significant for the in vivo simulations. The simulated signal intensity (SI) distribution within the sample confirm the distortions in the experimental images are caused by the complex interactions of the incident electromagnetic fields with tissue, which is heterogeneous in terms of conductivity and permittivity. The temperature distribution is likewise heterogeneous, raising concerns regarding hot spot generation in the sample that may exceed acceptable levels in future in vivo studies. As human imaging at 11.1 T is some time away, simulations are important in terms of predicting potential safety issues as well as evaluating practical concerns about the quality of images. Simulation on a whole human head at 11.1 T implies the wave behavior presents significant engineering challenges for ultra-high-field (UHF) MRI. Novel strategies will have to be employed in imaging technique and resonator design for UHF MRI to achieve the theoretical signal-to-noise ratio (SNR) improvements it offers over lower field systems. (C) 2005 Wiley Periodicals, Inc.

Adsorption in slit pores and pore-size distribution: A molecular layer structure theory

A new approach is developed to analyze the thermodynamic properties of a sub-critical fluid adsorbed in a slit pore of activated carbon. The approach is based on a representation that an adsorbed fluid forms an ordered structure close to a smoothed solid surface. This ordered structure is modelled as a collection of parallel molecular layers. Such a structure allows us to express the Helmholtz free energy of a molecular layer as the sum of the intrinsic Helmholtz free energy specific to that layer and the potential energy of interaction of that layer with all other layers and the solid surface. The intrinsic Helmholtz free energy of a molecular layer is a function (at given temperature) of its two-dimensional density and it can be readily obtained from bulk-phase properties, while the interlayer potential energy interaction is determined by using the 10-4 Lennard-Jones potential. The positions of all layers close to the graphite surface or in a slit pore are considered to correspond to the minimum of the potential energy of the system. This model has led to accurate predictions of nitrogen and argon adsorption on carbon black at their normal boiling points. In the case of adsorption in slit pores, local isotherms are determined from the minimization of the grand potential. The model provides a reasonable description of the 0-1 monolayer transition, phase transition and packing effect. The adsorption of nitrogen at 77.35 K and argon at 87.29 K on activated carbons is analyzed to illustrate the potential of this theory, and the derived pore-size distribution is compared favourably with that obtained by the Density Functional Theory (DFT). The model is less time-consuming than methods such as the DFT and Monte-Carlo simulation, and most importantly it can be readily extended to the adsorption of mixtures and capillary condensation phenomena.

Comparison of selected methods for measuring the virulence properties of Listeria spp.

The comparative ability of different methods to assess virulence of Listeria species was investigated in ten Listeria strains. All strains were initially subjected to pulsed-field gel electrophoresis analysis to determine their relatedness. Virulence characteristics were subsequently tested for by (i) determining the presence of six virulence genes by polymerase chain reaction; (ii) testing for the production of listeriolysin O, phosphatidylcholine phospholipase C, and phosphatidylinositol-specific phospholipase C; (iii) investigating the hydrophobicity of the strains; (iv) determining the strains ability to attach to, enter, and replicate within the Caco-2 cells. Variations in most of the virulence characteristics were obvious across the strains for the range of tests performed. A wide range of anomalous results among methods were apparent. In particular, the presence of virulence genes was found to be unrelated to the production of virulence-associated proteins in vitro, while virulence protein production and hydrophobicity in Listeria monocytogenes were found to be unrelated or marginally related, respectively, to the ability to invade the Caco-2 cell line. It was concluded that the methods investigated were unable to consistently and unequivocally measure the differences in the virulence properties of the strains.

Do different surrogate methods detect lateral genetic transfer events of different relative ages?

Non-tree-based ('surrogate') methods have been used to identify instances of lateral genetic transfer in microbial genomes but agreement among predictions of different methods can be poor. It has been proposed that this disagreement arises because different surrogate methods are biased towards the detection of certain types of transfer events. This conjecture is supported by a rigorous phylogenetic analysis of 3776 proteins in Escherichia coli K12 MG1655 to map the ages of transfer events relative to one another.

Molecular characterization of Thelastomatoidea (Nematoda : Oxyurida) from cockroaches in Australia

A molecular approach was used to genetically characterize 5 species (Aoruroides queenslandensis. Blattophila sphaerolaima, Cordonicola gibsoni, Desmicola ornato and Leidynemella fusiformis) belonging to the superfamily. Thelastomatoidea fi (Nematoda: Oxyurida), a group of pinworms that parasitizes terrestrial arthropods. The D3 domain of the large subunit Of nuclear ribosomal RNA (LSU) was sequenced for individual specimens, and the analysis of the sequence data allowed the genetic relationships of the 5 species to be studied dagger. The sequence variation in the D3 domain within individual species (0-1-8%) was significantly less than the differences among species (4(.)3-12(.)4%). Phylogenetic analyses, Using maximum parsimony, maximum likelihood, and neighbour-joining, tree-building methods, established relationships among the 5 species of Thelastomatoidea and Oxyuris equi (a species of the order Oxyurida). The molecular approach employed provides the prospect for developing DNA tools for the specific identification of the Thelastomatoidea, irrespective of developmental stage and sex, as a basis for systematic, ecological and/or population genetic investigations of members within this superfamily.

Molecular diversity in venom from the Australian Brown Snake, Pseudonaja textilis

Venom from the Australian elapid Pseudonaja textilis (Common or Eastern Brown snake), is the second most toxic snake venom known and is the most common cause of death from snake bite in Australia. This venom is known to contain a prothrombin activator complex, serine proteinase inhibitors, various phospholipase A(2)s, and pre-and postsynaptic neurotoxins. In this study, we performed a proteomic identification of the venom using two- dimensional gel electrophoresis, mass spectrometry, and de novo peptide sequencing. We identified most of the venom proteins including proteins previously not known to be present in the venom. In addition, we used immunoblotting and post-translational modification-specific enzyme stains and antibodies that reveal the complexity and regional diversity of the venom. Modifications observed include phosphorylation, gamma-carboxylation, and glycosylation. Glycoproteins were further characterized by enzymatic deglycosylation and by lectin binding specificity. The venom contains an abundance of glycoproteins with N-linked sugars that include glucose/mannose, N-acetylgalactosamine, N-acetylglucosamine, and sialic acids. Additionally there are multiple isoforms of mammalian coagulation factors that comprise a significant proportion of the venom. Indeed two of the identified proteins, a procoagulant and a plasmin inhibitor, are currently in development as human therapeutic agents.

Molecular dynamics characterization of n-octyl-beta-D-glucopyranoside micelle structure in aqueous solution

n-Octyl-beta-D-glueopyranoside (OG) is a non-ionic glycolipid, which is used widely in biotechnical and biochemical applications. All-atom molecular dynamics simulations from two different initial coordinates and velocities in explicit solvent have been performed to characterize the structural behaviour of an OG aggregate at equilibrium conditions. Geometric packing properties determined from the simulations and small angle neutron scattering experiment state that OG micelles are more likely to exist in a non-spherical shape, even at the concentration range near to the critical micelle concentration (0.025 M). Despite few large deviations in the principal moment of inertia ratios, the average micelle shape calculated from both simulations is a prolate ellipsoid. The deviations at these time scales are presumably the temporary shape change of a micelle. However, the size of the micelle and the accessible surface areas were constant during the simulations with the micelle surface being rough and partially elongated. Radial distribution functions computed for the hydroxyl oxygen atoms of an OG show sharper peaks at a minimum van der Waals contact distance than the acetal oxygen, ring oxygen, and anomeric carbon atoms. This result indicates that these atoms are pointed outwards at the hydrophilic/hydrophobic interface, form hydrogen bonds with the water molecules, and thus hydrate the micelle surface effectively. (c) 2005 Elsevier Inc. All rights reserved.

Molecular dynamics study of MscL interactions with a curved lipid bilayer

Mechanosensitivity is a ubiquitous sensory mechanism found in living organisms. The simplest known mechanotransducing mechanism is found in bacteria in the form of the mechanosensitive membrane channel of large conductance, MscL. This channel has been studied extensively using a variety of methods at a functional and structural level. The channel is gated by membrane tension in the lipid bilayer alone. It serves as a safety valve protecting bacterial cells against hypoosmotic shock. MscL of Escherichia coli embedded in bilayers composed of asymmetric amounts of single-tailed and double-tailed lipids has been shown to gate spontaneously, even in the absence of membrane tension. To gain insight into the effect of the lipid membrane composition and geometry on MscL structure, a fully solvated, all-atom model of MscL in a stress-free curved bilayer composed of double- and single-tailed lipids was studied using a 9.5-ns molecular dynamics simulation. The bilayer was modeled as a domed structure accommodating the asymmetric composition of the monolayers. During the course of the simulation a spontaneous restructuring of the periplasmic loops occurred, leading to interactions between one of the loops and phospholipid headgroups. Previous experimental studies of the role of the loops agree with the observation that opening starts with a restructuring of the periplasmic loop, suggesting an effect of the curved bilayer. Because of limited resources, only one simulation of the large system was performed. However, the results obtained suggest that through the geometry and composition of the bilayer the protein structure can be affected even on short timescales.

Periocular and Orbital Amyloidosis. Clinical Characteristics, Management, and Outcome

Objective: To present the clinical features and management outcome in a large series of patients with periocular and orbital amyloidosis. Design: Retrospective, noncomparative, interventional case series. Patients: All patients diagnosed with periocular and orbital amyloidosis in 6 oculoplastic and orbital units. Methods: Clinical records of all patients were reviewed. Main Outcome Measures: Clinical presentation, radiological and histological findings, treatment modalities, and outcome. Results. The study included 24 patients (15 female, 9 male) with a mean age of 57 17 years. Nineteen cases were unilateral, and 5 were bilateral. Clinical signs and symptoms included a visible or palpable periocular mass or tissue infiltration (95.8%), ptosis (54.2%), periocular discomfort or pain (25%), proptosis or globe displacement (21%), limitations in ocular motility (16.7%), recurrent periocular subcutaneous hemorrhages (12.5%), and diplopia (8.3%). Seven cases had orbital involvement, and 17 were periocular. Immunohistochemistry in 7 patients showed B cells or plasma cells producing monoclonal immunoglobulin chains that were deposited as amyloid light chains. Only 1 patient was diagnosed with systemic amyloid light chain amyloidosis. Treatment modalities were mainly observation and surgical debulking. During a mean follow-up period of 39 months, 21% showed significant progression after treatment, whereas 79% were stable or showed no recurrence after treatment. Conclusion: Periocular and orbital amyloidosis may present with a wide spectrum of clinical findings and result in significant ocular morbidity. Complete surgical excision is not feasible in many cases, and the goal of treatment is to preserve function and to prevent sight-threatening complications.

Idiopathic sclerosing orbital inflammation

Objective: To perform a multicenter review of the clinical features and treatment of 31 patients with idiopathic sclerosing orbital inflammation. Methods: We included all patients with histologically confirmed idiopathic sclerosing orbital inflammation from 5 regional orbital centers. We reviewed the case notes to determine the clinical presentation, diagnostic features, and response to treatment. The main outcome measures were duration and nature of symptoms, anatomical location of disease, histopathological findings, treatment modalities, treatment efficacy and adverse effects, and final clinical status. Results: We included 13 male and 18 female patients ranging in age from 7 to 83 years. The average duration of symptoms at presentation was 13.4 months. There was a predilection for the lateral and superior quadrants. Thirteen patients had apical disease, and 4 had extraorbital involvement. Histopathological findings invariably showed sclerosis associated with a sparse mixed cellular infiltrate. Twenty- seven patients were treated with oral prednisolone, response to which was good in 9 patients, partial in 11, and poor in 7. Six patients were treated with a second- line immunosuppressive agent, and 6 received radiotherapy. The response to radiotherapy was generally poor. Conclusions: Idiopathic sclerosing orbital inflammation is a rare condition that can be difficult to diagnose and manage. Early intervention with immunosuppression in the form of corticosteroids combinedwith secondline agents can result in control and even regression of the disease.

Accelerated leap methods for simulating discrete stochastic chemical kinetics

Biologists are increasingly conscious of the critical role that noise plays in cellular functions such as genetic regulation, often in connection with fluctuations in small numbers of key regulatory molecules. This has inspired the development of models that capture this fundamentally discrete and stochastic nature of cellular biology - most notably the Gillespie stochastic simulation algorithm (SSA). The SSA simulates a temporally homogeneous, discrete-state, continuous-time Markov process, and of course the corresponding probabilities and numbers of each molecular species must all remain positive. While accurately serving this purpose, the SSA can be computationally inefficient due to very small time stepping so faster approximations such as the Poisson and Binomial τ-leap methods have been suggested. This work places these leap methods in the context of numerical methods for the solution of stochastic differential equations (SDEs) driven by Poisson noise. This allows analogues of Euler-Maruyuma, Milstein and even higher order methods to be developed through the Itô-Taylor expansions as well as similar derivative-free Runge-Kutta approaches. Numerical results demonstrate that these novel methods compare favourably with existing techniques for simulating biochemical reactions by more accurately capturing crucial properties such as the mean and variance than existing methods.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

Evaluation of PCR in the molecular diagnosis of endocarditis

Objective. Infective endocarditis (IE) is diagnosed by the Duke criteria, which can be inconclusive particularly when blood cultures are negative. This study investigated the application of polymerase chain reaction (PCR) to identify bacterial DNA in excised valvular tissue, and its role in establishing the diagnosis of IE. Methods. Ninety-eight patients undergoing valve replacement surgery were studied. Twenty-eight patients were confirmed as definite for endocarditis by the Duke criteria; nine were considered as possible and 61 had no known or previous microbial infection of the endocardium. A broad-range PCR technique was used to amplify prokaryotic 16S rRNA genes present within homogenised heart valve tissue. Subsequent DNA sequencing of the PCR amplicon allowed identification of the infecting microorganism. Results. PCR results demonstrated the presence of bacterial DNA in the heart valves obtained from 14 out of 20 (70%) definite IE patients with positive blood cultures preoperatively. The causative microorganism for one patient with definite culture negative endocarditis was identified by PCR. Two out of nine (22%) of the valves from possible endocarditis patients also had bacterial DNA present converting them into the definite criteria whereas in the valves of seven out of nine (78%) of these patients no bacterial DNA was detected. Conclusion. The application of PCR to the explanted valves in patients with possible or confirmed diagnosis can augment the Duke criteria thereby improving post-surgical antimicrobial therapeutic options. © 2003 The British Infection Society. Published by Elsevier Ltd. All rights reserved.

Computational mechanics reveals nanosecond time correlations in molecular dynamics of liquid systems

Statistical complexity, a measure introduced in computational mechanics has been applied to MD simulated liquid water and other molecular systems. It has been found that statistical complexity does not converge in these systems but grows logarithmically without a limit. The coefficient of the growth has been introduced as a new molecular parameter which is invariant for a given liquid system. Using this new parameter extremely long time correlations in the system undetectable by traditional methods are elucidated. The existence of hundreds of picosecond and even nanosecond long correlations in bulk water has been demonstrated. © 2008 Elsevier B.V. All rights reserved.