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, function and intracellular dynamics of alphavirus replication complexes

Intracellular membrane alterations are hallmarks of positive-sense RNA (+RNA) virus replication. Strong evidence indicates that within these exotic compartments, viral replicase proteins engage in RNA genome replication and transcription. To date, fundamental questions such as the origin of altered membranes, mechanisms of membrane deformation and topological distribution and function of viral components, are still waiting for comprehensive answers. This study addressed some of the above mentioned questions for the membrane alterations induced during Semliki Forest virus (SFV) infection of mammalian cells. With the aid of electron and fluorescence microscopy coupled with radioactive labelling and immuno-cytochemistry techniques, our group and others showed that few hours after infection the four non structural proteins (nsP1-4) and newly synthesized RNAs of SFV colocalized in close proximity of small membrane invaginations, designated as spherules . These 50-70 nm structures were mainly detected in the perinuclear area, at the limiting membrane of modified endosomes and lysosomes, named CPV-I (cytopathic vacuoles type I). More rarely, spherules were also found at the plasma membrane (PM). In the first part of this study I present the first three-dimensional reconstruction of the CPV-I and the spherules, obtained by electron tomography after chemical or cryo-fixation. Different approaches for imaging these macromolecular assemblies to obtain better structure preservation and higher resolution are presented as unpublished data. This study provides insights into spherule organization and distribution of viral components. The results of this and other experiments presented in this thesis will challenge currently accepted models for virus replication complex formation and function. In a revisitation of our previous models, the second part of this work provides the first complete description of the biogenesis of the CPV-I. The results demonstrate that these virus-induced vacuoles, where hundreds of spherules accumulate at late stages during infection, represent the final phase of a journey initiated at the PM, which apparently serves as a platform for spherule formation. From the PM spherules were internalized by an endocytic event that required the activity of the class I PI3K, caveolin-1, cellular cholesterol and functional actin-myosin network. The resulting neutral endocytic carrier vesicle delivered the spherules to the membrane of pre-existing acidic endosomes via multiple fusion events. Microtubule based transport supported the vectorial transfer of these intermediates to the pericentriolar area where further fusions generated the CPV-I. A signal for spherule internalization was identified in one of the replicase proteins, nsP3. Infections of cells with viruses harbouring a deletion in a highly phosphorylated region of nsP3 did not result in the formation of CPV-Is. Instead, thousands of spherules remained at the PM throughout the infection cycle. Finally, the role of the replicase protein nsP2 during viral RNA replication and transcription was investigated. Three enzymatic activities, protease, NTPase and RNA-triphosphatase were studied with the aid of temperature sensitive mutants in vitro and, when possible, in vivo. The results highlighted the interplay of the different nsP2 functions during different steps of RNA replication and sub-genomic promoter regulation, and suggest that the protein could have different activities when participating in the replication complex or as a free enzyme.

Local structure in LaMnO3 and CaMnO3 perovskites: A quantitative structural refinement of mn K-edge XANES data

Hole-doped perovskites such as La1-xCaxMnO3 present special magnetic and magnetotransport properties, and it is commonly accepted that the local atomic structure around Mn ions plays a crucial role in determining these peculiar features. Therefore experimental techniques directly probing the local atomic structure, like x-ray absorption spectroscopy (XAS), have been widely exploited to deeply understand the physics of these compounds. Quantitative XAS analysis usually concerns the extended region [extended x-ray absorption fine structure (EXAFS)] of the absorption spectra. The near-edge region [x-ray absorption near-edge spectroscopy (XANES)] of XAS spectra can provide detailed complementary information on the electronic structure and local atomic topology around the absorber. However, the complexity of the XANES analysis usually prevents a quantitative understanding of the data. This work exploits the recently developed MXAN code to achieve a quantitative structural refinement of the Mn K-edge XANES of LaMnO3 and CaMnO3 compounds; they are the end compounds of the doped manganite series LaxCa1-xMnO3. The results derived from the EXAFS and XANES analyses are in good agreement, demonstrating that a quantitative picture of the local structure can be obtained from XANES in these crystalline compounds. Moreover, the quantitative XANES analysis provides topological information not directly achievable from EXAFS data analysis. This work demonstrates that combining the analysis of extended and near-edge regions of Mn K-edge XAS spectra could provide a complete and accurate description of Mn local atomic environment in these compounds.

Synthesis, structure, transformation studies and catalytic properties of open-framework cadmium thiosulfate compounds

