Dynamic load balancing for the distributed mining of molecular structures

In molecular biology, it is often desirable to find common properties in large numbers of drug candidates. One family of methods stems from the data mining community, where algorithms to find frequent graphs have received increasing attention over the past years. However, the computational complexity of the underlying problem and the large amount of data to be explored essentially render sequential algorithms useless. In this paper, we present a distributed approach to the frequent subgraph mining problem to discover interesting patterns in molecular compounds. This problem is characterized by a highly irregular search tree, whereby no reliable workload prediction is available. We describe the three main aspects of the proposed distributed algorithm, namely, a dynamic partitioning of the search space, a distribution process based on a peer-to-peer communication framework, and a novel receiverinitiated load balancing algorithm. The effectiveness of the distributed method has been evaluated on the well-known National Cancer Institute’s HIV-screening data set, where we were able to show close-to linear speedup in a network of workstations. The proposed approach also allows for dynamic resource aggregation in a non dedicated computational environment. These features make it suitable for large-scale, multi-domain, heterogeneous environments, such as computational grids.

Distributed mining of molecular fragments

In real world applications sequential algorithms of data mining and data exploration are often unsuitable for datasets with enormous size, high-dimensionality and complex data structure. Grid computing promises unprecedented opportunities for unlimited computing and storage resources. In this context there is the necessity to develop high performance distributed data mining algorithms. However, the computational complexity of the problem and the large amount of data to be explored often make the design of large scale applications particularly challenging. In this paper we present the first distributed formulation of a frequent subgraph mining algorithm for discriminative fragments of molecular compounds. Two distributed approaches have been developed and compared on the well known National Cancer Institute’s HIV-screening dataset. We present experimental results on a small-scale computing environment.

Self-assembled biomimetic [2Fe2S]-hydrogenase-based photocatalyst for molecular hydrogen evolution

The large-scale production of clean energy is one of the major challenges society is currently facing. Molecular hydrogen is envisaged as a key green fuel for the future, but it becomes a sustainable alternative for classical fuels only if it is also produced in a clean fashion. Here, we report a supramolecular biomimetic approach to form a catalyst that produces molecular hydrogen using light as the energy source. It is composed of an assembly of chromophores to a bis(thiolate)-bridged diiron ([2Fe2S]) based hydrogenase catalyst. The supramolecular building block approach introduced in this article enabled the easy formation of a series of complexes, which are all thoroughly characterized, revealing that the photoactivity of the catalyst assembly strongly depends on its nature. The active species, formed from different complexes, appears to be the [Fe-2(mu-pdt)(CO)(4){PPh2(4-py)}(2)] (3) with 2 different types of porphyrins (5a and 5b) coordinated to it. The modular supramolecular approach was important in this study as with a limited number of building blocks several different complexes were generated.

Constraint release in entangled binary blends of linear polymers: a molecular dynamics study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.

WAXS studies of global molecular orientation induced in nematic liquid crystals by simple shear flow

Global molecular orientation function coefficients for the nematic liquid crystal 4-cyano 4'-nn -pentylbiphenyl (5CB) in shear flow are presented, being extracted from 2-dimensional Wide-Angle X-ray Scattering data. A linear increase in orientation parameter P2 is observed with a logarithmic increase in shear rate. It is proposed that this arises from an increased number of LC directors aligning to the shear axis. Upon cessation of shear flow, the anisotropy is seen to relax away completely, over a time scale which is inversely proportional to the previously applied shear rate.

Genome scale reconstruction of a Salmonella metabolic model comparison of similarity and differences with a commensal Escherichia coli strain

Salmonella are closely related to commensal Escherichia coli but have gained virulence factors enabling them to behave as enteric pathogens. Less well studied are the similarities and differences that exist between the metabolic properties of these organisms that may contribute toward niche adaptation of Salmonella pathogens. To address this, we have constructed a genome scale Salmonella metabolic model (iMA945). The model comprises 945 open reading frames or genes, 1964 reactions, and 1036 metabolites. There was significant overlap with genes present in E. coli MG1655 model iAF1260. In silico growth predictions were simulated using the model on different carbon, nitrogen, phosphorous, and sulfur sources. These were compared with substrate utilization data gathered from high throughput phenotyping microarrays revealing good agreement. Of the compounds tested, the majority were utilizable by both Salmonella and E. coli. Nevertheless a number of differences were identified both between Salmonella and E. coli and also within the Salmonella strains included. These differences provide valuable insight into differences between a commensal and a closely related pathogen and within different pathogenic strains opening new avenues for future explorations.

