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.

Essential dynamics and sidechain hydrogen bond cluster studies on eosinophil cationic protein

Eosinophil Cationic Protein (ECP) is a member of RNase A superfamily which carries out the obligatory catalytic role of cleaving RNA. It is involved in a variety of biological functions. Molecular dynamics simulations followed by essential dynamics analysis on this protein are carried out with the goal of gaining insights into the dynamical properties at atomic level. The top essential modes contribute to subspaces and to the transition phase. Further, the sidechain-sidechain/sidechain-mainchain hydrogen bond clusters are analyzed in the top modes, and compared with those of crystal structure. The role of residues identified by these methods is discussed in the context of concerted motion, structure and stability of the protein.

Kinematic analysis and design of articulated manipulators with joint motion constraints

For an articulated manipulator with joint rotation constraints, we show that the maximum workspace is not necessarily obtained for equal link lengths but is also determined by the range and mean positions of the joint motions. We present expressions for sectional area, workspace volume, overlap volume and work area in terms of link ratios, mean positions and ranges of joint motion. We present a numerical procedure to obtain the maximum rectangular area that can be embedded in the workspace of an articulated manipulator with joint motion constraints. We demonstrate the use of analytical expressions and the numerical plots in the kinematic design of an articulated manipulator with joint rotation constraints.

Factors influencing stick-slip motion: effect of hardness, crystal structure and surface texture

In the present investigation, efforts were made to study the different frictional responses of materials with varying crystal structure and hardness during sliding against a relatively harder material of different surface textures and roughness. In the experiments, pins were made of pure metals and alloys with significantly different hardness values. Pure metals were selected based on different class of crystal structures, such as face centered cubic (FCC), body centered cubic (BCC), body centered tetragonal (BCT) and hexagonal close packed (HCP) structures. The surface textures with varying roughness were generated on the counterpart plate which was made of H-11 die steel. The experiments were conducted under dry and lubricated conditions using an inclined pin-on-plate sliding tester for various normal loads at ambient environment. In the experiments, it was found that the coefficient of friction is controlled by the surface texture of the harder mating surfaces. Further, two kinds of frictional response, namely steady-state and stick-slip, were observed during sliding. More specifically, stead-state frictional response was observed for the FCC metals, alloys and materials with higher hardness. Stick-slip frictional response was observed for the metals which have limited number of slip systems such as BCT and HCP. In addition, the stick-slip frictional response was dependent on the normal load, lubrication, hardness and surface texture of the counterpart material. However, for a given kind of surface texture, the roughness of the surface affects neither the average coefficient of friction nor the amplitude of stick-slip oscillation significantly.

Divide and Conquer: Efficient Large-Scale Structure from Motion Using Graph Partitioning

Despite significant advances in recent years, structure-from-motion (SfM) pipelines suffer from two important drawbacks. Apart from requiring significant computational power to solve the large-scale computations involved, such pipelines sometimes fail to correctly reconstruct when the accumulated error in incremental reconstruction is large or when the number of 3D to 2D correspondences are insufficient. In this paper we present a novel approach to mitigate the above-mentioned drawbacks. Using an image match graph based on matching features we partition the image data set into smaller sets or components which are reconstructed independently. Following such reconstructions we utilise the available epipolar relationships that connect images across components to correctly align the individual reconstructions in a global frame of reference. This results in both a significant speed up of at least one order of magnitude and also mitigates the problems of reconstruction failures with a marginal loss in accuracy. The effectiveness of our approach is demonstrated on some large-scale real world data sets.

Structure of tri-1-naphthylborane-benzene (1/1) complex, c30h21b.c6h6

M r= 470.46, rhombohedral, R3, a =8.710(4)A, a=91.10(3) o, V= 660.4 (9) A 3, Z= 1,D m= 1.170 (flotation in KI solution), D x=1.183 Mg m -a, Mo Kct, 2 = 0.7107/~,, /t =0.033 mm -1, F(000) - 248.0, T= 293 K, R -- 4.6%(481 unique reflections). The molecule has C a symmetry and is propeller shaped, the angle of twist about the B-C bond being 41.5 (7) °. The space group being chiral, this is yet another example of spontaneous resolution. The results of a thermal-motion analysis are discussed.

Structure, energetics and diffusion properties of isomers of trimethyl benzene in beta zeolite: Uptake and Monte Carlo simulation study

A Monte Carlo study along with experimental uptake measurements of 1,2,3-trimethyl benzene, 1,2,4-trimethyl benzene and 1,3,5-trimethyl benzene (TMB) in beta zeolite is reported. The TraPPE potential has been employed for hydrocarbon interaction and harmonic potential of Demontis for modeling framework of the zeolite. Structure, energetics and dynamics of TMB in zeolite beta from Monte Carlo runs reveal interesting information about the diameter, properties of these isomers on confinement. Of the three isomers, 135TMB is supposed to have the largest diameter. It is seen TraPPE with Demontis potential predicts a restricted motion of 135TMB in the channels of zeolite beta.Experimentally, 135TMB has the highest transport diffusivity whereas MID results suggest this has the lowest self diffusivity. (C) 2009 Elsevier Inc. Ail rights reserved.