Five new thiosulfate based inorganic-organic hybrid open-framework compounds have been synthesized employing mild reaction conditions. Of the five compounds, [Na-2(H2O)(8)][Cd(C10H8N2)( S2O3)(2)]center dot 2H(2)O, I and [Cd-2(C10H8N2)(2)(HS2O3)(2)(S2O3)(2)][(C10H9N2)(2)(C10H8N2)(2)]center dot 8H(2)O, II have one-dimensional (1D) structures and [Cd(C10H8N2)(H2O)(2)(S2O3)]center dot 2H(2)O, III, [Cd-2(C10H8N2)(3)(S2O3)(2)], IV and [Cd-2(C10H8N2)(2.5)(S2O3)(2)], V have three- dimensional (3D) structures. The 1D structures are somewhat related, formed by the bonding between tetrahedral Cd centers (CdN2S2) and 4,4'-bipyridine (bpy) units. The inter-chain spaces are occupied by the hanging thiosulfate units in both the cases along with Na(H2O)(6) chains in I and free bpy units in II. The three 3D structures have one-dimensional cadmium thiosulfate chains linked by bpy units. Interpenetration has been observed in all the 3D structures. The 3D structures appear to be related and can be derived from fgs net. Transformation studies on the 1D compound, [Na-2(H2O)(8)][Cd(C10H8N2)(S2O3)(2)]center dot 2H(2)O, I, indicated a facile formation of [Cd(C10H8N2)(H2O)(2)(S2O3)]center dot 2H(2)O, III. Prolonged heating of I gave rise to a 3D cadmium sulfate phase, [Cd-2(C10H8N2)(2)(H2O)(3)(SO4)(2)]center dot 2H(2)O, VI. Compound VI has one-dimensional cadmium sulfate chains formed by six-membered rings connected by bpy units to form a 3D structure, which appears to resemble the topological arrangement of III. Transformation studies of III indicates the formation of IV and V, and at a higher temperature a new 3D cadmium sulfate, [Cd(C10H8N2)(SO4)], VII. Compound VII has a 4 x 4 grid cadmium sulfate layers pillared by bpy units. All the compounds were characterized by PXRD, TGA, IR and UV-visible studies. Preliminary studies on the possible use of the 3D compounds (III-VII) in heterogeneous cyanosilylation of imines appear to be promising.

Topological defects in crystals: A density-wave theory

A new approach for describing dislocations and other topological defects in crystals, based on the density wave theory of Ramakrishnan and Yussouff is presented. Quantitative calculations are discussed in brief for the order parameter profiles, the atomic configuration and the free energy of a screw dislocation with Burgers vector b = (a/2, a/2,a/2 ) in a bcc solid. Our results for the free energy of the dislocation in a crystal of sizeR, when expressed as (λb 2/4π) ln (αR/|b|) whereλ is the shear elastic constant, yield, for example, the valueα ⋍ 1·85 for sodium at its freezing temperature (371°K). The density distribution in the presence of the dislocation shows that the dislocation core has a columnar character. To our knowledge, this study represents the first calculation of dislocation structure, including the core, within the framework of an order parameter theory incorporating thermal effects.

Charge Density Analysis of a Pentaborate Ion in an Ammonium Borate: Toward the Understanding of Topological Features in Borate Minerals

Structural and charge density distribution studies have been carried out on a single crystal data of an ammonium borate, [C(10)H(26)N(4)][B(5)O(6)(OH)(4)](2), synthesized by solvothermal method. Further, the experimentally observed geometry is used for the theoretical charge density calculations using the B3LYP/6-31G** level of theory, and the results are compared with the experimental values. Topological analysis of charge density based on the Atoms in Molecules approach for B-O bonds exhibit mixed covalent/ionic character. Detailed analysis of the hydrogen bonds in the crystal structure in the ammonium borate provides insights into the understanding of the reaction pathways that net atomic charges and electrostatic potential isosurfaces also give additional such systems. could result in the formation of borate minerals. The input to evaluate chemical and physical properties in such systems.

Use of Vertex Index in Structure-Activity Analysis and Design of Molecules

The last few decades have witnessed application of graph theory and topological indices derived from molecular graph in structure-activity analysis. Such applications are based on regression and various multivariate analyses. Most of the topological indices are computed for the whole molecule and used as descriptors for explaining properties/activities of chemical compounds. However, some substructural descriptors in the form of topological distance based vertex indices have been found to be useful in identifying activity related substructures and in predicting pharmacological and toxicological activities of bioactive compounds. Another important aspect of drug discovery e. g. designing novel pharmaceutical candidates could also be done from the distance distribution associated with such vertex indices. In this article, we will review the development and applications of this approach both in activity prediction as well as in designing novel compounds.

Interaction Signatures Stabilizing the NAD(P)-Binding Rossmann Fold: A Structure Network Approach