Scanning tunneling microscopy and molecular dynamics study of the Li2TiO3(001) surface

We have investigated the (001) surface structure of lithium titanate (Li2TiO3) using auger electron spectroscopy (AES), low-energy electron diffraction (LEED), and scanning tunneling microscopy (STM). Li2TiO3 is a potential fusion reactor blanket material. After annealing at 1200 K, LEED demonstrated that the Li2TiO3(001) surface was well ordered and not reconstructed. STM imaging showed that terraces are separated in height by about 0.3 nm suggesting a single termination layer. Moreover, hexagonal patterns with a periodicity of ∼0.4 nm are observed. On the basis of molecular dynamics (MD) simulations, these are interpreted as a dynamic arrangement of Li atoms.

Syntheses, spectroelectrochemical studies, and molecular and electronic structures of ferrocenyl ene-diynes

The readily available complex 1,1-dibromo-2-ferrocenylethylene provides a convenient entry point for the preparation of a wide range of cross-conjugated 1,1-bis(alkynyl)-2-ferrocenylethenes through simple Pd(0)/Cu(I)-mediated cross-coupling reactions with 1-alkynes. The ferrocene moiety in compounds of the general form FcCHC(CCR)2 is essentially electronically isolated from the cross-conjugated π system, as evidenced by IR and UV−vis spectroelectrochemical experiments and quantum chemical calculations. In contrast to the other examples which give stable ferrocenium derivatives upon electrochemical oxidation, the aniline derivatives [FcCHC(CCC6H4NH2-4)2]+ and [FcCHC(CCC6H4NMe2-4)2]+ proved to be unstable on the time scale of the spectroelectrochemical experiments, leading to passivation of the electrode surface over time. There is no significant thermodynamic stabilization of the radical anion [FcCHC(CCC6H4NO2-4)2]− relative to the neutral and dianionic analogues, although the dianion [FcCHC(CCC6H4NO2- 4)2]2− could be studied as a relatively chemically stable species and is well described in terms of two linked nitrophenyl radicals. The capacity to introduce a relatively isolated point charge at the periphery of the cross-conjugated π system appears to make these complexes useful templates for the construction of electrochemically gated quantum interference transistors.

Global scale transcriptome analysis of Arabidopsis embryogenesis in vitro

Background Somatic embryogenesis (SE) in plants is a process by which embryos are generated directly from somatic cells, rather than from the fused products of male and female gametes. Despite the detailed expression analysis of several somatic-to-embryonic marker genes, a comprehensive understanding of SE at a molecular level is still lacking. The present study was designed to generate high resolution transcriptome datasets for early SE providing the way for future research to understand the underlying molecular mechanisms that regulate this process. We sequenced Arabidopsis thaliana somatic embryos collected from three distinct developmental time-points (5, 10 and 15 d after in vitro culture) using the Illumina HiSeq 2000 platform. Results This study yielded a total of 426,001,826 sequence reads mapped to 26,520 genes in the A. thaliana reference genome. Analysis of embryonic cultures after 5 and 10 d showed differential expression of 1,195 genes; these included 778 genes that were more highly expressed after 5 d as compared to 10 d. Moreover, 1,718 genes were differentially expressed in embryonic cultures between 10 and 15 d. Our data also showed at least eight different expression patterns during early SE; the majority of genes are transcriptionally more active in embryos after 5 d. Comparison of transcriptomes derived from somatic embryos and leaf tissues revealed that at least 4,951 genes are transcriptionally more active in embryos than in the leaf; increased expression of genes involved in DNA cytosine methylation and histone deacetylation were noted in embryogenic tissues. In silico expression analysis based on microarray data found that approximately 5% of these genes are transcriptionally more active in somatic embryos than in actively dividing callus and non-dividing leaf tissues. Moreover, this identified 49 genes expressed at a higher level in somatic embryos than in other tissues. This included several genes with unknown function, as well as others related to oxidative and osmotic stress, and auxin signalling. Conclusions The transcriptome information provided here will form the foundation for future research on genetic and epigenetic control of plant embryogenesis at a molecular level. In follow-up studies, these data could be used to construct a regulatory network for SE; the genes more highly expressed in somatic embryos than in vegetative tissues can be considered as potential candidates to validate these networks.

