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.

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.

An asymmetric homodinuclear complex containing the mu-x-oxobis(mu-x-carboxylato)diruthenium(III) core: X-ray structure of [Ru20(o2ec6H4-p-OMe)3(en) (PPh3)2] (C104)" 3CHC13

Asymmetric tri-bridged diruthenium(III) complexes, [Ru2O(O(2)CR)(3)(en) (PPh(3))(2)](ClO4) (R = C6H4-p-X: X = OMe (1a), Me (1b); en=1,2-diaminoethane), were prepared and structurally characterized. Complex 1a 3CHCl(3), crystallizes in the triclinic space group P (1) over bar with a = 14.029(5), b = 14.205(5), c = 20.610(6) Angstrom, alpha= 107.26(3), beta = 101.84(3), gamma= 97.57(3)degrees, V= 3756(2) Angstrom(3) and Z = 2. The complex has an {Ru-2(mu-O)(mu-O(2)CR)(2)(2+)} core and exhibits [O4PRu(mu-O)RuPO2N2](+) coordination environments for the metal centers. The novel structural feature is the asymmetric arrangement of ligands at the terminal sites of the core which shows an Ru... Ru separation of 3.226(3) Angstrom and an Ru-O-Ru angle of 119.2(5)degrees. An intense visible band observed near 570 nm is assigned to a charge transfer transition involving the d pi-Ru(III) and p pi-mu-O Orbitals. Cyclic voltammetry of the complexes displays a reversible Ru-2(III,III) reversible arrow Ru-2(III,IV) couple near 0.8 V (versus SCE) in MeCN-0.1 M TBAP.

Hydrogen bonding in peptide helices - Analysis of two independent helices in the crystal structure of a peptide Boc-Val-Ala-Leu-Aib-Val-Ala-Phe-OMe

The crystal structure determination of the heptapeptide Boc-Val-Ala-Leu-Aib-Val-Ala-Phe-OMe reveals two peptide helices in the asymmetric unit, Crystal parameters are: space group P2(1), a = 10.356(2) Angstrom, b = 19.488(5) Angstrom, c = 23.756(6) Angstrom, beta = 102.25(2)degrees), V = 4685.4 Angstrom(3), Z = 4 and R = 5.7% for 7615 reflections [I>3 sigma(I)]. Both molecules adopt largely alpha-helical conformations with variations at the C-terminus, Helix type Is determined by analysing both 4-->1 and 5-->1 hydrogen-bond interactions and comparison with the results of analysis of protein structures. The presence of two 4-->1 hydrogen-bond interactions, besides four 5-->1 interact ions in both the conformations provides an opportunity to characterize bifurcated hydrogen bonds at high resolution, Comparison of the two helical conformations with related peptide structures suggests that distortions at the C-terminus are more facile than at the N-terminus.

Electronic structure of La1-xSrxCrO3

LaCrO3 is a wide-band-gap insulator which does not evolve to a metallic state even after hole doping. We report electronic structure of this compound and its Sr substituents investigated by photoemission and inverse photoemission spectroscopies in conjunction with various calculations. The results show that LaCrO 3 is close to the Mott-Hubbard insulating regime with a gap of about 2.8 eV. Analysis of Cr 2p core-level spectrum suggests that the intra-atomic Coulomb interaction strength and the charge-transfer energy to be 5.0 and 5.5 eV, respectively, We also estimate the intra-atomic exchange interaction strength and a crystal-field splitting of about 0.7 and 2.0 eV, respectively. Sr substitution leading to hole doping in this system decreases the charge-excitation gap, but never collapses it to give a metallic behavior. The changes in the occupied as well as unoccupied spectral features are discussed in terms of the formation of local Cr4+ configurations arising from strong electron-phonon interactions.

Axial binding of 2-methylimidazole to the tetraarylcarboxylatodiruthenium (II, III) core: X-ray crystal structure of [Ru-2(O2CC6H4-p-OMe)(4)(2-mimH)(2)](ClO4). 1.75CH(2)Cl(2).H2O

