Optical Doppler tomography based on a field programmable gate array

We report the design of and results obtained by using a field programmable gate array (FPGA) to digitally process optical Doppler tomography signals. The processor fits into the analog signal path in an existing optical coherence tomography setup. We demonstrate both Doppler frequency and envelope extraction using the Hilbert transform, all in a single FPGA. An FPGA implementation has certain advantages over general purpose digital signal processor (DSP) due to the fact that the processing elements operate in parallel as opposed to the DSP, which is primarily a sequential processor.     

Implementation of a field programmable gate array based transcranial Doppler ultrasound system using pulse compression techniques

A field programmable gate array (FPGA) based hardware platform that is used to implement a digital, multi-gate pulsed Doppler ultrasound system for transcranial Doppler (TCD) use is described. The Doppler audio signal is extracted from the digitised radio-frequency signal by matched filtering and suitable sampling. The system was configured to acquire Doppler signals from 16 depth locations and to display Doppler signal power versus depth as well as a sonogram from a user specified depth. The signal-to-noise performance of the system was comparable with that of a commercially available TCD unit for equivalent pulse repetition frequency and sample volume settings.The flexibility of the platform was used to demonstrate the feasibility of using coded transducer excitation and pulse compression techniques to improve axial resolution compared to a non-coded implementation. The axial resolution improvement was demonstrated using a flow phantom and measured using a vibrating wire phantom. The measured resolutions were 9.1 and 2.4 mm for the conventional and coded implementations, respectively. The reduction in signal-to-noise ratio of approximately 5 dB associated with the configuration using coded excitation was attributable to the frequency response characteristics of the transducer rather than the processing technique used. This work demonstrates both the flexibility associated with an FPGA implementation of a Doppler ultrasound system and the potential for coded excitation to improve axial resolution in TCD systems.     

Fast sequence clustering using a suffix array algorithm

MOTIVATION: Efficient clustering is important for handling the large amount of available EST sequences. Most contemporary methods are based on some kind of all-against-all comparison, resulting in a quadratic time complexity. A different approach is needed to keep up with the rapid growth of EST data. RESULTS: A new, fast EST clustering algorithm is presented. Sub-quadratic time complexity is achieved by using an algorithm based on suffix arrays. A prototype implementation has been developed and run on a benchmark data set. The produced clusterings are validated by comparing them to clusterings produced by other methods, and the results are quite promising. AVAILABILITY: The source code for the prototype implementation is available under a GPL license from http://www.ii.uib.no/~ketil/bio/.     

Protein sequence-similarity search acceleration using a heuristic algorithm with a sensitive matrix

Protein database search for public databases is a fundamental step in the target selection of proteins in structural and functional genomics and also for inferring protein structure, function, and evolution. Most database search methods employ amino acid substitution matrices to score amino acid pairs. The choice of substitution matrix strongly affects homology detection performance. We earlier proposed a substitution matrix named MIQS that was optimized for distant protein homology search. Herein we further evaluate MIQS in combination with LAST, a heuristic and fast database search tool with a tunable sensitivity parameter m, where larger m denotes higher sensitivity. Results show that MIQS substantially improves the homology detection and alignment quality performance of LAST across diverse m parameters. Against a protein database consisting of approximately 15 million sequences, LAST with m?=?10⁵ achieves better homology detection performance than BLASTP, and completes the search 20 times faster. Compared to the most sensitive existing methods being used today, CS-BLAST and SSEARCH, LAST with MIQS and m?=?10⁶ shows comparable homology detection performance at 2.0 and 3.9 times greater speed, respectively. Results demonstrate that MIQS-powered LAST is a time-efficient method for sensitive and accurate homology search.     

Display algorithm for space-filling molecular models using a video array processor

A new algorithm for the generation of coloured, space-filling molecular models is described. It relies on the parallel processing facilities of a DeAnza IP8500 raster graphics system to increase the speed of display generation. Using suitable storage of the information about the type as well as the colour of each atom, it becomes possible to switch rapidly between various schemes of colouring. Results are shown with DNA and protein molecules.     

