A novel harmony search-K means hybrid algorithm for clustering gene expression data

Recent progress in bioinformatics research has led to the accumulation of huge quantities of biological data at various data sources. The DNA microarray technology makes it possible to simultaneously analyze large number of genes across different samples. Clustering of microarray data can reveal the hidden gene expression patterns from large quantities of expression data that in turn offers tremendous possibilities in functional genomics, comparative genomics, disease diagnosis and drug development. The k- ¬means clustering algorithm is widely used for many practical applications. But the original k-¬means algorithm has several drawbacks. It is computationally expensive and generates locally optimal solutions based on the random choice of the initial centroids. Several methods have been proposed in the literature for improving the performance of the k-¬means algorithm. A meta-heuristic optimization algorithm named harmony search helps find out near-global optimal solutions by searching the entire solution space. Low clustering accuracy of the existing algorithms limits their use in many crucial applications of life sciences. In this paper we propose a novel Harmony Search-K means Hybrid (HSKH) algorithm for clustering the gene expression data. Experimental results show that the proposed algorithm produces clusters with better accuracy in comparison with the existing algorithms.     

Application of magnetic techniques in the field of drug discovery and biomedicine

Magnetic separation technology, using magnetic particles, is quick and easy method for sensitive and reliable capture of specific proteins, genetic material and other biomolecules. The technique offers an advantage in terms of subjecting the analyte to very little mechanical stress compared to other methods. Secondly, these methods are non-laborious, cheap and often highly scalable. Moreover, techniques employing magnetism are more amenable to automation and miniaturization. Now that the human genome is sequenced and about 30,000 genes are annotated, the next step is to identify the function of these individual genes, carrying out genotyping studies for allelic variation and SNP analysis, ultimately leading to identification of novel drug targets. In this post-genomic era, technologies based on magnetic separation are becoming an integral part of todays biology laboratory. This article briefly reviews the selected applications of magnetic separation techniques in the field of biotechnology, biomedicine and drug discovery.     

Phenotypic and Genetic Diversity of Aeromonas Species Isolated from Fresh Water Lakes in Malaysia

Gram-negative bacilli of the genus Aeromonas are primarily inhabitants of the aquatic environment. Humans acquire this organism from a wide range of food and water sources as well as during aquatic recreational activities. In the present study, the diversity and distribution of Aeromonas species from freshwater lakes in Malaysia was investigated using glycerophospholipid-cholesterol acyltransferase (GCAT) and RNA polymerase sigma-factor (rpoD) genes for speciation. A total of 122 possible Aeromonas strains were isolated and confirmed to genus level using the API20E system. The clonality of the isolates was investigated using ERIC-PCR and 20 duplicate isolates were excluded from the study. The specific GCAT-PCR identified all isolates as belonging to the genus Aeromonas, in agreement with the biochemical identification. A phylogenetic tree was constructed using the rpoD gene sequence and all 102 isolates were identified as: A. veronii 43%, A. jandaei 37%, A. hydrophila 6%, A. caviae 4%, A. salmonicida 2%, A. media 2%, A. allosaccharophila 1%, A. dhakensis 1% and Aeromonas spp. 4%. Twelve virulence genes were present in the following proportions—exu 96%, ser 93%, aer 87%, fla 83%, enolase 70%, ela 62%, act 54%, aexT 33%, lip 16%, dam 16%, alt 8% and ast 4%, and at least 2 of these genes were present in all 102 strains. The ascV, aexU and hlyA genes were not detected among the isolates. A. hydrophila was the main species containing virulence genes alt and ast either present alone or in combination. It is possible that different mechanisms may be used by each genospecies to demonstrate virulence. In summary, with the use of GCAT and rpoD genes, unambiguous identification of Aeromonas species is possible and provides valuable data on the phylogenetic diversity of the organism.     

Evidence on how a conserved glycine in the hinge region of HapR regulates its DNA binding ability: lessons from a natural variant