The fidelity of the folding pathways being encoded in the amino acid sequence is met with challenge in instances where proteins with no sequence homology, performing different functions and no apparent evolutionary linkage, adopt a similar fold. The problem stated otherwise is that a limited fold space is available to a repertoire of diverse sequences. The key question is what factors lead to the formation of a fold from diverse sequences. Here, with the NAD(P)-binding Rossmann fold domains as a case study and using the concepts of network theory, we have unveiled the consensus structural features that drive the formation of this fold. We have proposed a graph theoretic formalism to capture the structural details in terms of the conserved atomic interactions in global milieu, and hence extract the essential topological features from diverse sequences. A unified mathematical representation of the different structures together with a judicious concoction of several network parameters enabled us to probe into the structural features driving the adoption of the NAD(P)-binding Rossmann fold. The atomic interactions at key positions seem to be better conserved in proteins, as compared to the residues participating in these interactions. We propose a ``spatial motif'' and several ``fold specific hot spots'' that form the signature structural blueprints of the NAD(P)-binding Rossmann fold domain. Excellent agreement of our data with previous experimental and theoretical studies validates the robustness and validity of the approach. Additionally, comparison of our results with statistical coupling analysis (SCA) provides further support. The methodology proposed here is general and can be applied to similar problems of interest.

Exploiting the Structure of the Constraint Graph for Efficient Points-to Analysis

Points-to analysis is a key compiler analysis. Several memory related optimizations use points-to information to improve their effectiveness. Points-to analysis is performed by building a constraint graph of pointer variables and dynamically updating it to propagate more and more points-to information across its subset edges. So far, the structure of the constraint graph has been only trivially exploited for efficient propagation of information, e.g., in identifying cyclic components or to propagate information in topological order. We perform a careful study of its structure and propose a new inclusion-based flow-insensitive context-sensitive points-to analysis algorithm based on the notion of dominant pointers. We also propose a new kind of pointer-equivalence based on dominant pointers which provides significantly more opportunities for reducing the number of pointers tracked during the analysis. Based on this hitherto unexplored form of pointer-equivalence, we develop a new context-sensitive flow-insensitive points-to analysis algorithm which uses incremental dominator update to efficiently compute points-to information. Using a large suite of programs consisting of SPEC 2000 benchmarks and five large open source programs we show that our points-to analysis is 88% faster than BDD-based Lazy Cycle Detection and 2x faster than Deep Propagation. We argue that our approach of detecting dominator-based pointer-equivalence is a key to improve points-to analysis efficiency.

Bismuth carboxylates with Brucite- and Fluorite-related structures: synthesis structure and properties

Three new compounds of bismuth, C4N2H10]center dotBi(C7H4NO4)(C7H3NO4)]center dot H2O, I, Bi(C5H3N2O4) (C5H2N2O4)], II, and Bi(mu(2)-OH)(C7H3NO4)], III, have been prepared by the reaction between bismuth nitrate and heterocyclic aromatic dicarboxylic acids, 2,6-pyridinedicarboxylic acid, 4,5-imidazoledicarboxylic acid, and 3,4-pyridinedicarboxylic acid, respectively, under hydrothermal conditions. The structures of all the compounds have linkages between Bi2O2 and the corresponding dicarboxylate forming a simple molecular unit in I, a bilayer arrangement in II, and a three-dimensional extended structure in III. The topological arrangement of the nodal building units in the structures indicates that a brucite-related layer (II) and fluorite-related arrangement (III) can be realized in these structures. By utilizing the secondary interactions, one can correlate the structure of III to a Kagome-related one. The observation of such classical inorganic related structures in the bismuth carboxylates is noteworthy. Lewis acid catalytic studies on the formation of ketal suggest the possible participatory role of the lone pair of electrons. All the compounds are characterized employing elemental analysis, IR, UV-vis, and thermal studies.

An automated approach to network features of protein structure ensembles

Network theory applied to protein structures provides insights into numerous problems of biological relevance. The explosion in structural data available from PDB and simulations establishes a need to introduce a standalone-efficient program that assembles network concepts/parameters under one hood in an automated manner. Herein, we discuss the development/application of an exhaustive, user-friendly, standalone program package named PSN-Ensemble, which can handle structural ensembles generated through molecular dynamics (MD) simulation/NMR studies or from multiple X-ray structures. The novelty in network construction lies in the explicit consideration of side-chain interactions among amino acids. The program evaluates network parameters dealing with topological organization and long-range allosteric communication. The introduction of a flexible weighing scheme in terms of residue pairwise cross-correlation/interaction energy in PSN-Ensemble brings in dynamical/chemical knowledge into the network representation. Also, the results are mapped on a graphical display of the structure, allowing an easy access of network analysis to a general biological community. The potential of PSN-Ensemble toward examining structural ensemble is exemplified using MD trajectories of an ubiquitin-conjugating enzyme (UbcH5b). Furthermore, insights derived from network parameters evaluated using PSN-Ensemble for single-static structures of active/inactive states of 2-adrenergic receptor and the ternary tRNA complexes of tyrosyl tRNA synthetases (from organisms across kingdoms) are discussed. PSN-Ensemble is freely available from http://vishgraph.mbu.iisc.ernet.in/PSN-Ensemble/psn_index.html.