Neutron-diffraction study of structure and conformation of nucleotide uridine-5 phosphate (disodium salt)

Acta Crystallographica Section A: Foundations of Crystallography covers theoretical and fundamental aspects of the structure of matter. The journal is the prime forum for research in diffraction physics and the theory of crystallographic structure determination by diffraction methods using X-rays, neutrons and electrons. The structures include periodic and aperiodic crystals, and non-periodic disordered materials, and the corresponding Bragg, satellite and diffuse scattering, thermal motion and symmetry aspects. Spatial resolutions range from the subatomic domain in charge-density studies to nanodimensional imperfections such as dislocations and twin walls. The chemistry encompasses metals, alloys, and inorganic, organic and biological materials. Structure prediction and properties such as the theory of phase transformations are also covered.

Experiments on the fine structure of turbulence

Investigations have been carried out of some aspects of the fine-scale structure of turbulence in grid flows, in boundary layers in a zero pressure gradient and in a boundary layer in a strong favourable pressure gradient leading to relaminarization. Using a narrow-band filter with suitable mid-band frequencies, the properties of the fine-scale structure (appearing as high frequency pulses in the filtered signal) were analysed using the variable discriminator level technique employed earlier by Rao, Narasimha & Badri Narayanan (1971). It was found that, irrespective of the type of flow, the characteristic pulse frequency (say Np) defined by Rao et al. was about 0·6 times the frequency of the zero crossings. It was also found that, over the small range of Reynolds numbers tested, the ratio of the width of the fine-scale regions to the Kolmogorov scale increased linearly with Reynolds number in grid turbulence as well as in flat-plate boundarylayer flow. Nearly lognormal distributions were exhibited by this ratio as well as by the interval between successive zero crossings. The values of Np and of the zero-crossing rate were found to be nearly constant across the boundary layer, except towards its outer edge and very near the wall. In the zero-pressure-gradient boundary-layer flow, very near the wall the high frequency pulses were found to occur mostly when the longitudinal velocity fluctuation u was positive (i.e. above the mean), whereas in the outer part of the boundary layer the pulses more often occurred when u was negative. During acceleration this correlation between the fine-scale motion and the sign of u was less marked.

Orientational melting and reorientational motion in a cubane molecular crystal: A molecular simulation study

Detailed molecular simulations are carried out to investigate the effect of temperature on orientational order in cubane molecular crystal. We report a transition from an orientationally ordered to an orientationally disordered plastic crystalline phase in the temperature range 425-450 K. This is similar to the experimentally reported transition at 395 K. The nature of this transition is first order and is associated with a 4.8% increase in unit Cell volume that is comparable to the experimentally reported unit cell volume change of 5.4% (Phys. Rev. Lett. 1997, 78, 4938). An orientational order parameter, eta(T), has been defined in terms of average angle of libration of a molecular 3-fold axis and the orientational melting has been characterized by using eta(T). The orientational melting is associated with an anomaly in specific heat at constant pressure (C-p) and compressibility (kappa). The enthalpy of transition and entropy of transition associated with this orientational melting are 20.8 J mol(-1) and 0.046 J mol(-1) K-1, respectively. The structure of crystalline as well as plastic crystalline phases is characterized by using various radial distribution functions and orientational distribution functions. The coefficient of thermal expansion of the plastic crystalline phase is more than twice that of the crystalline phase.

A Geometric Algorithm for Learning Oblique Decision Trees

In this paper we present a novel algorithm for learning oblique decision trees. Most of the current decision tree algorithms rely on impurity measures to assess goodness of hyperplanes at each node. These impurity measures do not properly capture the geometric structures in the data. Motivated by this, our algorithm uses a strategy, based on some recent variants of SVM, to assess the hyperplanes in such a way that the geometric structure in the data is taken into account. We show through empirical studies that our method is effective.

Multiple support excitations of open-plane frames by a filtered white noise and soil-structure interaction

