Assessment Method for a Power Analysis to Identify Differentially Expressed Pathways

Gene expression data can provide a very rich source of information for elucidating the biological function on the pathway level if the experimental design considers the needs of the statistical analysis methods. The purpose of this paper is to provide a comparative analysis of statistical methods for detecting the differentially expression of pathways (DEP). In contrast to many other studies conducted so far, we use three novel simulation types, producing a more realistic correlation structure than previous simulation methods. This includes also the generation of surrogate data from two large-scale microarray experiments from prostate cancer and ALL. As a result from our comprehensive analysis of 41,004 parameter configurations, we find that each method should only be applied if certain conditions of the data from a pathway are met. Further, we provide method-specific estimates for the optimal sample size for microarray experiments aiming to identify DEP in order to avoid an underpowered design. Our study highlights the sensitivity of the studied methods on the parameters of the system. © 2012 Tripahti and Emmert-Streib.

Harnessing the complexity of gene expression data from cancer: from single gene to structural pathway methods

High-dimensional gene expression data provide a rich source of information because they capture the expression level of genes in dynamic states that reflect the biological functioning of a cell. For this reason, such data are suitable to reveal systems related properties inside a cell, e.g., in order to elucidate molecular mechanisms of complex diseases like breast or prostate cancer. However, this is not only strongly dependent on the sample size and the correlation structure of a data set, but also on the statistical hypotheses tested. Many different approaches have been developed over the years to analyze gene expression data to (I) identify changes in single genes, (II) identify changes in gene sets or pathways, and (III) identify changes in the correlation structure in pathways. In this paper, we review statistical methods for all three types of approaches, including subtypes, in the context of cancer data and provide links to software implementations and tools and address also the general problem of multiple hypotheses testing. Further, we provide recommendations for the selection of such analysis methods.

Gene Sets Net Correlations Analysis (GSNCA): a multivariate differential coexpression test for gene sets

Motivation: To date, Gene Set Analysis (GSA) approaches primarily focus on identifying differentially expressed gene sets (pathways). Methods for identifying differentially coexpressed pathways also exist but are mostly based on aggregated pairwise correlations, or other pairwise measures of coexpression. Instead, we propose Gene Sets Net Correlations Analysis (GSNCA), a multivariate differential coexpression test that accounts for the complete correlation structure between genes.

Results: In GSNCA, weight factors are assigned to genes in proportion to the genes' cross-correlations (intergene correlations). The problem of finding the weight vectors is formulated as an eigenvector problem with a unique solution. GSNCA tests the null hypothesis that for a gene set there is no difference in the weight vectors of the genes between two conditions. In simulation studies and the analyses of experimental data, we demonstrate that GSNCA, indeed, captures changes in the structure of genes' cross-correlations rather than differences in the averaged pairwise correlations. Thus, GSNCA infers differences in coexpression networks, however, bypassing method-dependent steps of network inference. As an additional result from GSNCA, we define hub genes as genes with the largest weights and show that these genes correspond frequently to major and specific pathway regulators, as well as to genes that are most affected by the biological difference between two conditions. In summary, GSNCA is a new approach for the analysis of differentially coexpressed pathways that also evaluates the importance of the genes in the pathways, thus providing unique information that may result in the generation of novel biological hypotheses.

On the structure and dynamics of ionic liquids

The structure and dynamics of the ionic liquid 1-ethyl-3-methylimidazolium nitrate is studied by molecular dynamics simulations. We find long-range spatial correlations between the ions and a three-dimensional local structure that reflects the asymmetry of the cations. The main contribution to the configurational energy comes from the electrostatic interactions which leads to charge-ordering effects. Radial screening and threedimensional distribution of charge are also analyzed. The motion of a single ion is studied via velocity and reorientational correlation functions. It is found that ions "rattle" in a long-lived cage, while the orientational structure relaxes on a time scale longer than 200 ps. As in a supercooled liquid, the mean square displacements reveal a subdiffusive dynamics. In addition, the presence of dynamic heterogeneities can be detected by analyzing the non-Gaussian behavior of the van Hove correlation function and the spatial arrangement of the most mobile ions. The short-time collective dynamics is also studied through the electric current time correlation function.

Structure and dynamics of a confined ionic liquid. topics of relevance to dye-sensitized solar cells

The behavior of a model ionic liquid (IL) confined between two flat parallel walls was studied at various interwall distances using computer simulations. The results focus both on structural and dynamical properties. Mass and charge density along the confinement axis reveal a structure of layers parallel to the walls that leads to an oscillatory profile in the electrostatic potential. Orientational correlation functions indicate that cations at the interface orient tilted with respect to the surface and that any other orientational order is lost thereafter. The diffusion coefficients of the ions exhibit a maximum as a function of the confinement distance, a behavior that results from a combination of the structure of the liquid as a whole and a faster molecular motion in the vicinity of the walls. We discuss the relevance of the present results and elaborate on topics that need further attention regarding the effects of ILs in the functioning of IL-based dye-sensitized solar cells.

Thermodynamics, structure, and dynamics in room temperature ionic liquids: The case of 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6])