HapR has been recognized as a quorum-sensing master regulator in Vibrio cholerae. Because it controls a plethora of disparate cellular events, the absence of a functional HapR affects the physiology of V. cholerae to a great extent. In the current study, we pursued an understanding of an observation of a natural protease-deficient non-O1, non-O139 variant V. cholerae strain V2. Intriguingly, a nonfunctional HapR (henceforth designated as HapR(V2)) harboring a substitution of glycine to aspartate at position 39 of the N-terminal hinge region has been identified. An in vitro gel shift assay clearly suggested the inability of HapR(V2) to interact with various cognate promoters. Reinstatement of glycine at position 39 restores DNA binding ability of HapR(V2) (HapR(V2G)), thereby rescuing the protease-negative phenotype of this strain. The elution profile of HapR(V2) and HapR(V2G) proteins in size-exclusion chromatography and their circular dichroism spectra did not reflect any significant differences to explain the functional discrepancies between the two proteins. To gain insight into the structure-function relationship of these two proteins, we acquired small/wide angle x-ray scattering data from samples of the native and G39D mutant. Although Guinier analysis and indirect Fourier transformation of scattering indicated only a slight difference in the shape parameters, structure reconstruction using dummy amino acids concluded that although HapR adopts a "Y" shape similar to its crystal structure, the G39D mutation in hinge drastically altered the DNA binding domains by bringing them in close proximity. This altered spatial orientation of the helix-turn-helix domains in this natural variant provides the first structural evidence on the functional role of the hinge region in quorum sensing-related DNA-binding regulatory proteins of Vibrio spp.     

Linear Mixed Effects Models under Inequality Constraints with Applications

Constraints arise naturally in many scientific experiments/studies such as in, epidemiology, biology, toxicology, etc. and often researchers ignore such information when analyzing their data and use standard methods such as the analysis of variance (ANOVA). Such methods may not only result in a loss of power and efficiency in costs of experimentation but also may result poor interpretation of the data. In this paper we discuss constrained statistical inference in the context of linear mixed effects models that arise naturally in many applications, such as in repeated measurements designs, familial studies and others. We introduce a novel methodology that is broadly applicable for a variety of constraints on the parameters. Since in many applications sample sizes are small and/or the data are not necessarily normally distributed and furthermore error variances need not be homoscedastic (i.e. heterogeneity in the data) we use an empirical best linear unbiased predictor (EBLUP) type residual based bootstrap methodology for deriving critical values of the proposed test. Our simulation studies suggest that the proposed procedure maintains the desired nominal Type I error while competing well with other tests in terms of power. We illustrate the proposed methodology by re-analyzing a clinical trial data on blood mercury level. The methodology introduced in this paper can be easily extended to other settings such as nonlinear and generalized regression models.     

Effect of boundary type and season on predatory arthropods associated with field margins on New Zealand farmland

Pitfall traps were used to monitor predatory arthropod numbers along two types of field boundary, a post and wire fence line and a Cupressus macrocarpa hedge, along the same paddock margin in Canterbury, New Zealand, over 24 months. The seven most abundant predator groups recorded were: Araneae > Phalangiidae > Staphylinidae > Coccinellidae > Chilopoda > Hemerobiidae > Carabidae. Araneae, Phalangiidae, Staphylinidae, Chilopoda and Hemerobiidae were found in larger numbers at the wire fence than at the hedge site, whereas the numbers of Carabidae and Coccinellidae adults exhibited no field margin preference. However, more species of Araneae and Staphylinidae were caught at the hedge site, whereas species richness of carabid beetles was greatest at the wire fence. Principal component analysis clearly separated the samples collected from the two habitats based on the assemblages of Araneae, Staphylinidae and Carabidae, and certain species in each of these taxonomic groups appeared to be particularly associated with one boundary type or the other. All the main taxonomic groups exhibited clear seasonal patterns, with distinct peaks in abundance occurring at certain times of the year. The results of the study reinforce the idea that management of field boundaries can be used to manipulate the type and abundance of particular groups of predatory arthropods, and that seasonal patterns should be taken into account in schemes of integrated pest management so that any adverse effects of biocide application on these beneficial species may be minimised.     

On the origin of the stereoselective affinity of Nutlin-3 geometrical isomers for the MDM2 protein