Selective photoreactions in a programmable array microscope (PAM): photoinitiated polymerization, photodecaging, and photochromic conversion.

BACKGROUND: Innovative thinking and experimentation were the hallmarks of Mack Fulwyler's approach to research. This report summarizes some of the ideas and their early realizations that he pursued in the field of imaging cytometry, work that was not published before his untimely death, although he composed the initial draft of this report. METHODS: Included are related experiments implemented in the programmable array microscope (PAM) devised for patterned illumination and detection, the instrument that Mack Fulwyler employed during a sabbatical leave in G?ttingen in 1998. Despite being the originator of instrumentation for flow cytometry and sorting, Mack Fulwyler was intensely interested in imaging systems, recognizing their ability to resolve cellular details obscured by the whole cell signals generally acquired in flow. At one point, these interests merged with those of two other authors (I.T.Y. and T.M.J.), leading to the Image Cytometry and Sorting (ICAS) strategy and project. A major goal was uncomplicated rare cell detection and isolation using a sequential process of cellular labeling via suitable probes, whole field imaging, and selective area-restricted photoinduced reactions designed to encapsulate and/or chemically or physically tag cells in a manner permitting subsequent fractionation by bulk techniques. RESULTS AND CONCLUSION: This publication features photoinduced polymerization, photodecaging, photoactivation, and photochromic conversion reactions carried out by Fulwyler and/or the other authors with the PAM, employing operator designated patterns and locations in various samples. Photopolymerization of polyethylene glycol-diacrylate to a gel-like structure allowing the specific selection of objects (cells) for further analysis and processing techniques was the approach explored personally by Mack Fulwyler in relation to the ICAS concept.     

A forward-backward fragment assembling algorithm for the identification of genomic amplification and deletion breakpoints using high-density single nucleotide polymorphism (SNP) array

Background  

DNA copy number aberration (CNA) is one of the key characteristics of cancer cells. Recent studies demonstrated the feasibility of utilizing high density single nucleotide polymorphism (SNP) genotyping arrays to detect CNA. Compared with the two-color array-based comparative genomic hybridization (array-CGH), the SNP arrays offer much higher probe density and lower signal-to-noise ratio at the single SNP level. To accurately identify small segments of CNA from SNP array data, segmentation methods that are sensitive to CNA while resistant to noise are required.     

Determination of best-fit potential parameters for a reactive force field using a genetic algorithm

The ReaxFF interatomic potential, used for organic materials, involves more than 600 adjustable parameters, the best-fit values of which must be determined for different materials. A new method of determining the set of best-fit parameters for specific molecules containing carbon, hydrogen, nitrogen and oxygen is presented, based on a parameter reduction technique followed by genetic algorithm (GA) minimization. This work has two novel features. The first is the use of a parameter reduction technique to determine which subset of parameters plays a significant role for the species of interest; this is necessary to reduce the optimization space to manageable levels. The second is the application of the GA technique to a complex potential (ReaxFF) with a very large number of adjustable parameters, which implies a large parameter space for optimization. In this work, GA has been used to optimize the parameter set to determine best-fit parameters that can reproduce molecular properties to within a given accuracy. As a test problem, the use of the algorithm has been demonstrated for nitromethane and its decomposition products.     

A "moonlighting" dizinc aminopeptidase from Streptomyces griseus: mechanisms for peptide hydrolysis and the 4 x 10(10)-fold acceleration of the alternative phosphodiester hydrolysis

