A linear programming approach for estimating the structure of a sparse linear genetic network from transcript profiling data

Background: A genetic network can be represented as a directed graph in which a node corresponds to a gene and a directed edge specifies the direction of influence of one gene on another. The reconstruction of such networks from transcript profiling data remains an important yet challenging endeavor. A transcript profile specifies the abundances of many genes in a biological sample of interest. Prevailing strategies for learning the structure of a genetic network from high-dimensional transcript profiling data assume sparsity and linearity. Many methods consider relatively small directed graphs, inferring graphs with up to a few hundred nodes. This work examines large undirected graphs representations of genetic networks, graphs with many thousands of nodes where an undirected edge between two nodes does not indicate the direction of influence, and the problem of estimating the structure of such a sparse linear genetic network (SLGN) from transcript profiling data. Results: The structure learning task is cast as a sparse linear regression problem which is then posed as a LASSO (l1-constrained fitting) problem and solved finally by formulating a Linear Program (LP). A bound on the Generalization Error of this approach is given in terms of the Leave-One-Out Error. The accuracy and utility of LP-SLGNs is assessed quantitatively and qualitatively using simulated and real data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) initiative provides gold standard data sets and evaluation metrics that enable and facilitate the comparison of algorithms for deducing the structure of networks. The structures of LP-SLGNs estimated from the INSILICO1, INSILICO2 and INSILICO3 simulated DREAM2 data sets are comparable to those proposed by the first and/or second ranked teams in the DREAM2 competition. The structures of LP-SLGNs estimated from two published Saccharomyces cerevisae cell cycle transcript profiling data sets capture known regulatory associations. In each S. cerevisiae LP-SLGN, the number of nodes with a particular degree follows an approximate power law suggesting that its degree distributions is similar to that observed in real-world networks. Inspection of these LP-SLGNs suggests biological hypotheses amenable to experimental verification. Conclusion: A statistically robust and computationally efficient LP-based method for estimating the topology of a large sparse undirected graph from high-dimensional data yields representations of genetic networks that are biologically plausible and useful abstractions of the structures of real genetic networks. Analysis of the statistical and topological properties of learned LP-SLGNs may have practical value; for example, genes with high random walk betweenness, a measure of the centrality of a node in a graph, are good candidates for intervention studies and hence integrated computational – experimental investigations designed to infer more realistic and sophisticated probabilistic directed graphical model representations of genetic networks. The LP-based solutions of the sparse linear regression problem described here may provide a method for learning the structure of transcription factor networks from transcript profiling and transcription factor binding motif data.

Morphological structure of bismuth-doped n-type amorphous germanium sulphide semiconductors

Recent observation of n-type conduction in amorphous Ge20Ss_xBix at large bismuth concentrations (x = 11), which otherwise shows p-type conduction, has aroused considerable interest in the international scientific community [1]. The mechanism of such impurity incorporation in a germanium chalcogenide glass is not understood and is a topic of current interest. In our recent publications [2-10] we have brought to light some hitherto unknown and interesting features of bismuth dopants in chalcogen-rich Ge-X (X -- S, Se) glassy compositions. In this communication we present our new results of investigations on vitreous semiconductors Ge20S80 Bi using electron microscopy, electron diffraction of as-prepared and annealed/pressure quenched compositions. Our results provide conclusive support to the formation of composite clusters containing all the three elements, germanium, sulphur and bismuth, which crystallize in simpler stoichiometric compounds Bi2S3 and GeS2.

Structure of N-[bis(methylthio)methylene]cinnamamide, C12H13NOS2

M r = 251.34, monoclinic, P21/n, a =14.626 (3), b= 7.144 (1), c= 11.996 (2)\AA \betat=90.03 (2) °, V= 1253.4 (6) \AA 3, Z = 4, Dm= 1.326 (3),Dx=1.331(3)gcm -3, MoKat, \lambda = 0.7107 )\AA , \mu=3.51 cm -3, F(000) = 528.0, T-- 293 K, R -- 3.5% for1455 significant reflections. Of particular interest is an intramolecular attractive interaction between the sulphur and oxygen atoms with an S...O distance of 2.658 (3)\AA, in which the oxygen atom appears to actas a nucleophile.