The stereoselective affinity of small-molecule binding to proteins is typically broadly explained in terms of the thermodynamics of the final bound complex. Using Brownian dynamics simulations, we show that the preferential binding of the MDM2 protein to the geometrical isomers of Nutlin-3, an effective anticancer lead that works by inhibiting the interaction between the proteins p53 and MDM2, can be explained by kinetic arguments related to the formation of the MDM2:Nutlin-3 encounter complex. This is a diffusively bound state that forms prior to the final bound complex. We find that the MDM2 protein stereoselectivity for the Nutlin-3a enantiomer stems largely from the destabilization of the encounter complex of its mirror image enantiomer Nutlin-3b, by the K70 residue that is located away from the binding site. On the other hand, the trans-Nutlin-3a diastereoisomer exhibits a shorter residence time in the vicinity of MDM2 compared with Nutlin-3a due to destabilization of its encounter complex by the collective interaction of pairs of charged residues on either side of the binding site: Glu25 and Lys51 on one side, and Lys94 and Arg97 on the other side. This destabilization is largely due to the electrostatic potential of the trans-Nutlin-3a isomer being largely positive over extended continuous regions around its structure, which are otherwise well-identified into positive and negative regions in the case of the Nutlin-3a isomer. Such rich insight into the binding processes underlying biological selectivity complements the static view derived from the traditional thermodynamic analysis of the final bound complex. This approach, based on an explicit consideration of the dynamics of molecular association, suggests new avenues for kinetics-based anticancer drug development and discovery.     

A geometric approach to determine association and coherence of the activation times of cell-cycling genes under differing experimental conditions

Differing arresting agents and protocols can be used to synchronize cells in cultures to specific phases of the cell when studying cell-cycle gene expressions. Often, data derived from individual experiments are analyzed separately, since no appropriate statistical methodology is available at the moment to analyze the data from all such experiments simultaneously. The focus of this paper is to determine the association and coherence of the relative activation times of cell-cycling genes under different experimental conditions. Using a circular-circular regression model, we define two parameters, a rotation parameter for the angular difference between cells' arresting times (phases) in two cell-cycle experiments, and an association parameter to describe the correspondence between the cycle times of maximal expression (phase angles) for a set of genes studied in two experiments. Further, we propose a procedure to assess coherence across multiple experiments, i.e. to what extent the circular ordering of the phase angles of genes is maintained across multiple experiments. Coherence of genes across experiments suggests that functionally these genes tend to respond in a stereotypically sequenced way under different experimental conditions. Our proposed methodology is illustrated by applying it to a HeLa cell-cycle gene-expression data.     

Male sex hormones promote vagally mediated reflex airway responsiveness to cholinergic stimulation

A sex disparity in airway responsiveness to cholinergic stimulation has been observed in laboratory mice in that males are considerably more responsive than females, but the basis for this difference is unclear. In this report, we demonstrate that male sex hormones promote murine airway responsiveness to cholinergic stimulation via vagus nerve-mediated reflex mechanisms. In tissue bath preparations, no sex-based differences were observed in the contractile responses of isolated tracheal and bronchial ring segments to carbachol, indicating that the mechanism(s) responsible for the in vivo sex difference is (are) absent ex vivo. Bilateral cervical vagotomy was found to abolish in vivo airway responsiveness to methacholine in male mice, whereas it did not alter the responses of females, suggesting a regulatory role for male sex hormones in promoting reflex airway constriction. To test this possibility, we next studied mice with altered circulating male sex hormone levels. Castrated male mice displayed airway responsiveness equivalent to that observed in intact females, whereas administration of exogenous testosterone to castrated males restored responsiveness, albeit not to the level observed in intact males. Administration of exogenous testosterone to intact female mice similarly enhanced responsiveness. Importantly, the promotive effects of exogenous testosterone in castrated male and intact female mice were absent when bilateral vagotomy was performed. Together, these data indicate that male sex hormones promote cholinergic airway responsiveness via a vagally mediated reflex mechanism that may be important in the regulation of airway tone in the normal and diseased lung.