Speed and accuracy improvement of higher-order epistasis detection on CUDA-enabled GPUs

The discovery of higher-order epistatic interactions is an important task in the field of genome wide association studies which allows for the identification of complex interaction patterns between multiple genetic markers. Some existing bruteforce approaches explore the whole space of k-interactions in an exhaustive manner resulting in almost intractable execution times. Computational cost can be reduced drastically by restricting the search space with suitable preprocessing filters which prune unpromising candidates. Other approaches mitigate the execution time by employing massively parallel accelerators in order to benefit from the vast computational resources of these architectures. In this paper, we combine a novel preprocessing filter, namely SingleMI, with massively parallel computation on modern GPUs to further accelerate epistasis discovery. Our implementation improves both the runtime and accuracy when compared to a previous GPU counterpart that employs mutual information clustering for prefiltering. SingleMI is open source software and publicly available at: https://github.com/sleeepyjack/singlemi/.     

CUDASW++2.0: enhanced Smith-Waterman protein database search on CUDA-enabled GPUs based on SIMT and virtualized SIMD abstractions

Background

Ca²⁺-binding proteins are important for the transduction of Ca²⁺ signals into physiological outcomes. As in calmodulin many of the Ca²⁺-binding proteins bind Ca²⁺ through EF-hand motifs. Amongst the large number of EF-hand containing Ca²⁺-binding proteins are a subfamily expressed in neurons and retinal photoreceptors known as the CaBPs and the related calneuron proteins. These were suggested to be vertebrate specific but exactly which family members are expressed outside of mammalian species had not been examined.

Findings

We have carried out a bioinformatic analysis to determine when members of this family arose and the conserved aspects of the protein family. Sequences of human members of the family obtained from GenBank were used in Blast searches to identify corresponding proteins encoded in other species using searches of non-redundant proteins, genome sequences and mRNA sequences. Sequences were aligned and compared using ClustalW. Some families of Ca²⁺-binding proteins are known to show a progressive expansion in gene number as organisms increase in complexity. In contrast, the results for CaBPs and calneurons showed that a full complement of CaBPs and calneurons are present in the teleost fish Danio rerio and possibly in cartilaginous fish. These findings suggest that the entire family of genes may have arisen at the same time during vertebrate evolution. Certain members of the family (for example the short form of CaBP1 and calneuron 1) are highly conserved suggesting essential functional roles.

Conclusions

The findings support the designation of the calneurons as a distinct sub-family. While the gene number for CaBPs/calneurons does not increase, a distinctive evolutionary change in these proteins in vertebrates has been an increase in the number of splice variants present in mammals.     

Smoldyn on graphics processing units: massively parallel Brownian dynamics simulations

Space is a very important aspect in the simulation of biochemical systems; recently, the need for simulation algorithms able to cope with space is becoming more and more compelling. Complex and detailed models of biochemical systems need to deal with the movement of single molecules and particles, taking into consideration localized fluctuations, transportation phenomena, and diffusion. A common drawback of spatial models lies in their complexity: models can become very large, and their simulation could be time consuming, especially if we want to capture the systems behavior in a reliable way using stochastic methods in conjunction with a high spatial resolution. In order to deliver the promise done by systems biology to be able to understand a system as whole, we need to scale up the size of models we are able to simulate, moving from sequential to parallel simulation algorithms. In this paper, we analyze Smoldyn, a widely diffused algorithm for stochastic simulation of chemical reactions with spatial resolution and single molecule detail, and we propose an alternative, innovative implementation that exploits the parallelism of Graphics Processing Units (GPUs). The implementation executes the most computational demanding steps (computation of diffusion, unimolecular, and bimolecular reaction, as well as the most common cases of molecule-surface interaction) on the GPU, computing them in parallel on each molecule of the system. The implementation offers good speed-ups and real time, high quality graphics output     

Second international conference on accelerating biopharmaceutical development

The Second International Conference on Accelerating Biopharmaceutical Development was held in Coronado, California. The meeting was organized by the Society for Biological Engineering (SBE) and the American Institute of Chemical Engineers (AIChE); SBE is a technological community of the AIChE. Bob Adamson (Wyeth) and Chuck Goochee (Centocor) were co‐chairs of the event, which had the theme “Delivering cost‐effective, robust processes and methods quickly and efficiently.” The first day focused on emerging disruptive technologies and cutting‐edge analytical techniques. Day two featured presentations on accelerated cell culture process development, critical quality attributes, specifications and comparability, and high throughput protein formulation development. The final day was dedicated to discussion of technology options and new analysis methods provided by emerging disruptive technologies; functional interaction, integration and synergy in platform development; and rapid and economic purification process development. © This meeting report was written for and published by Landes Bioscience in the journal, mAbs. Reichert JM, Jacob N, Amanullah A. Second International Conference on Accelerating Biopharmaceutical Development: March 9–12, 2009, Coronado, CA USA. mAbs 2009 1(3); http://www. landesbioscience.com/journals/mabs/article/8491 Biotechnol. Prog., 2009     

Plant growth using EMCS hardware on the ISS

Under separate contracts with ESA (FUMO and ERM Study) and as a link in the development of the European Modular Cultivation System's (EMCS) functionality and biocompatibility, plant studies have been performed at The Plant Biocentre in Trondheim, Norway. The main goal was to test whether the breadboards containing the major components planned for use in the EMCS would be optimal for space experiments with plant material. The test plans and the experimental set-up for the verification of biocompatibility and biological functionality included the use of a few model plant species including cress (Lepidium sativum L.) and Arabidopsis thaliana. The plants were tested at different developmental levels of morphological and physiological complexity (illumination, life support, humidity control, water supply, observation, short- and long-term plant growth experiments and contamination prevention). Results from the tests show that the EMCS concept is useful for long duration plant growth on the ISS.     

Energy efficiency in cloud computing data center: a survey on hardware technologies

Cluster Computing - The internet is expanding its viewpoint into each conceivable part of the cutting-edge economy. Unshackled from our web programs today, the internet is characterizing our way of...