A detailed investigation of the phase diagram of 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6]) is presented on the basis of a wide set of experimental data accessing thermodynamic, structural, and dynamical properties of this important room temperature ionic liquid (RTIL). The combination of quasi adiabatic, continuous calorimetry, wide angle neutron and X-ray diffraction, and quasi elastic neutron scattering allows the exploration of many novel features of this material. Thermodynamic and microscopic structural information is derived on both glassy and crystalline states and compared with results that recently appeared in the literature allowing direct information to be obtained on the existence of two crystalline phases that were not previously characterized and confirming the view that RTILs show a substantial degree of order (even in their amorphous states), which resembles the crystalline order. We highlight a strong connection between structure and dynamics, showing the existence of three temperature ranges in the glassy state across which both the spatial correlation and the dynamics change. The complex crystalline polymorphism in [bmim][PF6] also is investigated; we compare our findings with the corresponding findings for similar RTILs. These results provide a strong experimental basis for the exploration of the features of the phase diagram of RTILs and for the further study of longer alkyl chain salts.

Liquid structure of the ionic liquid 1,3-dimethylimidazolium bis{(trifluoromethyl)sulfonyl}amide

Neutron diffraction has been used to determine the liquid structure of 1,3-dimethylimidazolium bis{( trifluoromethyl) sulfonyl} amide ([dmim][NTf2]). Significantly smaller charge ordering is found in this liquid compared with analogous chloride and hexafluorophosphate salts due to the diffuse charge density and size of the [NTf2](-) anion. This is manifested in a much larger cation-cation and cation-anion separation and an overlap of the cation-cation and cation-anion shells. Comparison of the liquid structure with the crystal structure reported by Holbrey et al. ( Dalton Trans. 2004, 2267) indicates little correlation, for example, the [NTf2](-) anion adopts a trans orientation predominantly in the liquid whereas a cis orientation is found in the solid phase.

Development of a QSPR correlation for the parachor of 1,3-dialkyl imidazolium based ionic liquids

Ionic liquids have received significant interest from both academia and industry for a wide range of applications which often requires knowledge of their thermophysical properties. Quantitative structure-property relationship correlations and group contribution methods for thermophysical properties of ionic liquids are a basic necessity for the development of computer aided molecular design approaches for these liquids and subsequently offer the potential for designing an ionic liquid having a desirable set of thermophysical properties. However, the limited availability of experimental thermophysical data and their quality have prevented the development of such tools. Based on previously reported experimental surface tension data, a correlation of the parachors with the molar volume of the ionic liquids has been developed. The predicted parachor values have been shown to be in good agreement with the experimental data. A maximum deviation of

Thermodynamic properties and atomistic structure of the dry amorphous silica surface from a reactive force field model

A force field model of the Keating type supplemented by rules to break, form, and interchange bonds is applied to investigate thermodynamic and structural properties of the amorphous SiO2 surface. A simulated quench from the liquid phase has been carried out for a silica sample made of 3888 silicon and 7776 oxygen atoms arranged on a slab similar to 40 angstrom thick, periodically repeated along two directions. The quench results into an amorphous sample, exposing two parallel square surfaces of similar to 42 nm(2) area each. Thermal averages computed during the quench allow us to determine the surface thermodynamic properties as a function of temperature. The surface tension turns out to be gamma=310 +/- 20 erg/cm(2) at room temperature and gamma=270 +/- 30 at T=2000 K, in fair agreement with available experimental estimates. The entropy contribution Ts-s to the surface tension is relatively low at all temperatures, representing at most similar to 20% of the surface energy. Almost without exceptions, Si atoms are fourfold coordinated and oxygen atoms are twofold coordinated. Twofold and threefold rings appear only at low concentration and are preferentially found in proximity of the surface. Above the glass temperature T-g=1660 +/- 50 K, the mobility of surface atoms is, as expected, slightly higher than that of bulk atoms. The computation of the height-height correlation function shows that the silica surface is rough in the equilibrium and undercooled liquid phase, becoming smooth below the glass temperature T-g.

Structure observed in measured electron energy distribution functions in capacitively coupled radio frequency hydrogen plasmas

Electron energy probability functions measured with a passively compensated Langmuir probe in asymmetric capacitively coupled hydrogen and deuterium plasmas exhibit structure. The otherwise relatively continuous distribution appears to have an abrupt peak in electron density near 5 eV. This structure occurs at a higher energy in deuterium than hydrogen and there is a correlation between floating potential and the voltage at which the structure is observed in the second derivative of the I(V) characteristic. While the cause of the structure has yet to be clarified, spectroscopic observations and computer-based hydrogen models indicate that the high energy tail of the distribution is strongly modulated during the radio frequency cycle. The effect of this modulation on plasma properties and probe measurements has yet to be explored. (C) 1999 American Institute of Physics. [S0003-6951(99)00819-0].

Structure, stereochemistry and synthesis of enantiopure cyclohexenone cis-diol bacterial metabolites derived from phenols