The Crystal Structure of Pivaloyl- D- prolyl-L- prolyl- L-alanyl-N- methylamide

Pivaloyl-D-prolyl-L-prolyl-L-analyl-N-methylam~de (I), C1UH32N40c4r,y stallizes in the orthorhombic space group P21212,w ith four molecules in a unit cell of dimensions a = 9.982 (l),b = 10.183 (3), c = 20.746 (2)A . The structure has been refined to R 0.048 for 1 745 observed reflections. All the peptide bonds in the molecule are trans and both the prolyl residues are in the CY-exo-conformation. The molecule assumes a highly folded conformation in which a Type II' DL bend is followed by a Type I LL bend, both stabilised by intramolecular 4 + 1 hydrogen bonds. This conformation, which has been observed for the first time, is of interest in relation to the structure of gramicidin S.

An uncommon nucleotide conformation shown by molecular structure of deoxyuridine-5-phosphate and nucleic acid stereochemistry

CRYSTAL structure determinations of nucleic acid fragments have shown that several of the conformational features found in the monomeric building blocks are also manifested at the nucleic acid level. Stereochemical variations between thymine and uracil nucleotides are therefore of interest as they can provide a structural basis for some of the differences between the conformations of DNA and RNA. X-ray studies have so far not shown any major dissimilarities between these two nucleotide species although the sugar ring of deoxyribonucleotides is found to possess greater flexibility than that in ribonucleotides. We report here the molecular structure of deoxyuridine-5'-phosphate (dUMP-5') which is not a common monomer unit of DNAs as it is replaced by its 5-methyl analogue deoxythymidine-5'-phosphate (dTMP-5'). The investigation was undertaken to help determine whether or not this implied a fundamental difference between the geometries of these two molecules.

Solution structure of O-glycosylated C-terminal leucine zipper domain of human salivary mucin (MUC7)

Solution structures of a 23 residue glycopeptide II (KIS* RFLLYMKNLLNRIIDDMVEQ, where * denotes the glycan Gal-beta-(1-3)-alpha-GalNAc) and its deglycosylated counterpart I derived from the C-terminal leucine zipper domain of low molecular weight human salivary mucin (MUC7) were studied using CD, NMR spectroscopy and molecular modeling. The peptide I was synthesized using the Fmoc chemistry following the conventional procedure and the glycopeptide II was synthesized incorporating the O-glycosylated building block (N alpha-Fmoc-Ser-[Ac-4,-beta-D-Gal-(1,3)-Ac(2)alpha-D-GalN(3)]-OPfp) at the appropriate position in stepwise assembly of peptide chain. Solution structures of these glycosylated and nonglycosylated peptides were studied in water and in the presence of 50% of an organic cosolvent, trifluoroethanol (TFE) using circular dichroism (CD), and in 50% TFE using two-dimensional proton nuclear magnetic resonance (2D H-1 NMR) spectroscopy. CD spectra in aqueous medium indicate that the apopeptide I adapts, mostly, a beta-sheet conformation whereas the glycopeptide II assumes helical structure. This transition in the secondary structure, upon glycosylation, demonstrates that the carbohydrate moiety exerts significant effect on the peptide backbone conformation. However, in 50% TFE both the peptides show pronounced helical structure. Sequential and medium range NOEs, C alpha H chemical shift perturbations, (3)J(NH:C alpha H) couplings and deuterium exchange rates of the amide proton resonances in water containing 50% TFE indicate that the peptide I adapts alpha-helical structure from Ile2-Val21 and the glycopeptide II adapts alpha-helical structure from Ser3-Glu22. The observation of continuous stretch of helix in both the peptides as observed by both NMR and CD spectroscopy strongly suggests that the C-terminal domain of MUC7 with heptad repeats of leucines or methionine residues may be stabilized by dimeric leucine zipper motif. The results reported herein may be invaluable in understanding the aggregation (or dimerization) of MUC7 glycoprotein which would eventually have implications in determining its structure-function relationship.

