A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, first a Fortran code is developed for three dimensional linear elastostatics using constant boundary elements; the code is based on a MATLAB code developed by the author earlier. Next, the code is parallelized using BLACS, MPI, and ScaLAPACK. Later, the parallelized code is used to demonstrate the usefulness of the Boundary Element Method (BEM) as applied to the realtime computational simulation of biological organs, while focusing on the speed and accuracy offered by BEM. A computer cluster is used in this part of the work. The commercial software package ANSYS is used to obtain the `exact' solution against which the solution from BEM is compared; analytical solutions, wherever available, are also used to establish the accuracy of BEM. A pig liver is the biological organ considered. Next, instead of the computer cluster, a Graphics Processing Unit (GPU) is used as the parallel hardware. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature. Also, a serial MATLAB code, and both serial and parallel versions of a Fortran code, which can solve three dimensional (3D) linear elastostatic problems using constant boundary elements, are provided as supplementary files that can be freely downloaded.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated, with specific reference to the speed and the accuracy offered by BEM. First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used. A pig liver is the biological organ considered. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.

A Study of Speed of the Boundary Element Method as applied to the Realtime Computational Simulation of Biological Organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated. Biological organs are assumed to follow linear elastostatic material behavior, and constant boundary element is the element type used. First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used. Results indicate that BEM is fast enough to provide for realtime graphics if biological organs are assumed to follow linear elastostatic material behavior. Although the present work does not conduct any simulation using nonlinear material models, results from using the linear elastostatic material model imply that it would be difficult to obtain realtime performance if highly nonlinear material models that properly characterize biological organs are used. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.

A method to improve size estimates of walleye pollock (Theragra chalcogramma) Atka mackerel (Pleurogrammus monopterygius) consumed by pinnipeds: digestion correction factors applied to bones and otoliths recovered in scats

The lengths of otoliths and other skeletal structures recovered from the scats of pinnipeds, such as Steller sea lions (Eumetopias jubatus), correlate with body size and can be used to estimate the length of prey consumed. Unfortunately, otoliths are often found in too few scats or are too digested to usefully estimate prey size. Alternative diagnostic bones are frequently recovered, but few bone-size to prey-size correlations exist and bones are also reduced in size by various degrees owing to digestion. To prevent underestimates in prey sizes consumed techniques are required to account for the degree of digestion of alternative bones prior to estimating prey size. We developed a method (using defined criteria and photo-reference material) to assign the degree of digestion for key cranial structures of two prey species: walleye pollock (Theragra chalcogramma) and Atka mackerel (Pleurogrammus monopterygius). The method grades each structure into one of three condition categories; good, fair or poor. We also conducted feeding trials with captive Steller sea lions, feeding both fish species to determine the extent of erosion of each structure and to derive condition-specific digestion correction factors to reconstruct the original sizes of the structures consumed. In general, larger structures were relatively more digested than smaller ones. Mean size reduction varied between different types of structures (3.3−26.3%), but was not influenced by the size of the prey consumed. Results from the observations and experiments were combined to be able to reconstruct the size of prey consumed by sea lions and other pinnipeds. The proposed method has four steps: 1) measure the recovered structures and grade the extent of digestion by using defined criteria and photo-reference collection; 2) exclude structures graded in poor condition; 3) multiply measurements of structures in good and fair condition by their appropriate digestion correction factors to derive their original size; and 4) calculate the size of prey from allometric regressions relating corrected structure measurements to body lengths. This technique can be readily applied to piscivore dietary studies that use hard remains of fish.

Artificial neural network applied to diagnosis of lung cancer

Artificial neural network(ANN) approach was applied to classification of normal persons and lung cancer patients based on the metal content of hair and serum samples obtained by inductively coupled plasma atomic emission spectrometry (ICP-AES) for the two groups. This method was verified with independent prediction samples and can be used as an aiding means of the diagnosis of lung cancer. The case of predictive classification with one element missing in the prediction samples was studied in details, The significance of elements in hair and serum samples for classification prediction was also investigated.

Density-viscosity product of small-volume ionic liquid samples using quartz crystal impedance analysis

Quartz crystal impedance analysis has been developed as a technique to assess whether room-temperature ionic liquids are Newtonian fluids and as a small-volume method for determining the values of their viscosity-density product, rho eta. Changes in the impedance spectrum of a 5-MHz fundamental frequency quartz crystal induced by a water-miscible room-temperature ionic liquid, 1-butyl-3-methylimidazolium. trifluoromethylsulfonate ([C(4)mim][OTf]), were measured. From coupled frequency shift and bandwidth changes as the concentration was varied from 0 to 100% ionic liquid, it was determined that this liquid provided a Newtonian response. A second water-immiscible ionic liquid, 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide [C(4)mim][NTf2], with concentration varied using methanol, was tested and also found to provide a Newtonian response. In both cases, the values of the square root of the viscosity-density product deduced from the small-volume quartz crystal technique were consistent with those measured using a viscometer and density meter. The third harmonic of the crystal was found to provide the closest agreement between the two measurement methods; the pure ionic liquids had the largest difference of similar to 10%. In addition, 18 pure ionic liquids were tested, and for 11 of these, good-quality frequency shift and bandwidth data were obtained; these 12 all had a Newtonian response. The frequency shift of the third harmonic was found to vary linearly with square root of viscosity-density product of the pure ionic liquids up to a value of root(rho eta) approximate to 18 kg m(-2) s(-1/2), but with a slope 10% smaller than that predicted by the Kanazawa and Gordon equation. It is envisaged that the quartz crystal technique could be used in a high-throughput microfluidic system for characterizing ionic liquids.

Conditional probability applied to the definition of limestone geochemical facies: A new Methodological Approach

A new method, based on linear correlation and phase diagrams was successfully developed for processes like the sedimentary process, where the deposition phase can have different time duration - represented by repeated values in a series - and where the erosion can play an important rule deleting values of a series. The sampling process itself can be the cause of repeated values - large strata twice sampled - or deleted values: tiny strata fitted between two consecutive samples. What we developed was a mathematical procedure which, based upon the depth chemical composition evolution, allows the establishment of frontiers as well as the periodicity of different sedimentary environments. The basic tool isn't more than a linear correlation analysis which allow us to detect the existence of eventual evolution rules, connected with cyclical phenomena within time series (considering the space assimilated to time), with the final objective of prevision. A very interesting discovery was the phenomenon of repeated sliding windows that represent quasi-cycles of a series of quasi-periods. An accurate forecast can be obtained if we are inside a quasi-cycle (it is possible to predict the other elements of the cycle with the probability related with the number of repeated and deleted points). We deal with an innovator methodology, reason why it's efficiency is being tested in some case studies, with remarkable results that shows it's efficacy. Keywords: sedimentary environments, sequence stratigraphy, data analysis, time-series, conditional probability.

Quantification of clenbuterol at trace level in human urine by ultra-high pressure liquid chromatography-tandem mass spectrometry.

Clenbuterol is a β2 agonist agent with anabolic properties given by the increase in the muscular mass in parallel to the decrease of the body fat. For this reason, the use of clenbuterol is forbidden by the World Anti-Doping Agency (WADA) in the practice of sport. This compound is of particular interest for anti-doping authorities and WADA-accredited laboratories due to the recent reporting of risk of unintentional doping following the eating of meat contaminated with traces of clenbuterol in some countries. In this work, the development and the validation of an ultra-high pressure liquid chromatography coupled to electrospray ionization tandem mass spectrometry (UHPLC-ESI-MS/MS) method for the quantification of clenbuterol in human urine is described. The analyte was extracted from urine samples by liquid-liquid extraction (LLE) in basic conditions using tert butyl-methyl ether (TBME) and analyzed by UHPLC-MS/MS with a linear gradient of acetonitrile in 9min only. The simple and rapid method presented here was validated in compliance with authority guidelines and showed a limit of quantification at 5pg/mL and a linearity range from 5pg/mL to 300pg/mL. Good trueness (85.8-105%), repeatability (5.7-10.6% RSD) and intermediate precision (5.9-14.9% RSD) results were obtained. The method was then applied to real samples from eighteen volunteers collecting urines after single oral doses administration (1, 5 and 10μg) of clenbuterol-enriched yogurts.

Survival analysis statistics applied to threshold data obtained from the ascending forced-choice method of limits

Sensory thresholds are often collected through ascending forced-choice methods. Group thresholds are important for comparing stimuli or populations; yet, the method has two problems. An individual may correctly guess the correct answer at any concentration step and might detect correctly at low concentrations but become adapted or fatigued at higher concentrations. The survival analysis method deals with both issues. Individual sequences of incorrect and correct answers are adjusted, taking into account the group performance at each concentration. The technique reduces the chance probability where there are consecutive correct answers. Adjusted sequences are submitted to survival analysis to determine group thresholds. The technique was applied to an aroma threshold and a taste threshold study. It resulted in group thresholds similar to ASTM or logarithmic regression procedures. Significant differences in taste thresholds between younger and older adults were determined. The approach provides a more robust technique over previous estimation methods.

Coupling liquid MALDI MS to liquid chromatography

Matrix-assisted laser desorption/ionisation (MALDI) coupled with time-of-flight (TOF) mass spectrometry (MS) is a powerful tool for the analysis of biological samples, and nanoflow high-performance liquid chromatography (nanoHPLC) is a useful separation technique for the analysis of complex proteomics samples. The off-line combination of MALDI and nanoHPLC has been extensively investigated and straightforward techniques have been developed, focussing particularly on automated MALDI sample preparation that yields sensitive and reproducible spectra. Normally conventional solid MALDI matrices such as α-cyano-4-hydroxycinnamic acid (CHCA) are used for sample preparation. However, they have limited usefulness in quantitative measurements and automated data acquisition because of the formation of heterogeneous crystals, resulting in highly variable ion yields and desorption/ ionization characteristics. Glycerol-based liquid support matrices (LSM) have been proposed as an alternative to the traditional solid matrices as they provide increased shot-to-shot reproducibility, leading to prolonged and stable ion signals and therefore better results. This chapter focuses on the integration of the liquid LSM MALDI matrices into the LC-MALDI MS/MS approach in identifying complex and large proteomes. The interface between LC and MALDI consists of a robotic spotter, which fractionates the eluent from the LC column into nanoliter volumes, and co-spots simultaneously the liquid matrix with the eluent fractions onto a MALDI target plate via sheath flow. The efficiency of this method is demonstrated through the analysis of trypsin digests of both bovine serum albumin (BSA) and Lactobacillus plantarum WCFS1 proteins.

Spectrophotometric determination of zinc in pharmaceutical and biological samples with di-2-pyridyl ketone benzoylhydrazone: A simple, fast, accurate and sensitive method

A simple, fast, accurate, and sensitive spectrophotometric method was developed to determine zinc(II). This method is based on the reaction of Zn(II) with di-2-pyridyl ketone benzoylhydrazone (DPKBH), at pH=5.5 and 50% (v/v) ethanol. Beers law was obeyed in the range 0.020-1.82 mu g mL(-1) with a molar apsorptivity of 3.64 x 10(4) L mol(-1) cm(-1), and a detection limit (3) of 2.29 mu g L-1. The action of some interfering ions was verified and the developed method applied to pharmaceutical and biological samples. The results were then compared with those obtained by using a flame atomic absorption technique.

Micro-Raman spectroscopy and SEM/EDX applied to improve the zircon fission track method used for dating geological formations

The zircon mineral is widely studied in geochronology. In the case of the fission track method (FTM), the age is determined by the density of fission tracks at the zircon surface, which can be observed with an optical microscope after an appropriate chemical treatment (etching). The etching must be isotropic at the zircon grain surface to be used in the FTM, which leads those zircon grains whose etching is anisotropic to be discarded. The only reason for this discarding is the nonuniform morphology of the surface grain seen by optical microscopy, that is, no further physicochemical analysis is performed. In this work, combining micro-Raman and scanning electron microscopy (SEM) to study the etching anisotropy, it was shown that zircon grains that present at least one area at the surface where the density of fission track is uniform can be used in the FTM. The micro-Raman showed characteristic spectra of the standard zircon sample either from the areas where there are tracks or from where there are not. The only difference found was in the Raman bandwidths, which were broader for the areas with higher density of fission tracks. This suggests simply a decrease in the relative percentage of the crystalline/amorphous phases at these areas. The SEM/energy dispersive spectrometry (EDX) showed that there were no significant differences in the principal chemical composition at the areas with and without fission tracks. Copyright (c) 2008 John Wiley & Sons, Ltd.

Staining methods applied to glycol methacrylate embedded tissue sections

The use of glycol methacrylate (GMA) avoids some technical artifacts, which are usually observed in paraffin-embedded sections, providing good morphological resolution. on the other hand, weak staining have been mentioned during the use of different methods in plastic sections. In the present study, changes in the histological staining procedures have been assayed during the use of staining and histochemical methods in different GMA-embedded tissues.Samples of tongue, submandibular and sublingual glands, cartilage, portions of respiratory tract and nervous ganglion were fixed in 4% formaldehyde and embedded in glycol methacrylate. The sections of tongue and nervous ganglion were stained by H&E. Picrosirius, Toluidine Blue and Sudan Black B methods were applied, respectively, for identification of collagen fibers in submandibular gland, sulfated glycosaminoglycans in cartilage (metachromasia) and myelin lipids in nervous ganglion. Periodic Acid-Schiff (PAS) method was used for detection of glycoconjugates in submandibular gland and cartilage while AB/PAS combined methods were applied for detection of mucins in the respiratory tract. In addition, a combination of Alcian Blue (AB) and Picrosirius methods was also assayed in the sublingual gland sections.The GMA-embedded tissue sections showed an optimal morphological integrity and were favorable to the staining methods employed in the present study. In the sections of tongue and nervous ganglion, a good contrast of basophilic and acidophilic structures was obtained by H&E. An intense eosinophilia was observed either in the striated muscle fibers or in the myelin sheaths in which the lipids were preserved and revealed by Sudan Black B. In the cartilage matrix, a strong metachromasia was revealed by Toluidine Blue in the negatively-charged glycosaminoglycans. In the chondrocytes, glycogen granules were intensely positive to PAS method. Extracellular glycoproteins were also PAS positive in the basal membrane and in the region occupied by the lamina externa and reticular fibers surrounding each smooth muscle cells of the blood vessels. In the epithelial cells of the respiratory tract, acid and neutral mucins were histochemically detected by AB and PAS methods, respectively. Moreover, granules containing acid and neutral mucins were revealed in purple by AB and PAS concomitantly. In the sublingual gland sections, a distinct affinity of acid mucins by AB (in turquoise-blue) and collagen fibers by Picrosirius (in red) was obtained when these methods were combined. Although some routine dyes used in paraffin sections have showed a weak stain in historesin sections, our results showed that different dyes could be applied in GMA sections if modified staining procedures were assayed. Therefore, appropriate staining contrast and, thus, detection of one or different substances in a same section can be acquired in association to the good morphological resolution provided by GMA. (C) 2003 Elsevier Ltd. All rights reserved.

Phase diagrams of a CTAB/organic solvent/buffer system applied to extraction of enzymes by reverse micelles

A partial pseudo-ternary phase diagram has been studied for the cethyltrimethylammonium bromide/isooctane:hexanol:butanol/potassium phosphate buffer system, where the two-phase diagram consisting of the reverse micelle phase (L-2) in equilibrium with the solvent is indicated. Based on these diagrams two-phase systems of reverse micelles were prepared with different compositions of the compounds and used for extraction and recovery of two enzymes, and the percentage of enzyme recovery yield monitored. The enzymes glucose-6-phosphate dehydrogenase (G6PD) and xylose redutase (XR) obtained from Candida guilliermondii yeast were used in the extraction procedures. The recovery yield data indicate that micelles having different composition give selective extraction of enzymes. The method can thus be used to optimize enzyme extraction processes. (c) 2007 Elsevier B.V. All rights reserved.

Screen-Printed Carbon Electrode Modified with Poly-L-histidine Applied to Gold(III) Determination

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)