Seismic structural design is essentially the estimation of structural response to a forced motion, which may be deterministic or stochastic, imposed on the ground. The assumption that the same ground motion acts at every point of the base of the structure (or at every support) is not always justifiable; particularly in case of very large structures when considerable spatial variability in ground motion can exist over significant distances example long span bridges. This variability is partly due to the delay in arrival of the excitation at different supports (which is called the wave passage effect) and due to heterogeneity in the ground medium which results in incoherency and local effects. The current study examines the influence of the wave passage effect (in terms of delay in arrival of horizontal ground excitation at different supports and neglecting transmission through the structure) on the response of a few open-plane frame building structures with soil-structure interaction. The ground acceleration has been modeled by a suitably filtered white noise. As a special case, the ground excitation at different supports has also been treated as statistically independent to model the extreme case of incoherence due to local effects and due to modifications to the ground motion resulting from wave reflections and refractions in heterogeneous soil media. The results indicate that, even for relatively short spanned building frames, wave passage effect can be significant. In the absence of soil-structure interaction, it can significantly increase the root mean square (rms) value of the shear in extreme end columns for the stiffer frames but has negligible effect on the flexible frames when total displacements are considered. It is seen that pseudo-static displacements increasingly contribute to the rms value of column shear as the time delay increases both for the stiffer and for the more flexible frames. When soil-structure interaction is considered, wave passage effect (in terms of total displacements) is significant only for low soil shear modulus, G. values (where soil-structure interaction significantly lowers the fundamental frequency) and for stiff frames. The contribution of pseudo-static displacement to these rms values is found to decrease with increase in G. In general, wave passage effect for most interactive frames is insignificant compared to the attenuating effect a decrease in G, has on the response of the interactive structure to uniform support excitation. When the excitations at different supports are statistically independent, it is seen that for both the stiff and flexible frames, the rms value of the column shear in extreme end columns is several times larger (more for the stiffer frames) than the value corresponding to uniform base excitation with the pseudo-static displacements contributing over 99% of the rms value of column shear. Soil-structure interaction has an attenuating effect on the rms value of the column shear, the effect decreasing with increase in G,. Here too, the pseudo-static displacements contribute very largely to the column shear. The influence of the wave passage effect on the response of three 2-bay frames with and without soil-structure interaction to a recorded horizontal accelerogram is also examined. (C) 2010 Elsevier Ltd. All rights reserved.

Structure and dynamics of benzene in one-dimensional channels

Molecular dynamics investigation of benzene in one-dimensional channel systems A1PO(4)-5, VPI-5, and carbon nanotube is reported. The results suggest that, in all the three host systems, the plane of benzene is almost perpendicular to the channel axis when the molecule is near the center of the channel and the plane of benzene is parallel to the channel axis when the molecule is near the wall of the channel. The density distribution of benzene as a function of channel length, z and the radial distance, r, from the channel axis is also different in the three host structures. Anisotropy in translational diffusion coefficient, calculated in body-fixed frame of benzene, suggests that benzene prefers to move with its plane parallel to the direction of motion in A1PO(4)-5 and VPI-5 whereas in carbon nanotube the motion occurs predominantly with the plane of the benzene perpendicular to the direction of motion.;Anisotropy associated with the rotational motion is seen to alter significantly in confinement as compared to liquid benzene. In A1PO(4)-5, the rotational anisotropy is reversed as compared to liquid benzene thereby suggesting that anisotropy arising out of molecular geometry can be reduced. Reorientational correlation times for C-6 and C-2 axes Of benzene are reported. Apart from the inertial decay of reorientational correlation function due to free, rotation, two other distinct regimes of decay are observed in narrower channels (AIPO(4)-5 and carbon nanotube): (i) an initial fast decay (0.5-2 ps) and (ii) a slower decay (>2 ps) of reorientational correlation function where C-6 decays slower than C-2 Similar to what is observed in liquid benzene. In the initial fast decay, it is seen that the decay for C-6 is faster than C-2 which is in contrast to what is observed in liquid benzene or for benzene confined in VPI-5.

An Efficient Random Walk Strategy for Sampling Based Robot Motion Planners

Sampling based planners have been successful in path planning of robots with many degrees of freedom, but still remains ineffective when the configuration space has a narrow passage. We present a new technique based on a random walk strategy to generate samples in narrow regions quickly, thus improving efficiency of Probabilistic Roadmap Planners. The algorithm substantially reduces instances of collision checking and thereby decreases computational time. The method is powerful even for cases where the structure of the narrow passage is not known, thus giving significant improvement over other known methods.

Coupled electro-mechanical response of an electroactive polymer cantilever structure and its application in energy harvesting

Ionic polymer-metal composites (IPMC), piezoelectric polymer composites and nematic elastomer composites are materials, which exhibit characteristics of both sensors and actuators. Large deformation and curvature are observed in these systems when electric potential is applied. Effects of geometric non-linearity due to the chargeinduced motion in these materials are poorly understood. In this paper, a coupled model for understanding the behavior of an ionic polymer beam undergoing large deformation and large curvature is presented. Maxwell's equations and charge transport equations are considered which couple the distribution of the ion concentration and the pressure gradient along length of a cantilever beam with interdigital electrodes. A nonlinear constitutive model is derived accounting for the visco-elasto-plastic behavior of these polymers and based on the hypothesis that the presence of electrical charge stretches/contracts bonds, which give rise to electrical field dependent softening/hardening. Polymer chain orientation in statistical sense plays a role on such softening or hardening. Elementary beam kinematics with large curvature is considered. A model for understanding the deformation due to electrostatic repulsion between asymmetrical charge distributions across the cross-sections is presented. Experimental evidence that Silver(Ag) nanoparticle coated IPMCs can be used for energy harvesting is reported. An IPMC strip is vibrated in different environments and the electric power against a resistive load is measured. The electrical power generated was observed to vary with the environment with maximum power being generated when the strip is in wet state. IPMC based energy harvesting systems have potential applications in tidal wave energy harvesting, residual environmental energy harvesting to power MEMS and NEMS devices.