Searching for stable fullerenes in space with computational chemistry

We report a computational study of the stability and infrared (IR) vibrational spectra of neutral and singly ionised fullerene cages containing between 44 and 70 carbon atoms. The stability is characterised in terms of the standard enthalpy of formation per CC bond, the HOMO-LUMO gap, and the energy required to eliminate a C$_2$ fragment. We compare the simulated IR spectra of these fullerene species to the observed emission spectra of several planetary nebulae (Tc 1, SMP SMC 16, and SMP LMC 56) where strong C$_{60}$ emission has been detected. Although we could not conclusively identify fullerenes other than C$_{60}$ and C$_{70}$, our results point to the possible presence of smaller (44, 50, and 56-atom) cages in those astronomical objects. Observational confirmation of our prediction should become possible when the James Webb Space Telescope comes online.Comment: 11 pages, 13 figures, 1 table. Accepted for publication on MNRA

Metric Features of a Dipolar Model

The lattice spin model, with nearest neighbor ferromagnetic exchange and long range dipolar interaction, is studied by the method of time series for observables based on cluster configurations and associated partitions, such as Shannon entropy, Hamming and Rohlin distances. Previous results based on the two peaks shape of the specific heat, suggested the existence of two possible transitions. By the analysis of the Shannon entropy we are able to prove that the first one is a true phase transition corresponding to a particular melting process of oriented domains, where colored noise is present almost independently of true fractality. The second one is not a real transition and it may be ascribed to a smooth balancing between two geometrical effects: a progressive fragmentation of the big clusters (possibly creating fractals), and the slow onset of a small clusters chaotic phase. Comparison with the nearest neighbor Ising ferromagnetic system points out a substantial difference in the cluster geometrical properties of the two models and in their critical behavior.Comment: 20 pages, 15 figures, submitted to JPhys

Investigating genome wide DNA methylation in bronchial and lung fibroblasts from healthy individuals and individuals with COPD

Rationale: Lung fibroblasts are implicated in respiratory disease pathology including chronic obstructive pulmonary disease (COPD). Phenotypic differences between fibroblasts isolated from the bronchi versus the lung parenchyma have been described but no studies have compared the cell types on a genome wide scale. DNA methylation is a reversible modification of the DNA structure with the ability to affect cell function via the alteration of gene expression. Here we compared genome wide DNA methylation profiles from bronchial and lung fibroblasts and assessed modification to these profiles in cells isolated from individuals with COPD. Methods: DNA was isolated from lung (LgF) and bronchial fibroblasts (BrF) at passage 4 and bisulphite treated. Site specific, quantitative genome wide methylation was determined using the Illumina 450K Infinium Methylation BeadChip array. Linear modelling and DMRcate functions identified differentially methylated sites and regions respectively between BrF and LgF and from cells isolated from healthy individuals versus those with COPD. Results: 3980 CpG (methylation) sites significantly differed, following Bonferroni correction, between BrF and LgF isolated from healthy individuals. These sites had a broad distribution of effect size, with 240 CpG sites displaying a difference in methylation of >50%. 78 of these sites were validated in a second cohort of matched BrF and LgF isolated from the same individuals. There was genomic proximity to these sites and DMRcate was used to refine the individual CpG sites to 5 regions of interest associated with 5 genes; HLX, TWIST1, CREB5, SKAP2 and PRDM16. Differences in methylation were less pronounced when comparing cells isolated from healthy individuals to those with COPD. In BrF 47 DMRcate regions were identified with a maximum difference in methylation of at least 20%. In LgF 3 DMRcate regions were identified with a maximum difference in methylation of at least 20%. Conclusions: DNA methylation profiles are significantly different between BrF and LgF but only small modifications are associated with COPD. Future work will focus on validating a methylation based marker of lung versus bronchial fibroblasts to differentiate cell types by validating our differential DNA methylation observations with gene/protein expression

Dipolar interaction between two-dimensional magnetic particles