Diruthenium (II. III) complexes of the type [Ru-2(O2CAr)(4) (2-mimH)(2)](ClO4) (Ar = C6H4-p-X : X=OMe,1, X=Me, 2, 2-mimH=2-methylimidazole) have been isolated from the reaction of Ru2Cl(O2CAr)(4) with 2-mimH in CH2Cl2 followed by the addition of NaClO4. The crystal structure of 1.1.75CH(2)Cl(2).H2O has been determined. The crystal belongs to the monoclinic space group p2(1)/c with the following unit cell dimensions for the C40H40N4O16ClRu2.1.75CH(2)Cl(2).H2O (M = 1237.0) : a = 12.347(3)Angstrom, b = 17.615(5)Angstrom, c = 26.148(2)Angstrom,beta = 92.88(1)degrees. v = 5679(2)Angstrom(3). Z=4, D-c = 1.45 g cm(-3). lambda(Mo-K-alpha) = 0.7107 Angstrom, mu(Mo-K-alpha) = 8.1 cm(-1), T = 293 K, R = 0.0815 (wR(2) = 0.2118) for 5834 reflections with 1 > 2 sigma(I). The complex has a tetracarboxylatodiruthenium (II, III) core and two axially bound 2-methylimidazole ligands. The Ru-Ru bond length is 2.290(1)Angstrom. The Ru-Ru bond order is 2.5 and the complex is three-electron paramagnetic. The complex shows an irreversible Ru-2(II,III)-->Ru-2(Il,II) reduction near -0.2 V vs SCE in CH2Cl2-0. 1 MTBAP. The complexes exemplify the first adduct of the tetracarboxylatodiruthenium (II,III) core having N-donor ligands

Enantiospecific approach to the tricyclic core structure of tricycloillicinone, ialibinones, and takaneones via ring-closing metathesis reaction

Enantiospecific synthesis of the tricyclic core structure present in the biologically active natural products tricycloillicinone, ialibinones, and takaneones, starting from the readily available campholenaldehyde employing a transannular RCM reaction as the key step, has been accomplished.

Importance of amino-terminus in maintenance of oligomeric structure of cytosolic serine hydroxymethyltransferase

The role of the amino and carboxyl-terminal regions of cytosolic serine hydroxymethyltransferase (SHMT) in subunit assembly and catalysis was studied using six amino-terminal (lacking the first 6, 14, 30, 49, 58, and 75 residues) and two carboxyl-terminal (lacking the last 49 and 185 residues) deletion mutants. These mutants were constructed from a full length cDNA clone using restriction enzyme/PCR-based methods and overexpressed in Escherichia coli. The overexpressed proteins, des-(A1-K6)-SHMT and des-(A1- W14)-SHMT were present in the soluble fraction and they were purified to homogeneity. The deletion clones, for des-(A1–V30)-SHMT and des-(A1–L49)-SHMT were expressed at very low levels, whereas des-(A1–R58)-SHMT, des-(A1–G75)-SHMT, des-(Q435–F483)-SHMT and des-(L299-F483)-SHMT mutant proteins were not soluble and formed inclusion bodies. Des-(A1–K6)-SHMT and des-(A1–W14)-SHMT catalyzed both the tetrahydrofolate-dependent and tetrahydrofolate-independent reactions, generating characteristic spectral intermediates with glycine and tetrahydrofolate. The two mutants had similar kinetic parameters to that of the recombinant SHMT (rSHMT). However, at 55 °C, the des-(A1–W14)-SHMT lost almost all the activity within 5 min, while at the same temperature rSHMT and des-(A1–K6)-SHMT retained 85% and 70% activity, respectively. Thermal denaturation studies showed that des-(A1–W14)-SHMT had a lower apparent melting temperature (52°C) compared to rSHMT (56°C) and des-(A1–K6)-SHMT (55 °C), suggesting that N-terminal deletion had resulted in a decrease in the thermal stability of the enzyme. Further, urea induced inactivation of the enzymes revealed that 50% inactivation occurred at a lower urea concentration (1.2 ± 0.1 M) in the case of des-(A1–W14)-SHMT compared to rSHMT (1.8 ±0.1 M) and des-(A1–K6)-SHMT (1.7 ±0.1 M). The apoenzyme of des-(A1- W14)-SHMT was present predominantly in the dimer form, whereas the apoenzymes of rSHMT and des-(A1–K6)-SHMT were a mixture of tetramers (≈75% and ≈65%, respectively) and dimers. While, rSHMT and des-(A1–K6)-SHMT apoenzymes could be reconstituted upon the addition of pyridoxal-5'-phosphate to 96% and 94% enzyme activity, respectively, des-(A1–W14)-SHMT apoenzyme could be reconstituted only upto 22%. The percentage activity regained correlated with the appearance of visible CD at 425 nm and with the amount of enzyme present in the tetrameric form upon reconstitution as monitored by gel filtration. These results demonstrate that, in addition to the cofactor, the N-terminal arm plays an important role in stabilizing the tetrameric structure of SHMT.