Effect of laser processing parameters on the structure of ductile iron

Laser processing of structure sensitive hypereutectic ductile iron, a cast alloy employed for dynamically loaded automative components, was experimentally investigated over a wide range of process parameters: from power (0.5-2.5 kW) and scan rate (7.5-25 mm s(-1)) leading to solid state transformation, all the way through to melting followed by rapid quenching. Superfine dendritic (at 10(5) degrees C s(-1)) or feathery (at 10(4) degrees C s(-1)) ledeburite of 0.2-0.25 mu m lamellar space, gamma-austenite and carbide in the laser melted and martensite in the transformed zone or heat-affected zone were observed, depending on the process parameters. Depth of geometric profiles of laser transformed or melt zone structures, parameters such as dendrile arm spacing, volume fraction of carbide and surface hardness bear a direct relationship with the energy intensity P/UDb2, (10-100 J mm(-3)). There is a minimum energy intensity threshold for solid state transformation hardening (0.2 J mm(-3)) and similarly for the initiation of superficial melting (9 J mm(-3)) and full melting (15 J mm(-3)) in the case of ductile iron. Simulation, modeling and thermal analysis of laser processing as a three-dimensional quasi-steady moving heat source problem by a finite difference method, considering temperature dependent energy absorptivity of the material to laser radiation, thermal and physical properties (kappa, rho, c(p)) and freezing under non-equilibrium conditions employing Scheil's equation to compute the proportion of the solid enabled determination of the thermal history of the laser treated zone. This includes assessment of the peak temperature attained at the surface, temperature gradients, the freezing time and rates as well as the geometric profile of the melted, transformed or heat-affected zone. Computed geometric profiles or depth are in close agreement with the experimental data, validating the numerical scheme.

Crystal structure of nonadentate tricompartmental ligand derived from pyridine-2,6-dicarboxylic acid: Spectroscopic, electrochemical and thermal investigations of its transition metal(II) complexes

The coordinating behavior of a new dihydrazone ligand, 2,6-bis(3-methoxysalicylidene) hydrazinocarbonyl]pyridine towards manganese(II), cobalt(II), nickel(II), copper(II), zinc(II) and cadmium(II) has been described. The metal complexes were characterized by magnetic moments, conductivity measurements, spectral (IR, NMR, UV-Vis, FAB-Mass and EPR) and thermal studies. The ligand crystallizes in triclinic system, space group P-1, with alpha=98.491(10)degrees, beta=110.820(10)degrees and gamma=92.228(10)degrees. The cell dimensions are a=10.196(7)angstrom, b=10.814(7)angstrom, c=10.017(7)angstrom, Z=2 and V=1117.4(12). IR spectral studies reveal the nonadentate behavior of the ligand. All the complexes are neutral in nature and possess six-coordinate geometry around each metal center. The X-band EPR spectra of copper(II) complex at both room temperature and liquid nitrogen temperature showed unresolved broad signals with g(iso) = 2.106. Cyclic voltametric studies of copper(II) complex at different scan rates reveal that all the reaction occurring are irreversible. (C) 2011 Elsevier B.V. All rights reserved.

Long-term effects of Macroternies termites, herbivores and annual early fire on woody undergrowth community in Sudanian woodland, Burkina Faso