We determine the effective dipolar interaction between single domain two-dimensional ferromagnetic particles (islands or dots), taking into account their finite size. The first correction term decays as 1/D^5, where D is the distance between particles. If the particles are arranged in a regular two-dimensional array and are magnetized in plane, we show that the correction term reinforces the antiferromagnetic character of the ground state in a square lattice, and the ferromagnetic one in a triangular lattice. We also determine the dipolar spin-wave spectrum and evaluate how the Curie temperature of an ensemble of magnetic particles scales with the parameters defining the particle array: height and size of each particle, and interparticle distance. Our results show that dipolar coupling between particles might induce ferromagnetic long range order at experimentally relevant temperatures. However, depending on the size of the particles, such a collective phenomenon may be disguised by superparamagnetism.Comment: 11 pages, 5 figure

A systematic review of strategies to recruit and retain primary care doctors

Background There is a workforce crisis in primary care. Previous research has looked at the reasons underlying recruitment and retention problems, but little research has looked at what works to improve recruitment and retention. The aim of this systematic review is to evaluate interventions and strategies used to recruit and retain primary care doctors internationally. Methods A systematic review was undertaken. MEDLINE, EMBASE, CENTRAL and grey literature were searched from inception to January 2015.Articles assessing interventions aimed at recruiting or retaining doctors in high income countries, applicable to primary care doctors were included. No restrictions on language or year of publication. The first author screened all titles and abstracts and a second author screened 20%. Data extraction was carried out by one author and checked by a second. Meta-analysis was not possible due to heterogeneity. Results 51 studies assessing 42 interventions were retrieved. Interventions were categorised into thirteen groups: financial incentives (n=11), recruiting rural students (n=6), international recruitment (n=4), rural or primary care focused undergraduate placements (n=3), rural or underserved postgraduate training (n=3), well-being or peer support initiatives (n=3), marketing (n=2), mixed interventions (n=5), support for professional development or research (n=5), retainer schemes (n=4), re-entry schemes (n=1), specialised recruiters or case managers (n=2) and delayed partnerships (n=2). Studies were of low methodological quality with no RCTs and only 15 studies with a comparison group. Weak evidence supported the use of postgraduate placements in underserved areas, undergraduate rural placements and recruiting students to medical school from rural areas. There was mixed evidence about financial incentives. A marketing campaign was associated with lower recruitment. Conclusions This is the first systematic review of interventions to improve recruitment and retention of primary care doctors. Although the evidence base for recruiting and care doctors is weak and more high quality research is needed, this review found evidence to support undergraduate and postgraduate placements in underserved areas, and selective recruitment of medical students. Other initiatives covered may have potential to improve recruitment and retention of primary care practitioners, but their effectiveness has not been established

New Dynamic Monte Carlo Renormalization Group Method

The dynamical critical exponent of the two-dimensional spin-flip Ising model is evaluated by a Monte Carlo renormalization group method involving a transformation in time. The results agree very well with a finite-size scaling analysis performed on the same data. The value of $z = 2.13 \pm 0.01$ is obtained, which is consistent with most recent estimates

Analytical and computational study of magnetization switching in kinetic Ising systems with demagnetizing fields

An important aspect of real ferromagnetic particles is the demagnetizing field resulting from magnetostatic dipole-dipole interaction, which causes large particles to break up into domains. Sufficiently small particles, however, remain single-domain in equilibrium. This makes such small particles of particular interest as materials for high-density magnetic recording media. In this paper we use analytic arguments and Monte Carlo simulations to study the effect of the demagnetizing field on the dynamics of magnetization switching in two-dimensional, single-domain, kinetic Ising systems. For systems in the ``Stochastic Region,'' where magnetization switching is on average effected by the nucleation and growth of fewer than two well-defined critical droplets, the simulation results can be explained by the dynamics of a simple model in which the free energy is a function only of magnetization. In the ``Multi-Droplet Region,'' a generalization of Avrami's Law involving a magnetization-dependent effective magnetic field gives good agreement with our simulations.Comment: 29 pages, REVTeX 3.0, 10 figures, 2 more figures by request. Submitted Phys. Rev.