Rainfastness of poly(vinyl alcohol) deposits on Vicia faba leaf surfaces: from laboratory-scale washing to simulated rain

Rainfastness is the ability of agrochemical deposits to resist wash-off by rain and other related environmental phenomena. This work reports laboratory-scale and raintower studies of the rainfastness of fluorescently labeled poly(vinyl alcohol) (PVA) using fluorescent microscopy combined with image analysis. Samples of hydrolyzed PVA exhibit improved rainfastness over a threshold molecular weight, which correlates with PVA film dissolution, swelling, and crystalline properties. It was also established that the rainfastness of PVA scaled with the molecular weight over this threshold. These PVA samples were further characterized in order to determine the effect of the crystallinity on rainfastness. The quantification of rainfastness is of great interest to the field of agrochemical formulation development in order to improve the efficacy of pesticides and their adjuvants.

Histomorphometric Evaluation of Bioceramic Molecular Impregnated and Dual Acid-Etched Implant Surfaces in the Human Posterior Maxilla

Background: Physical and bioceramic incorporation surface treatments at the nanometer scale showed higher means of bone-to-implant contact (BIC) and torque values compared with surface topography at the micrometer scale; however, the literature concerning the effect of nanometer scale parameters is sparse. Purpose: The aim of this study was to evaluate the influence of two different implant surfaces on the percentage bone-to-implant contact (BIC%) and bone osteocyte density in the human posterior maxilla after 2 months of unloaded healing. Materials and Methods: The implants utilized presented dual acid-etched (DAE) surface and a bioceramic molecular impregnated treatment (Ossean(R), Intra-Lock International, Boca Raton, FL, USA) serving as control and test, respectively. Ten subjects (59 1 9 years of age) received two implants (one of each surface) during conventional implant surgery in the posterior maxilla. After the non-loaded period of 2 months, the implants and the surrounding tissue were removed by means of a trephine and were non-decalcified processed for ground sectioning and analysis of BIC%, bone density in threaded area (BA%), and osteocyte index (Oi). Results: Two DAE implants were found to be clinically unstable at time of retrieval. Histometric evaluation showed significantly higher BIC% and Oi for the test compared to the control surface (p < .05), and that BA% was not significantly different between groups. Wilcoxon matched pairs test was used to compare the differences of histomorphometric variables between implant surfaces. The significance test was conducted at a 5% level of significance. Conclusion: The histological data suggest that the bioceramic molecular impregnated surface-treated implants positively modulated bone healing at early implantation times compared to the DAE surface.

Screening of ARHSP-TCC Patients Expands the Spectrum of SPG11 Mutations and Includes a Large Scale Gene Deletion

Autosomal recessive spastic paraplegia with thinning of corpus callosum (ARHSP-TCC) is a complex form of HSP initially described in Japan but subsequently reported to have a worldwide distribution with a particular high frequency in multiple families from the Mediterranean basin. We recently showed that ARHSP-TCC is commonly associated with mutations in SPG11/KIAA1840 on chromosome 15q. We have now screened a collection of new patients mainly originating from Italy and Brazil, in order to further ascertain the spectrum of mutations in SPG11, enlarge the ethnic origin of SPG11 patients, determine the relative frequency at the level of single Countries (i.e., Italy), and establish whether there is one or more common mutation. In 25 index cases we identified 32 mutations; 22 are novel, including 9 nonsense, 3 small deletions, 4 insertions, 1 in/del, 1 small duplication, 1 missense, 2 splice-site, and for the first time a large genomic rearrangement. This brings the total number of SPG11 mutated patients in the SPATAX collection to 111 cases in 44 families and in 17 isolated cases, from 16 Countries, all assessed using homogeneous clinical criteria. While expanding the spectrum of mutations in SPG11, this larger series also corroborated the notion that even within apparently homogeneous population a molecular diagnosis cannot be achieved without full gene sequencing. (C) 2008 Wiley-Liss, Inc.