Topological analysis of electron density and the electrostatic properties of isoniazid: an experimental and theoretical study

Isoniazid (isonicotinohydrazide) is an important first-line antitubercular drug that targets the InhA enzyme which synthesizes the critical component of the mycobacterial cell wall. An experimental charge-density analysis of isoniazid has been performed to understand its structural and electronic properties in the solid state. A high-resolution single-crystal X-ray intensity data has been collected at 90 K. An aspherical multipole refinement was carried out to explore the topological and electrostatic properties of the isoniazid molecule. The experimental results were compared with the theoretical charge-density calculations performed using CRYSTAL09 with the B3LYP/6-31G** method. A topological analysis of the electron density reveals that the Laplacian of electron density of the N-N bond is significantly less negative, which indicates that the charges at the b.c.p. (bond-critical point) of the bond are least accumulated, and so the bond is considered to be weak. As expected, a strong negative electrostatic potential region is present in the vicinity of the O1, N1 and N3 atoms, which are the reactive locations of the molecule. The C-H center dot center dot center dot N, C-H center dot center dot center dot O and N-H center dot center dot center dot N types of intermolecular hydrogen-bonding interactions stabilize the crystal structure. The topological analysis of the electron density on hydrogen bonding shows the strength of intermolecular interactions.

Confining Dirac electrons on a topological insulator surface using potentials and a magnetic field

We study the effects of extended and localized potentials and a magnetic field on the Dirac electrons residing at the surface of a three-dimensional topological insulator like Bi2Se3. We use a lattice model to numerically study the various states; we show how the potentials can be chosen in a way which effectively avoids the problem of fermion doubling on a lattice. We show that extended potentials of different shapes can give rise to states which propagate freely along the potential but decay exponentially away from it. For an infinitely long potential barrier, the dispersion and spin structure of these states are unusual and these can be varied continuously by changing the barrier strength. In the presence of a magnetic field applied perpendicular to the surface, these states become separated from the gapless surface states by a gap, thereby giving rise to a quasi-one-dimensional system. Similarly, a magnetic field along with a localized potential can give rise to exponentially localized states which are separated from the surface states by a gap and thereby form a zero-dimensional system. Finally, we show that a long barrier and an impurity potential can produce bound states which are localized at the impurity, and an ``L''-shaped potential can have both bound states at the corner of the L and extended states which travel along the arms of the potential. Our work opens the way to constructing wave guides for Dirac electrons.

Impact of Stone-Wales and lattice vacancy defects on the electro-thermal transport of the free standing structure of metallic ZGNR

We report the effect of topological as well as lattice vacancy defects on the electro-thermal transport properties of the metallic zigzag graphene nano ribbons at their ballistic limit. We employ the density function theory-Non equilibrium green's function combination to calculate the transmission details. We then present an elaborated study considering the variation in the electrical current and the heat current transport with the change in temperature as well as the voltage gradient across the nano ribbons. The comparative analysis shows, that in the case of topological defects, such as the Stone-Wales defect, the electrical current transport is minimum. Besides, for the voltage gradient of 0.5 Volt and the temperature gradient of 300 K, the heat current transport reduces by similar to 62 % and similar to 50% for the cases of Stones-Wales defect and lattice vacancy defect respectively, compared to that of the perfect one.

Pressure-induced phase transition in Bi2Se3 at 3 GPa: electronic topological transition or not?

In recent years, a low pressure transition around P similar to 3 GPa exhibited by the A(2)B(3)-type 3D topological insulators is attributed to an electronic topological transition (ETT) for which there is no direct evidence either from theory or experiments. We address this phase transition and other transitions at higher pressure in bismuth selenide (Bi2Se3) using Raman spectroscopy at pressure up to 26.2 GPa. We see clear Raman signatures of an isostructural phase transition at P similar to 2.4 GPa followed by structural transitions at similar to 10 GPa and 16 GPa. First-principles calculations reveal anomalously sharp changes in the structural parameters like the internal angle of the rhombohedral unit cell with a minimum in the c/a ratio near P similar to 3 GPa. While our calculations reveal the associated anomalies in vibrational frequencies and electronic bandgap, the calculated Z(2) invariant and Dirac conical surface electronic structure remain unchanged, showing that there is no change in the electronic topology at the lowest pressure transition.