A unique "enzyme catalytic promiscuity" has recently been observed, wherein a phosphodiester and a phosphonate ester are hydrolyzed by a dinuclear aminopeptidase and its metal derivatives from Streptomyces griseus (SgAP) [Park, H. I., Ming, L.-J. (1999) Angew. Chem., Int. Ed. Engl. 38, 2914-2916 and Ercan, A., Park, H. I., Ming, L.-J. (2000) Chem. Commun. 2501-2502]. Because tetrahedral phosphocenters often serve as transition-state inhibitors toward the hydrolysis of the peptide, phosphoester hydrolysis by peptidases is thus not expected to occur effectively and must take place through a unique mechanism. Owing to the very different structures and mechanistic requirements between phosphoesters and peptides during hydrolysis, the study of this effective phosphodiester hydrolysis by SgAP may provide further insight into the action of this enzyme that is otherwise not obtainable from regular peptide substrates. We present herein a detailed investigation of both peptide and phosphodiester hydrolyses catalyzed by SgAP. The latter exhibits a first-order rate enhancement of 4 x 10(10)-fold compared to the uncatalyzed reaction at pH 7.0 and 25 degrees C. The results suggest that peptide and phosphodiester hydrolyses by SgAP may share a common reaction mechanism to a certain extent. However, their differences in pH dependence, phosphate and fluoride inhibition patterns, and proton inventory reflect that they must follow different pathways. Mechanisms for the two hydrolyses are drawn on the basis of the results, which provide the foundation for further investigation of the catalytic promiscuity of this enzyme by means of physical and molecular biology methods. The catalytic versatility of SgAP suggests that this enzyme may serve as a unique "natural model system" for further investigation of dinuclear hydrolysis. A better understanding of enzyme catalytic promiscuity is also expected to shed light on the evolution and action of enzymes.     

FPGA acceleration of the phylogenetic likelihood function for Bayesian MCMC inference methods

Background  

Likelihood (ML)-based phylogenetic inference has become a popular method for estimating the evolutionary relationships among species based on genomic sequence data. This method is used in applications such as RAxML, GARLI, MrBayes, PAML, and PAUP. The Phylogenetic Likelihood Function (PLF) is an important kernel computation for this method. The PLF consists of a loop with no conditional behavior or dependencies between iterations. As such it contains a high potential for exploiting parallelism using micro-architectural techniques. In this paper, we describe a technique for mapping the PLF and supporting logic onto a Field Programmable Gate Array (FPGA)-based co-processor. By leveraging the FPGA's on-chip DSP modules and the high-bandwidth local memory attached to the FPGA, the resultant co-processor can accelerate ML-based methods and outperform state-of-the-art multi-core processors.     

Considerations when using the significance analysis of microarrays (SAM) algorithm

Background  

Users of microarray technology typically strive to use universally acceptable data analysis strategies to determine significant expression changes in their experiments. One of the most frequently utilised methods for gene expression data analysis is SAM (significance analysis of microarrays). The impact of selection thresholds, on the output from SAM, may critically alter the conclusion of a study, yet this consideration has not been systematically evaluated in any publication.     

Effect of the magnetic field on the resonant particle acceleration

The acceleration of charged particles trapped by a potential wave in a magnetic field is investigated as applied to the problem of the generation of fast particles in a laser plasma. The conditions for unlimited particle acceleration are determined, and the spectra of fast particles are found.     

Development of a DNA sensor using a molecular logic gate

This communication reports the increase in fluorescence resonance energy transfer (FRET) efficiency between two laser dyes in the presence of deoxyribonucleic acid (DNA). Two types of molecular logic gates have been designed where DNA acts as input signal and fluorescence intensity of different bands are taken as output signal. Use of these logic gates as a DNA sensor has been demonstrated.     

Protein structure prediction in a 210-type lattice model: parameter optimization in the genetic algorithm using orthogonal array

We have applied the orthogonal array method to optimize the parameters in the genetic algorithm of the protein folding problem. Our study employed a 210-type lattice model to describe proteins, where the orientation of a residue relative to its neighboring residue is described by two angles. The statistical analysis and graphic representation show that the two angles characterize protein conformations effectively. Our energy function includes a repulsive energy, an energy for the secondary structure preference, and a pairwise contact potential. We used orthogonal array to optimize the parameters of the population, mating factor, mutation factor, and selection factor in the genetic algorithm. By designing an orthogonal set of trials with representative combinations of these parameters, we efficiently determined the optimal set of parameters through a hierarchical search. The optimal parameters were obtained from the protein crambin and applied to the structure prediction of cytochrome B562. The results indicate that the genetic algorithm with the optimal parameters reduces the computing time to reach a converged energy compared to nonoptimal parameters. It also has less chance to be trapped in a local energy minimum, and predicts a protein structure which is closer to the experimental one. Our method may also be applicable to many other optimization problems in computational biology.     