Termites, herbivores and fire are recognized as major guilds that structure woody plant communities in African savanna and woodland ecosystems. An understanding of their interaction is crucial to design appropriate management regimes. The aim of this study was to evaluate the long-term impacts of herbivore, fire and termite activities on regeneration of trees. Permanent experimental quadrats were established in 1992 in the Sudanian woodland of Burkina Faso subjected to grazing by livestock and annual early fire and the control. Within the treatment quadrats, an inventory of the woody undergrowth community was conducted on termitaria occupied by Macrotermes subhyalinus, extended termitosphere (within 5 m radius from the mound base) and adjacent area (beyond 5 m from the mound base). Hierarchical analysis was performed to determine significant differences in species richness, abundance and diversity indices among vegetation patches within fire and herbivory treatments. Grazed quadrats had significantly (P < 0.001) more species and stem density of woody undergrowth than non-grazed quadrats but maintained similar level of species richness and stem density of woody undergrowth on termitaria. There were not significant differences (P>0.05) in species richness and stem density between burnt and unburnt quadrats. Termitaria supported a highly diverse woody undergrowth with higher stem density than either the extended termitosphere or rest of quadrats. The density of woody undergrowth was significantly related with mature trees of selected species on termitaria (R-2 = 0.593; P<0.001) than that on the extended termitosphere (R-2 = 0.333; P<0.001) and adjacent area (R-2 = 0.197; P<0.001). It can be concluded that termites facilitate the regeneration of woody species while grazing and annual early fire play a minor role in the regeneration of woody species. The current policy that prohibits grazing should be revised to accommodate the interests of livestock herders. (C) 2014 Elsevier GmbH. All rights reserved.

Synthesis, and crystal and electronic structure of sodium metal phosphate for use as a hybrid capacitor in non-aqueous electrolyte

Energy storage devices based on sodium have been considered as an alternative to traditional lithium based systems because of the natural abundance, cost effectiveness and low environmental impact of sodium. Their synthesis, and crystal and electronic properties have been discussed, because of the importance of electronic conductivity in supercapacitors for high rate applications. The density of states of a mixed sodium transition metal phosphate (maricite, NaMn1/3Co1/3Ni1/3PO4) has been determined with the ab initio generalized gradient approximation (GGA)+Hubbard term (U) method. The computed results for the mixed maricite are compared with the band gap of the parent NaFePO4 and the electrochemical experimental results are in good agreement. A mixed sodium transition metal phosphate served as an active electrode material for a hybrid supercapacitor. The hybrid device (maricite versus carbon) in a nonaqueous electrolyte shows redox peaks in the cyclic voltammograms and asymmetric profiles in the charge-discharge curves while exhibiting a specific capacitance of 40 F g(-1) and these processes are found to be quasi-reversible. After long term cycling, the device exhibits excellent capacity retention (95%) and coulombic efficiency (92%). The presence of carbon and the nanocomposite morphology, identified through X-ray photoelectron spectroscopy (XPS) and transmission electron microscopy (TEM) studies, ensures the high rate capability while offering possibilities to develop new cathode materials for sodium hybrid devices.

X-Ray Studies On The Bilayer Structure Of Trypsin-Treated Rat-Brain Myelin

Trypsin-treated rat brain myelin was subjected to biochemical and X-ray studies. Untreated myelin gave rise to a pattern of three rings with a fundamental repeat period of 155 Angstrom consisting of two bilayers per repeat period, whereas myelin treated with trypsin showed a fundamental repeat period of 75 Angstrom with one bilayer per repeat period. The integrated raw intensity of the h=4 reflection with respect to the h=2 reflection is 0.38 for untreated myelin. The corresponding value reduced to 0.23, 0.18, 0.17 for myelin treated with 5, 10, 40 units of trypsin per mg of myelin, respectively, for 30 min at 30 degrees C. The decrease in relative raw intensity of the higher-order reflection relative to the lower-order reflection is suggestive of a disordering of the phosphate groups upon trypsin treatment or an increased mosaicity of the membrane or a combination of both these effects, However, trypsin treatment does not lead to a complete breakdown of the membrane, The integrated intensity of the h=1 reflection, though weak, is above the measurable threshold for untreated myelin, whereas the corresponding intensity is below the measurable threshold for trypsin-treated myelin, indicating a possible asymmetric to symmetric transition of the myelin bilayer structure about its centre after trypsin treatment.

X-ray Structure of Gelonin at 1.8 Å Resolution