Importance of the amino terminus in maintenance of oligomeric structure of sheep liver cytosolic serine hydroxymethyltransferase

The Role Of The Amino And Carboxyl-Terminal Regions Of Cytosolic Serine Hydroxymethyltransferase (SHMT) In Subunit Assembly And Catalysis Was Studied Using Sis Amino-Terminal (Lacking The First 6, 14, 30, 49, 58, And 75 Residues) And Two Carboxyl-Terminal (Lacking The Last 49 And 185 Residues) Deletion Mutants. These Mutants Were Constructed From A Full Length Cdna Clone Using Restriction Enzyme/PCR-Based Methods And Overexpressed In Escherichia Coli. The Overexpressed Proteins, Des-(A1-K6) SHMT And Des-(A1-W14)-SHMT Were Present In The Soluble Fraction And They Were Purified To Homogeneity. The Deletion Clones, For Des-(A1-V30)-SHMT And Des-(A1-L49)-SHMT Were Expressed At Very Low Levels, Whereas Des-(A1-R58)-SHMT, Des-/A1-G75)-SHMT, Des-(Q435-F483)-SHMT And Des-(L299-F483)-SHMT Mutant Proteins Were Not Soluble And Formed Inclusion Bodies. Des-(A1-K6)-SHMT And Des-(A1-W14)-SHMT Catalyzed Both The Tetrahydrofolate-Dependent And Tetrahydrofolate-Independent Reactions, Generating Characteristic Spectral Intermediates With Glycine And Tetrahydrofolate. The Two Mutants Had Similar Kinetic Parameters To That Of The Recombinant SHMT (Rshmt). However, At 55 Degrees C, The Des-(A1-W14)-SHMT Lost Almost All The Activity Within 5 Min, While At The Same Temperature Rshmt And Des-(A1-K6)-SHMT Retained 85% And 70% Activity, Respectively. Thermal Denaturation Studies Showed That Des-(A1-W14)-SHMT Had A Lower Apparent Melting Temperature (52 Degrees C) Compared To Rshmt (56 Degrees C) And Des-(A1-K6)-SHMT (55 Degrees C), Suggesting That N-Terminal Deletion Had Resulted In A Decrease In The Thermal Stability Of The Enzyme. Further Urea Induced Inactivation Of The Enzymes Revealed That 50% Inactivation Occurred At A Lower Urea Concentration (1.2+/-0.1 M) In The Case Of Des-(A1-W14)-SHMT Compared To Rshmt (1.8+/-0.1 M) And Des-(A1 -K6)-SHMT (1.7+/-0.1 M). The Apoenzyme Of Des-/A1-K6)-SHMT Was Present Predominantly In The Dimer Form, Whereas The Apoenzymes Of Rshmt And Des-(A1-K6)-SHMT Were A Mixture Of Tetramers (Approximate To 75% And Approximate To 65%, Respectively) And Dimers. While, Rshmt And Des-(A1-K6)-SHMT Apoenzymes Could Be Reconstituted Upon The Addition Of Pyridoxal-5'-Phosphate To 96% And 94% Enzyme Activity, Respectively Des-(A1-W14)-SHMT Apoenzyme Could Be Reconstituted Only Upto 22%. The Percentage Activity Refined Correlated With The Appearance Of Visible CD At 425 Nm And With The Amount Of Enzyme Present In The Tetrameric Form Upon Reconstitution As Monitored By Gel Filtration. These Results Demonstrate That, In Addition To The Cofactor, The N-Terminal Arm Plays An Important Role In Stabilizing The Tetrameric Structure Of SHMT.

The crystal structure of eulytite Pb3BiV3O12

The crystal structure of Pb3BiV3O12 was solved using single-crystal X-ray diffraction technique. The compound crystallizes in the cubic system View the MathML source (No. 220) with eulytite structure with a = 10.7490(7) Å, V = 1241.95(14) Å3 and Z = 4. The final R1 value of 0.0198 (wR2=0.0384) was achieved for 359 independent reflections during the structure refinement. The Pb2+ and Bi3+ cations occupy the special position (16c) while the oxygen anions occupy the general position (48e) in the crystal structure. Unlike many other eulytite compounds, all the crystallographic positions are fully occupied. The structure consists of edge-shared Pb/Bi octahedra linked at the corners to independent [VO4]3− tetrahedra units, generating a eulytite-type network in the crystal lattice.

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.