Differential dephosphorylation of the insulin receptor and its 160-kDa substrate (pp160) in rat adipocytes.

A permeabilized rat adipocyte model was developed which permitted an examination of: 1) insulin receptor autophosphorylation, 2) phosphorylation of a putative insulin receptor substrate of 160 kDa, pp160, and 3) the dephosphorylation reactions associated with each of these phosphoproteins. Rat adipocytes, preincubated with [32P]orthophosphate for 2 h, were exposed to insulin (10(-7) M) at the time of digitonin permeabilization. Phosphorylation of pp160 and autophosphorylation of the insulin receptor increased as a function of Mn2+ concentration in the media with near maximum responses at 10 mM. Maximum response was at least as large as the intact cell response to 10(-7) M insulin. In contrast, magnesium did not increase phosphorylation of pp160 although an increase in receptor autophosphorylation was observed. Autophosphorylation was preserved at digitonin concentrations of 20-100 micrograms/ml, but pp160 phosphorylation was negligible beyond 40 micrograms/ml. Our previous work demonstrated that the insulin receptor was associated with a phosphotyrosine phosphatase activity in permeabilized adipocytes (Mooney, R., and Anderson, D. (1989) J. Biol. Chem. 264, 6850-6857). The current permeabilized adipocyte model made possible an examination of the effects of phosphotyrosine phosphatase inhibitors, including several divalent metal cations (Zn2+, Co2+, and Ni2+), vanadate, and molybdate on both net phosphorylation of pp160 and autophosphorylation of the insulin receptor. Zn2+ at 100 microM, Ni2+ at 1 mM, and Co2+ at 1 or 5 mM increased insulin-dependent phosphorylation of pp160 at least 5-fold and autophosphorylation 2-fold. At higher concentrations of Zn2+ (1 mM) and Ni2+ (5 mM), however, no increase in phosphorylation of pp160 was observed and autophosphorylation was inhibited. Vanadate (1 mM) and molybdate (100 microM) increased insulin-dependent phosphorylation of pp160 by 3-fold when tested separately and 7-fold in combination. Insulin receptor autophosphorylation was increased 50% by each and 3-fold when the agents were combined. Dephosphorylation of pp160 and the insulin receptor was analyzed directly by permeabilizing prelabeled insulin-treated adipocytes in the presence of EDTA (10 mM). Dephosphorylation of pp160 was especially rapid with a t1/2 of approximately 10 s. The t1/2 for the insulin receptor was 37 s. Zn2+ at 1 mM (a concentration that inhibited the insulin receptor kinase) was a strong inhibitor of dephosphorylation, prolonging the rate of pp160 dephosphorylation more than 12-fold and insulin receptor dephosphorylation 3-fold.(ABSTRACT TRUNCATED AT 400 WORDS)     

Measurement and analysis of gas exchange during exercise using a programmable calculator

Although exercise testing is useful in the diagnosis and management of cardiovascular and pulmonary diseases, a rapid comprehensive method for measurement of ventilation and gas exchange has been limited to expensive complex computer-based systems. We devised a relatively inexpensive, technically simple, and clinically oriented exercise system built around a desktop calculator. This system automatically collects and analyzes data on a breath-by-breath basis. Our calculator system overcomes the potential inaccuracies of gas exchange measurement due to water vapor dilution and mismatching of expired flow and gas concentrations. We found no difference between the calculator-derived minute ventilation, CO2 production, O2 consumption, and respiratory exchange ratio and the values determined from simultaneous mixed expired gas collections in 30 constant-work-rate exercise studies. Both tabular and graphic displays of minute ventilation, CO2 production, O2 consumption, respiratory exchange ratio, heart rate, end-tidal O2 tension, end-tidal CO2 tension, and arterial blood gas value are included for aid in the interpretation of clinical exercise tests.