Resilience of protein-protein interaction networks as determined by their large-scale topological features

The relationship between the structure and function of biological networks constitutes a fundamental issue in systems biology. Particularly, the structure of protein-protein interaction networks is related to important biological functions. In this work, we investigated how such a resilience is determined by the large scale features of the respective networks. Four species are taken into account, namely yeast Saccharomyces cerevisiae, worm Caenorhabditis elegans, fly Drosophila melanogaster and Homo sapiens. We adopted two entropy-related measurements (degree entropy and dynamic entropy) in order to quantify the overall degree of robustness of these networks. We verified that while they exhibit similar structural variations under random node removal, they differ significantly when subjected to intentional attacks (hub removal). As a matter of fact, more complex species tended to exhibit more robust networks. More specifically, we quantified how six important measurements of the networks topology (namely clustering coefficient, average degree of neighbors, average shortest path length, diameter, assortativity coefficient, and slope of the power law degree distribution) correlated with the two entropy measurements. Our results revealed that the fraction of hubs and the average neighbor degree contribute significantly for the resilience of networks. In addition, the topological analysis of the removed hubs indicated that the presence of alternative paths between the proteins connected to hubs tend to reinforce resilience. The performed analysis helps to understand how resilience is underlain in networks and can be applied to the development of protein network models.

PACKMOL: A Package for Building Initial Configurations for Molecular Dynamics Simulations

Adequate initial configurations for molecular dynamics simulations consist of arrangements of molecules distributed in space in such a way to approximately represent the system`s overall structure. In order that the simulations are not disrupted by large van der Waals repulsive interactions, atoms from different molecules Must keep safe pairwise distances. Obtaining Such a molecular arrangement can be considered it packing problem: Each type molecule must satisfy spatial constraints related to the geometry of the system, and the distance between atoms of different molecules Must be greater than some specified tolerance. We have developed a code able to pack millions of atoms. grouped in arbitrarily complex molecules, inside a variety of three-dimensional regions. The regions may be intersections of spheres, ellipses, cylinders, planes, or boxes. The user must provide only the structure of one molecule of each type and the geometrical constraints that each type of molecule must satisfy. Building complex mixtures, interfaces, solvating biomolecules in water, other solvents, or mixtures of solvents, is straight forward. In addition. different atoms belonging to the same molecule may also be restricted to different spatial regions, in Such a way that more ordered molecular arrangements call be built, as micelles. lipid double-layers, etc. The packing time for state-of-the-art molecular dynamics systems varies front a few seconds to a few Minutes in a personal Computer. The input files are simple and Currently compatible with PDB, Tinker, Molden, or Moldy coordinate files. The package is distributed as free software and call be downloaded front http://www.ime.unicamp.br/similar to martinez/packmol/. (C) 2009 Wiley Periodicals. Inc. J Comput Chem 30: 2157-2164, 2009

Cellular automaton model for molecular traffic jams

We consider the time evolution of an exactly solvable cellular automaton with random initial conditions both in the large-scale hydrodynamic limit and on the microscopic level. This model is a version of the totally asymmetric simple exclusion process with sublattice parallel update and thus may serve as a model for studying traffic jams in systems of self-driven particles. We study the emergence of shocks from the microscopic dynamics of the model. In particular, we introduce shock measures whose time evolution we can compute explicitly, both in the thermodynamic limit and for open boundaries where a boundary-induced phase transition driven by the motion of a shock occurs. The motion of the shock, which results from the collective dynamics of the exclusion particles, is a random walk with an internal degree of freedom that determines the jump direction. This type of hopping dynamics is reminiscent of some transport phenomena in biological systems.