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.

Structure and Rheology of Cationic Molecular Hydrogels of Quinuclidine Grafted Bile Salts. Influence of the Ionic Strength and Counter-Ion type

Quinuclidine grafted cationic bile salts are forming salted hydrogels. An extensive investigation of the effect of the electrolyte and counterions on the gelation has been envisaged. The special interest of the quinuclidine grafted bile salt is due to its broader experimental range of gelation to study the effect of electrolyte. Rheological features of the hydrogels are typical of enthalpic networks exhibiting a scaling law of the elastic shear modulus with the concentration (scaling exponent 2.2) modeling cellular solids in which the bending modulus is the dominant parameter. The addition of monovalent salt (NaCl) favors the formation of gels in a first range (0.00117 g cm-3 (0.02 M) < TNaCl < 0.04675 g cm-3 (0.8 M)). At larger salt concentrations, the gels become more heterogeneous with nodal zones in the micron scale. Small-angle neutron scattering experiments have been used to characterize the rigid fibers ( ≈ 68 Å) and the nodal zones. Stress sweep and creeprecovery measurements are used to relate the lack of linear viscoelastic domain to a mechanism of disentanglement of the fibers from their associations into fagots. The electrostatic interactions can be screened by addition of salt to induce a progressive evolution toward flocculation. SEM, UV absorbance, and SAXS study of the Bragg peak at large Q-values complete the investigation.

Micro-Raman imaging of Te precipitates in CdZnTe (Zn similar to 4%) crystals

Micro-Raman imaging of the distribution of Te precipitates in CdZnTe crystals in different phases is reported. For the normal phase of Te precipitates, the Raman modes appear centered around 121(A1), 141(E)/TO(CdTe) cm−1 and a weak mode around 92(E) cm−1 in CdZnTe indicating the presence of trigonal lattice of Te. Under high pressure phase, the volume of Te precipitates collapses, giving more bond energy resulting in the blueshift of the corresponding Raman bands. Also, the spatial distribution of the area ratio of 121 to 141 cm−1 Raman modes is used to quantify Te precipitates. Further, near-infrared microscopy images support these results.

trans-Dichloro-bis-(arylazoimidazole)palladium(II): Synthesis, Structure, Photoisomerization, and DFT Calculation

Reaction between PdCl2 and 1-alkyl-2-(arylazo)imidazole (RaaiR') or 1-alkyl-2-(naphthyl-alpha/beta-azo)imidazole (alpha/beta-NaiR') under reflux in ethanol has isolated complexes of compositions Pd(RaaiR')(2)Cl-2 (5, 6) and Pd(alpha/beta-NaiR')(2)Cl-2 (7, 8). The X-ray structure determination of one of the molecules, Pd(alpha-NaiBz)(2)Cl-2 (7c), has reported a trans-PdCl2 configuration, and alpha-NaiBz acts as monodentate N(imidazole) donor ligand. The spectral (IR, UV-vis, H-1 NMR) data support the structure. UV light irradiation (light source: Perkin-Elmer LS 55 spectrofluorimeter, Xenon discharge lamp, lambda = 360-396 nm) in a MeCN solution of the complexes shows E-to-Z isomerization of the coordinated azoimidazole unit. The reverse transformation, Z-to-E, is very slow with visible light irradiation. Quantum yields (phi(E-Z)) of E-to-Z isomerization are calculated, and phi is lower than that of the free ligand but comparable with those of Cd(II) and Hg(II) complexes of the same ligand. The Z-to-E isomerization is a thermally induced process. The activation energy (E-a) of Z-to-E isomerization is calculated by controlled-temperature experimentation. cis-Pd(azoimidazole)Cl-2 complexes (azomidazole acts as N(imidazole) and N(azo) Chelating ligand) do not respond upon light irradiation, which supports the idea that the presence of noncoordinated azo-N to make free azo (-N=N-) function is important to reveal photochromic activity. DFT calculation of Pd(alpha-NaiBz)(2)Cl-2 (7c) has suggested that the HOMO of the molecule is constituted of Pd (32%) and Cl (66%), and hence photo excitation may use the energy of Pd and Cl instead of that of the photofunctional -N=N-Ar motif; thus, the rate of photoisomerization and quantum yield decrease versus the free ligand values.