Wide-Scale Analysis of Human Functional Transcription Factor Binding Reveals a Strong Bias towards the Transcription Start Site

We introduce a novel method to screen the promoters of a set of genes with shared biological function, against a precompiled library of motifs, and find those motifs which are statistically over-represented in the gene set. The gene sets were obtained from the functional Gene Ontology (GO) classification; for each set and motif we optimized the sequence similarity score threshold, independently for every location window (measured with respect to the TSS), taking into account the location dependent nucleotide heterogeneity along the promoters of the target genes. We performed a high throughput analysis, searching the promoters (from 200bp downstream to 1000bp upstream the TSS), of more than 8000 human and 23,000 mouse genes, for 134 functional Gene Ontology classes and for 412 known DNA motifs. When combined with binding site and location conservation between human and mouse, the method identifies with high probability functional binding sites that regulate groups of biologically related genes. We found many location-sensitive functional binding events and showed that they clustered close to the TSS. Our method and findings were put to several experimental tests. By allowing a "flexible" threshold and combining our functional class and location specific search method with conservation between human and mouse, we are able to identify reliably functional TF binding sites. This is an essential step towards constructing regulatory networks and elucidating the design principles that govern transcriptional regulation of expression. The promoter region proximal to the TSS appears to be of central importance for regulation of transcription in human and mouse, just as it is in bacteria and yeast.Comment: 31 pages, including Supplementary Information and figure

Cancer somatic mutations cluster in a subset of regulatory sites predicted from the ENCODE data

Background: Transcriptional regulation of gene expression is essential for cellular differentiation and function, and defects in the process are associated with cancer. The ENCODE project has mapped potential regulatory sites across the complete genome in many cell types, and these regions have been shown to harbour many of the somatic mutations that occur in cancer cells, suggesting that their effects may drive cancer initiation and development. The ENCODE data suggests a very large number of regulatory sites, and methods are needed to identify those that are most relevant and to connect them to the genes that they control. Methods: Predictive models of gene expression were developed by integrating the ENCODE data for regulation, including transcription factor binding and DNase1 hypersensitivity, with RNA-seq data for gene expression. A penalized regression method was used to identify the most predictive potential regulatory sites for each transcript. Known cancer somatic mutations from the COSMIC database were mapped to potential regulatory sites, and we examined differences in the mapping frequencies associated with sites chosen in regulatory models and other (rejected) sites. The effects of potential confounders, for example replication timing, were considered. Results: Cancer somatic mutations preferentially occupy those regulatory regions chosen in our models as most predictive of gene expression. Conclusion: Our methods have identified a significantly reduced set of regulatory sites that are enriched in cancer somatic mutations and are more predictive of gene expression. This has significance for the mechanistic interpretation of cancer mutations, and the understanding of genetic regulation

Relationship between propagule pressure and colonization pressure in invasion ecology: a test with ships' ballast

Increasing empirical evidence indicates the number of released individuals (i.e. propagule pressure) and number of released species (i.e. colonization pressure) are key determinants of the number of species that successfully invade new habitats. In view of these relationships, and the possibility that ships transport whole communities of organisms, we collected 333 ballast water and sediment samples to investigate the relationship between propagule and colonization pressure for a variety of diverse taxonomic groups (diatoms, dinoflagellates and invertebrates). We also reviewed the scientific literature to compare the number of species transported by ships to those reported in nature. Here, we show that even though ships transport nearly entire local communities, a strong relationship between propagule and colonization pressure exists only for dinoflagellates. Our study provides evidence that colonization pressure of invertebrates and diatoms may fluctuate widely irrespective of propagule pressure. We suggest that the lack of correspondence is explained by reduced uptake of invertebrates into the transport vector and the sensitivity of invertebrates and diatoms to selective pressures during transportation. Selection during transportation is initially evident through decreases in propagule pressure, followed by decreased colonization pressure in the most sensitive taxa