The eMinerals project: developing the concept of the virtual organisation to support collaborative work on molecular-scale environmental simulations

The eMinerals project has established an integrated compute and data minigrid infrastructure together with a set of collaborative tools,. The infrastructure is designed to support molecular simulation scientists working together as a virtual organisation aiming to understand a number of strategic processes in environmental science. The eMinerals virtual organisation is now working towards applying this infrastructure to tackle a new generation of scientific problems. This paper describes the achievements of the eMinerals virtual organisation to date, and describes ongoing applications of the virtual organisation infrastructure.

New insights on growth mechanisms of protein clusters at surfaces: an AFM and simulation study

Despite its relevance to a wide range of technological and fundamental areas, a quantitative understanding of protein surface clustering dynamics is often lacking. In inorganic crystal growth, surface clustering of adatoms is well described by diffusion-aggregation models. In such models, the statistical properties of the aggregate arrays often reveal the molecular scale aggregation processes. We investigate the potential of these theories to reveal hitherto hidden facets of protein clustering by carrying out concomitant observations of lysozyme adsorption onto mica surfaces, using atomic force microscopy. and Monte Carlo simulations of cluster nucleation and growth. We find that lysozyme clusters diffuse across the substrate at a rate that varies inversely with size. This result suggests which molecular scale mechanisms are responsible for the mobility of the proteins on the substrate. In addition the surface diffusion coefficient of the monomer can also be extracted from the comparison between experiments and simulations. While concentrating on a model system of lysozyme-on-mica, this 'proof of concept' study successfully demonstrates the potential of our approach to understand and influence more biomedically applicable protein-substrate couples.

Self-organisation in the assembly of gels from mixtures of different dendritic peptide building blocks

This paper investigates dendritic peptides capable of assembling into nanostructured gels, and explores the effect on self-assembly of mixing different molecular building blocks. Thermal measurements, small angle Xray scattering (SAXS) and circular dichroism (CD) spectroscopy are used to probe these materials on macroscopic, nanoscopic and molecular length scales. The results from these investigations demonstrate that in this case, systems with different "size" and "chirality" factors can self-organise, whilst systems with different "shape" factors cannot. The "size" and "chirality" factors are directly connected with the molecular information programmed into the dendritic peptides, whilst the shape factor depends on the group linking these peptides together-this is consistent with molecular recognition hydrogen bond pathways between the peptidic building blocks controlling the ability of these systems to self-recognise. These results demonstrate that mixtures of relatively complex peptides, with only subtle differences on the molecular scale, can self-organise into nanoscale structures, an important step in the spontaneous assembly of ordered systems from complex mixtures.

Self-assembly of two-component gels: stoichiometric control and component selection

Two-component systems capable of self-assembling into soft gel-phase materials are of considerable interest due to their tunability and versatility. This paper investigates two-component gels based on a combination of a L-lysine-based dendron and a rigid diamine spacer (1,4-diaminobenzene or 1,4-diaminocyclohexane). The networked gelator was investigated using thermal measurements, circular dichroism, NMR spectroscopy and small angle neutron scattering (SANS) giving insight into the macroscopic properties, nanostructure and molecular-scale organisation. Surprisingly, all of these techniques confirmed that irrespective of the molar ratio of the components employed, the "solid-like" gel network always consisted of a 1:1 mixture of dendron/diamine. Additionally, the gel network was able to tolerate a significant excess of diamine in the "liquid-like" phase before being disrupted. In the light of this observation, we investigated the ability of the gel network structure to evolve from mixtures of different aromatic diamines present in excess. We found that these two-component gels assembled in a component-selective manner, with the dendron preferentially recognising 1,4-diaminobenzene (>70%). when similar competitor diamines (1,2- and 1,3-diaminobenzene) are present. Furthermore, NMR relaxation measurements demonstrated that the gel based oil 1,4-diaminobenzene was better able to form a selective ternary complex with pyrene than the gel based oil 1,4-diaminocyclohexane, indicative of controlled and selective pi-pi interactions within a three-component assembly. As such, the results ill this paper demonstrate how component selection processes in two-component gel systems call control hierarchical self-assembly.

Terahertz spectroscopy and imaging at the nanoscale for biological and security applications

The chemical specificity of terahertz spectroscopy, when combined with techniques for sub-wavelength sensing, is giving new understanding of processes occurring at the nanometre scale in biological systems and offers the potential for single molecule detection of chemical and biological agents and explosives. In addition, terahertz techniques are enabling the exploration of the fundamental behaviour of light when it interacts with nanoscale optical structures, and are being used to measure ultrafast carrier dynamics, transport and localisation in nanostructures. This chapter will explain how terahertz scale modelling can be used to explore the fundamental physics of nano-optics, it will discuss the terahertz spectroscopy of nanomaterials, terahertz near-field microscopy and other sub-wavelength techniques, and summarise recent developments in the terahertz spectroscopy and imaging of biological systems at the nanoscale. The potential of using these techniques for security applications will be considered.

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.

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.