3D numerical simulation of avascular tumour growth: effect of hypoxic micro-environment in host tissue

A three-dimensional (3D) mathematical model of tumour growth at the avascular phase and vessel remodelling in host tissues is proposed with emphasis on the study of the interactions of tumour growth and hypoxic micro-environment in host tissues. The hybrid based model includes the continuum part, such as the distributions of oxygen and vascular endothelial growth factors (VEGFs), and the discrete part of tumour cells (TCs) and blood vessel networks. The simulation shows the dynamic process of avascular tumour growth from a few initial cells to an equilibrium state with varied vessel networks. After a phase of rapidly increasing numbers of the TCs, more and more host vessels collapse due to the stress caused by the growing tumour. In addition, the consumption of oxygen expands with the enlarged tumour region. The study also discusses the effects of certain factors on tumour growth, including the density and configuration of preexisting vessel networks and the blood oxygen content. The model enables us to examine the relationship between early tumour growth and hypoxic micro-environment in host tissues, which can be useful for further applications, such as tumour metastasis and the initialization of tumour angiogenesis.

Syntheses and characterization of trans-[NiL2(NCS)2][L = 2-(aminomethyl) pyridine], trans-[NiL02(NSC)2][L0 = 2-(2-aminoethyl)pyridine] and trans-[NiL00 2(NSC)2] [L00 = 2-(2-methylaminoethyl)pyridine] complexes: X-ray single crystal structure of trans-[NiL02(NSC)2] [L0 = 2-(2-aminoethyl)pyridine]