Gelonin is a single chain ribosome inactivating protein (RIP) with potential application in the treatment of cancer and AIDS. Diffraction quality crystals grown using PEG3350, belong to the space group P2(1), with it a = 49.4 Angstrom b = 44.9 Angstrom, c = 137.4 Angstrom and beta = 98.4 degrees, and contain two molecules in the asymmetric unit. Diffraction data collected to 1.8 Angstrom resolution has a R(m) value of 7.3%. Structure of gelonin has been solved by the molecular replacement method, using ricin A chain as the search model. Crystallographic refinement using X-PLOR resulted in a model for which the r.m.s deviations from ideal bond lengths and bond angles are 0.012 Angstrom and 2.7 degrees, respectively The final R-factor is 18.4% for 39,806 reflections for which I > 1.0 sigma(I).The C-alpha atoms of the two molecules in the asymmetric unit superpose to within 0.38 Angstrom for 247 atom pairs. The overall fold of gelonin is similar to that of other RIPs such as ricin A chain and alpha-momorcharin, the r.m.s.d. for C-alpha superpositions being 1.3 and 1.4 Angstrom, respectively The-catalytic residues (Glu166, Arg169 and Tyr113) in the active site form a hydrogen bond scheme similar to that observed in other RIPs. The conformation of Tyr74 in the active site, however, is significantly different from that in alpha-momorcharin. Three well defined water molecules are located in the active site cavity and one of them, X319, superposes to within 0.2 Angstrom of a corresponding water molecule in the structure of alpha-momorcharin. Any of the three could be the substrate water molecule in the hydrolysis reaction catalysed by gelonin.Difference electron density for a N-linked sugar moiety has been observed near only one of the two potential glycosylation sites in the sequence. The amino acid at position 239 has been established as Lys by calculation of omit electron density maps.The two cysteine residues in the sequence, Cys44 and Cys50, form a disulphide bond, and are therefore not available for disulphide conjugation with antibodies. Based on the structure, the region of the molecule that is involved in intradimer interactions is suggested to be suitable for introducing a Cys residue for purposes of conjugation with an antibody to produce useful immunotoxins.

Electronic Structure of Vacancy Ordered Spinels, $GaMo_{4}S_{8}$ and $GaV_{4}S_{8}$, from ab Initio Calculations

We report ab initio calculations for the band dispersions and total as well as partial densities of states for vacancy ordered, clustered spinels, GaMo4S8 and GaV4S8. Results are presented for the high temperature cubic phase for both compounds. Additionally, we discuss results of similar calculations for GaMo4S8 in an idealized cubic structure, as well as the nonmagnetic and the ferromagnetic states of the low temperature rhombohedral structure. Comparison of these results allows us to discuss the unusual aspects of the electronic structure of this interesting class of compounds, and provide estimates of the crystal-field and exchange splitting strengths.

Structure of KNb0.5OPO4, a KTiOPO4 analog

The crystal structure of KNb0.5V0.5OPO4, a new KTiOPO4 isomorph, has been refined from powder X-ray diffraction data by Rietveld refinement. The structure is orthorhombic, space group Pna2(1), with a = 12.933(1), b = 6.4713(8), and c = 10.7273(6) Angstrom, Z = 8. There is a preferential distribution of Nb(V) and V(III) atoms in the octahedral M(1) [0.806Nb, 0.194V] and M(2) [0.194Nb, 0.806V] sites, the M(1)O-6 octahedra being more distorted than the M(2)O-6 octahedra. The results are compared with other KTiOPO4 derivatives.

Structure of the benzene center dot center dot center dot ICI complex: A UVPES and ab initio molecular orbital study

UVPES studies and ab initio and DFT computations have been done on the benzene...ICl complex; electron spectral data and computed orbital energies show that donor orbitals are stabilized and acceptor orbitals are destabilized due to complexation. Calculations predict an oblique structure for the complex in which the interacting site is a C=C bond center in the donor and iodine atom in the acceptor, in full agreement with earlier experimental reports. BSSE-corrected binding energies closely match the enthalpy of complexation reported, and the NBO analysis clearly reveals the involvement of the pi orbital of benzene and the sigma* orbital of ICl in the complex.