Biotransformation of 3-substituted and 2,5-disubstituted phenols, using whole cells of P. putida UV4, yielded cyclohexenone cis-diols as single enantiomers; their structures and absolute configurations have been determined by NMR and ECD spectroscopy, X-ray crystallography, and stereochemical correlation involving a four step chemoenzymatic synthesis from the corresponding cis-dihydrodiol metabolites. An active site model has been proposed, to account for the formation of enantiopure cyclohexenone cis-diols with opposite absolute configurations.

Biosynthesis and structure of the Burkholderia cenocepacia K56-2 lipopolysaccharide core oligosaccharide:truncation of the core oligosaccharide leads to increased binding and sensitivity to polymyxin B

Burkholderia cenocepacia is an opportunistic pathogen that displays a remarkably high resistance to antimicrobial peptides. We hypothesize that high resistance to antimicrobial peptides in these bacteria is because of the barrier properties of the outer membrane. Here we report the identification of genes for the biosynthesis of the core oligosaccharide (OS) moiety of the B. cenocepacia lipopolysaccharide. We constructed a panel of isogenic mutants with truncated core OS that facilitated functional gene assignments and the elucidation of the core OS structure in the prototypic strain K56-2. The core OS structure consists of three heptoses in the inner core region, 3-deoxy-d-manno-octulosonic acid, d-glycero-d-talo-octulosonic acid, and 4-amino-4-deoxy-l-arabinose linked to d-glycero-d-talo-octulosonic acid. Also, glucose is linked to heptose I, whereas heptose II carries a second glucose and a terminal heptose, which is the site of attachment of the O antigen. We established that the level of core truncation in the mutants was proportional to their increased in vitro sensitivity to polymyxin B (PmB). Binding assays using fluorescent 5-dimethylaminonaphthalene-1-sulfonyl-labeled PmB demonstrated a correlation between sensitivity and increased binding of PmB to intact cells. Also, the mutant producing a heptoseless core OS did not survive in macrophages as compared with the parental K56-2 strain. Together, our results demonstrate that a complete core OS is required for full PmB resistance in B. cenocepacia and that resistance is due, at least in part, to the ability of B. cenocepacia to prevent binding of the peptide to the bacterial cell envelope.

Electron-impact excitation of Sc II: collision strengths and effective collision strengths for fine-structure transitions.

Accurate fine-structure atomic data for the Fe-peak elements are essential for interpreting astronomical spectra. There is a severe paucity of data available for Sc II, highlighted by the fact that no collision strengths are readily available for this ion. We present electron-impact excitation collision strengths and Maxwellian averaged effective collision strengths for Sc II. The collision strengths were calculated for all 3916 transitions amongst 89 jj levels (arising from the 3d4s, 3d2, 4s2, 3d4p, 4s4p, 3d5s, 3d4d, 3d5p, 4p2 and 3d4f configurations), resulting in a 944 coupled channel problem. The R-matrix package RMATRXII was utilized, along with the transformation code FINE and the external region code PSTGF, to calculate the collision strengths for a range of incident electron energies in the 0 to 8.3 Rydberg region. Maxwellian averaged effective collision strengths were then produced for 27 temperatures lying within the astrophysically significant range of 30 to 105 K.
The collision strengths and effective collision strengths were produced for two different target models. The purpose was to systematically examine the effect of including open 3p correlation terms into the configuration interaction expansion for the wavefunction. The first model consisted of all 36 CI terms that could be generated with the 3p core closed. The second model incorporated an additional six configurations which allowed for single-electron excitations from within the 3p core. Comparisons are made between the two models and the results of Bautista et al., obtained by private communication. It is concluded that the first model produced the most reliable set of collision and effective collision strengths for use in astrophysical and plasma applications.

Quantitative structure-toxicity relationships for chlorophenols to bioluminescent lux-marked bacteria using atom-based semi-empirical molecular-orbital descriptors

Literature data on the toxicity of chlorophenols for three luminescent bacteria (Vibrio fischeri, and the lux-marked Pseudomonas fluorescens 10586s pUCD607 and Burkholderia spp. RASC c2 (Tn4431)) have been analyzed in relation to a set of computed molecular physico-chemical properties. The quantitative structure-toxicity relationships of the compounds in each species showed marked differences when based upon semi-empirical molecular-orbital molecular and atom based properties. For mono-, di- and tri-chlorophenols multiple linear regression analysis of V. fischeri toxicity showed a good correlation with the solvent accessible surface area and the charge on the oxygen atom. This correlation successfully predicted the toxicity of the heavily chlorinated phenols, suggesting in V. fischeri only one overall mechanism is present for all chlorophenols. Good correlations were also found for RASC c2 with molecular properties, such as the surface area and the nucleophilic super-delocalizability of the oxygen. In contrast the best QSTR for P. fluorescens contained the 2nd order connectivity index and ELUMO suggesting a different, more reactive mechanism. Cross-species correlations were examined, and between V. fischeri and RASC c2 the inclusion of the minimum value of the nucleophilic susceptibility on the ring carbons produced good results. Poorer correlations were found with P. fluorescens highlighting the relative similarity of V. fischeri and RASC c2, in contrast to that of P. fluorescens.