[NiL2(NCS)2] (1) [L = 2-(aminomethyl)pyridine], [NiL02(NCS)2] (2) [(L0) = 2-(2-aminoethyl)pyridine and [NiL00 2(NCS)2] (3) [L00 = 2-(2-methylaminoethyl)pyridine] have been synthesized from solution. All the complexes possess trans geometry as is evident from solid state UV–Vis spectral study and X-ray single crystal structure analysis of complex 2 unambiguously proves trans geometry of the species.

Stress analysis of carotid atheroma in a transient ischaemic attack patient using the MRI-based fluid-structure interaction method

Rupture of atherosclerotic plaque is a major cause of mortality. Plaque stress analysis, based on patient-specific multisequence in vivo MRI, can provide critical information for the understanding of plaque rupture and could eventually lead to plaque rupture prediction. However, the direct link between stress and plaque rupture is not fully understood. In the present study, the plaque from a patient who recently experienced a transient ischaemic attack (TIA) was studied using a fluid-structure interaction method to quantify stress distribution in the plaque region based on in vivo MR images. The results showed that wall shear stress is generally low in the artery with a slight increase at the plaque throat owing to minor luminal narrowing. The oscillatory shear index is much higher in the proximal part of the plaque. Both local wall stress concentrations and the relative stress variation distribution during a cardiac cycle indicate that the actual plaque rupture site is collocated with the highest rupture risk region in the studied patient.

Fluid-structure interaction study on human arterial plaque with patient specific geometry and boundary conditions

Atherosclerotic plaque rupture has been extensively considered as the leading cause of death in western countries. It is believed that high stresses within plaque can be an important factor on triggering the rupture of the plaque. Stress analysis in the coronary and carotid arteries with plaque have been developed by many researchers from 2D to 3-D models, from structure analysis only to the Fluid-Structure Interaction (FSI) models[1].

Stress analysis on human arterial plaques by fluid structure interactions - Multi-case study

Atherosclerotic plaque rupture has been extensively considered as the leading cause of death in the world. It is believed that high stress within plaque can be an important factor which can trigger the rupture of the plaque. High resolution multi-spectral magnetic resonance imaging (MRI) has allowed the plaque components (arterial wall, lipids, and fibrous cap) to be visualized in vivo [1]. The patient specific finite element model can be generated from the image data to perform stress analysis and provide critical information on understanding plaque rupture mechanisms [2]. The present work is to apply the procedure to a total of 14 patients (S1 ∼ S14), to study the stress distributions on carotid artery plaque reconstructed from multi-spectral magnetic resonance images, and the possible relationships between stress and plaque burdens.

Carotid arterial plaque stress analysis using fluid-structure interactive simulation based on in-vivo magnetic resonance images of four patients

The rupture of atherosclerotic plaques is known to be associated with the stresses that act on or within the arterial wall. The extreme wall tensile stress (WTS) is usually recognized as a primary trigger for the rupture of vulnerable plaque. The present study used the in-vivo high-resolution multi-spectral magnetic resonance imaging (MRI) for carotid arterial plaque morphology reconstruction. Image segmentation of different plaque components was based on the multi-spectral MRI and co-registered with different sequences for the patient. Stress analysis was performed on totally four subjects with different plaque burden by fluid-structure interaction (FSI) simulations. Wall shear stress distributions are highly related to the degree of stenosis, while the level of its magnitude is much lower than the WTS in the fibrous cap. WTS is higher in the luminal wall and lower at the outer wall, with the lowest stress at the lipid region. Local stress concentrations are well confined in the thinner fibrous cap region, and usually locating in the plaque shoulder; the introduction of relative stress variation during a cycle in the fibrous cap can be a potential indicator for plaque fatigue process in the thin fibrous cap. According to stress analysis of the four subjects, a risk assessment in terms of mechanical factors could be made, which may be helpful in clinical practice. However, more subjects with patient specific analysis are desirable for plaque-stability study.

Local structure of Sr2FeMoxW1-xO6 double perovskites across the composition-driven metal to insulator transition

Sr2FeMoO6 oxides exhibit a half-metallic ferromagnetic (HM-FM) ground state and peculiar magnetic and magnetotransport properties, which are interesting for applications in the emerging field of spintronics and attractive for fundamental research in the field of heavily correlated electron systems. Sr2FeWO6 is an insulator with an antiferromagnetic (I-AFM) ground state. The solid solutions Sr2FeMoxW1-xO6 also have peculiar properties-W doping enhances chemical order which allows stabilization of the HM-FM state; as the W content exceeds a certain value a metal to insulator transition (MIT) occurs. The role of W in determining the physical properties of Sr2FeMoxW1-xO6 systems has been a matter of intense investigation. This work deals with the problem of the structural and electronic changes related to the MIT from a local perspective by means of x-ray absorption spectroscopy (XAS). This technique allows one to probe in detail the local structure and electronic modifications around selected absorber ions (W, Mo, Fe and Sr in our case). The results of XAS analysis in the whole composition range (0 <= x <= 1), in the near edge (XANES) and extended (EXAFS) regions, demonstrate an abrupt change of the local structure around the Fe and Mo sites at the critical composition, x(c). This change represents the microstructural counterpart associated with the MIT. Conversely, the local structure and electronic configuration of W ions remain unaltered in the whole composition range, suggesting indirect participation of W in the MIT.

Similarities in the Genomic Sequence and Coat Protein-Structure of Plant Viruses

The genomic sequences of several RNA plant viruses including cucumber mosaic virus, brome mosaic virus, alfalfa mosaic virus and tobacco mosaic virus have become available recently. The former two viruses are icosahedral while the latter two are bullet and rod shaped, respectively in particle morphology. The non-structural 3a proteins of cucumber mosaic virus and brome mosaic virus have an amino acid sequence homology of 35% and hence are evolutionarily related. In contrast, the coat proteins exhibit little homology, although the circular dichroism spectrum of these viruses are similar. The non-coding regions of the genome also exhibit variable but extensive homology. Comparison of the brome mosaic virus and alfalfa mosaic virus sequences reveals that they are probably related although with a much larger evolutionary distance. The polypeptide folds of the coat protein of three biologically distinct isometric plant viruses, tomato bushy stunt virus, southern bean mosaic virus and satellite tobacco necrosis virus have been shown to display a striking resemblance. All of them consist of a topologically similar 8-standard β-barrel. The implications of these studies to the understanding of the evolution of plant viruses will be discussed.

Structure of nonactin-calcium perchlorate, C40H64O12.Ca(ClO4)2, and a comparative study of metal-nonactin complexes

M r= 975.9, orthorhombic, Pnna, a = 20.262 (3), b= 15.717 (2), c= 15.038 (1)A, V= 4788.97 A 3, z = 4, D x = 1.35 Mg m -3, Cu Kct radiation, 2 = 1.5418 A, /t = 2.79 mm -1, F(000) -= 2072, T = 293 K, R = 0.08, 3335 observed reflections. The molecular structure and the crystal packing are similar to those observed in the nonactin complexes of sodium thiocyanate and potassium thiocyanate. The eight metal-O distances are nearly the same in the potassium complex whereas the four distances involving carbonyl O atoms are shorter than the remaining four involving the tetrahydrofuran-ring O atoms in the Na and the Ca complexes. This observation can be explained in terms of the small ionic radii of Na + and Ca 2+, and leads to a plausible structural rationale for the stronger affinity of nonactin for K + than for the other two metal ions.

Structure of the Monosodium Salt of D-Glucose 6-Hydrogenphosphate

Na+.C6HI209 P-, Mr=282.1, monoclinic, e2~, a=5-762(1), b=7.163(2), c=12.313(1)A, fl= 99.97 (1) °, U= 500.5 A 3, Z= 2, D m = 1.86, D x = 1.87 Mg m -s, Cu Ka, 2 = 1.5418 A, /a = 3-3 mm -1, F(000) = 292, T= 300 K, final R for 922 observed reflections is 0-042. The phosphate ester bond, P-O(6), is 1.575 (5)A, slightly shorter than the P~O bond in monopotassium phosphoenolpyruvate [1.612 (6) A] [Hosur & Viswamitra (1981). Acta Cryst. B37, 839-843]. The pyranose sugar ring takes a 4C 1 chair conformation. The conformation about the exocyclic C(5)-C(6) bond is gauche-trans. The endocyclic C-O bonds in the glucose ring are nearly equal with C(5)-O(5) = 1.435 (8) and C(1)-O(5) = 1.436 (9) A. The sodium ion has seven near neighbours within a distance of 2.9 A. The crystal structure is stabilized by hydrogen bonds between the O atoms of symmetryrelated molecules.

Conformational flexibility in androgenic steroids - the structure of a new form of (+)-17-beta-hydroxy-19-nor-4-androsten-3-one (19-nortestosterone), c18h26o2

The structure and conformation of a second crystalline modification of 19-nortestosterone has been determined by X-ray methods. M r = 274, monoclinic P2 l, a=9.755(2), b= 11.467(3), c= 14.196(3)/L fl=101.07(2) ° , V=1558.4 (8) A 3, Z=4, Ox= I. 168 g cm -3, Mo Ka, 2 = 0.7107 ,/k, ~ = 0.80 cm -l, F(000) = 600, T= 300 K. R = 0.060 for 2158 observed reflections. The two molecules in the asymmetric unit show significant differences in the A-ring conformation from that of the previously reported form of the title compound [Precigoux, Busetta, Courseille & Hospital (1975). Acta Cryst. B31, 1527-1532]. The l a,2fl-half-chair conformation of the A ring increases its conformational freedom compared with testosterone.