How epigenetic evolution can guide genetic evolution (abstract)

The expression level of a gene in future generations can be modified both by genetic mutations and by the attachment of methyl groups to the DNA. Since the DNA methylation pattern along a genome is inherited, methylation patterns constitute a significant epigenetic inheritance mechanism that is subject to evolution by natural selection. The variation rate of methylation patterns is generally higher than that of DNA which suggests that evolution of methylation patterns might be more rapid than that of genetic evolution. But, common consequences of methylation, such as reduced expression of methylated genes, could also be produced by genetic changes and these would have higher heritability. The question we address in this work is how the evolution of epigenetic methylation-dependent phenotypes might interact with the evolution of genetic DNA-determined phenotypes. There is no biological mechanism known to directly transfer methyl groups into equivalent DNA changes. However, in principle an indirect mechanism could cause evolved methylation patterns to enable the subsequent evolution of equivalent genetic patterns in a manner analogous to the Baldwin effect (Baldwin, Am. Nat., 30:441-451, 1896; Jablonka et al, TREE, 13:206-210, 1998). The Baldwin effect describes how non-heritable acquired characteristics can influence the evolution of equivalent genetic characteristics without any direct Lamarckian inheritance of acquired characters. This occurs because the ability to acquire or learn a new behaviour changes the selective pressures acting on genetic changes. Specifically, genetic changes that support this behaviour, e.g. by reducing learning time by making a small part of the behaviour genetically innate, may be selected for when the learning mechanism is present even though these same genetic changes may not be selected for when the learning mechanism is absent. Over generations, the modified selection pressures so produced can cause genetic assimilation of a phenotype that was previously acquired, even to the extent of making the acquisition mechanism subsequently redundant. Thus a learned behaviour can guide the evolution of an equivalent innate behaviour (Hinton & Nowlan, Complex Systems, 1: 495-502, 1987). In the Baldwin effect a rapid mechanism of lifetime adaptation guides the relatively slow genetic evolution of the same behaviour. By analogy, Jablonka et al have suggested that “genetic adaptations may be guided by heritable induced or learnt phenotypic adaptations”. Here we hypothesise that “inherited epigenetic variations may be able to ‘hold’ an adapted state for long enough to allow similar genetic variations to catch up”, as they put it, even if the epigenetic variations are not induced or learnt but simply evolved by natural selection on methylation patterns. We assume that an individual may only express one phenotype in its lifetime, but that a given genome will persist relatively unchanged on a timescale that allows its methylome to adapt by natural selection. Thus, in contrast to the Baldwin effect, in this case two mechanisms of evolution by natural selection are coupled — one acting at a different variation rate from the other. We present a simple model to illustrate how a rapidly evolving methylome can guide a slowly evolving but highly-heritable genome. This is used to show that methylome evolution can enable genetic evolution to cross fitness valleys that would otherwise require multiple genetic changes that were each selected against. This finding suggests that the relatively rapid evolution of methylation patterns can produce novel phenotypes that are subsequently genetically assimilated in DNA evolution without direct transfer or appeal to induced phenotypes. This can enable the genetic evolution of new phenotypes that would not be found by genetic evolution alone, even if methylation is not significant in the ultimate phenotype

Report on the sixth blind test of organic crystal structure prediction methods.

The sixth blind test of organic crystal structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal and a bulky flexible molecule. This blind test has seen substantial growth in the number of participants, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and `best practices' for performing CSP calculations. All of the targets, apart from a single potentially disordered Z' = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms.The organisers and participants are very grateful to the crystallographers who supplied the candidate structures: Dr. Peter Horton (XXII), Dr. Brian Samas (XXIII), Prof. Bruce Foxman (XXIV), and Prof. Kraig Wheeler (XXV and XXVI). We are also grateful to Dr. Emma Sharp and colleagues at Johnson Matthey (Pharmorphix) for the polymorph screening of XXVI, as well as numerous colleagues at the CCDC for assistance in organising the blind test. Submission 2: We acknowledge Dr. Oliver Korb for numerous useful discussions. Submission 3: The Day group acknowledge the use of the IRIDIS High Performance Computing Facility, and associated support services at the University of Southampton, in the completion of this work. We acknowledge funding from the EPSRC (grants EP/J01110X/1 and EP/K018132/1) and the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC through grant agreements n. 307358 (ERC-stG- 2012-ANGLE) and n. 321156 (ERC-AG-PE5-ROBOT). Submission 4: I am grateful to Mikhail Kuzminskii for calculations of molecular structures on Gaussian 98 program in the Institute of Organic Chemistry RAS. The Russian Foundation for Basic Research is acknowledged for financial support (14-03-01091). Submission 5: Toine Schreurs provided computer facilities and assistance. I am grateful to Matthew Habgood at AWE company for providing a travel grant. Submission 6: We would like to acknowledge support of this work by GlaxoSmithKline, Merck, and Vertex. Submission 7: The research was financially supported by the VIDI Research Program 700.10.427, which is financed by The Netherlands Organisation for Scientific Research (NWO), and the European Research Council (ERC-2010-StG, grant agreement n. 259510-KISMOL). We acknowledge the support of the Foundation for Fundamental Research on Matter (FOM). Supercomputer facilities were provided by the National Computing Facilities Foundation (NCF). Submission 8: Computer resources were provided by the Center for High Performance Computing at the University of Utah and the Extreme Science and Engineering Discovery Environment (XSEDE), supported by NSF grant number ACI-1053575. MBF and GIP acknowledge the support from the University of Buenos Aires and the Argentinian Research Council. Submission 9: We thank Dr. Bouke van Eijck for his valuable advice on our predicted structure of XXV. We thank the promotion office for TUT programs on advanced simulation engineering (ADSIM), the leading program for training brain information architects (BRAIN), and the information and media center (IMC) at Toyohashi University of Technology for the use of the TUT supercomputer systems and application software. We also thank the ACCMS at Kyoto University for the use of their supercomputer. In addition, we wish to thank financial supports from Conflex Corp. and Ministry of Education, Culture, Sports, Science and Technology. Submission 12: We thank Leslie Leiserowitz from the Weizmann Institute of Science and Geoffrey Hutchinson from the University of Pittsburgh for helpful discussions. We thank Adam Scovel at the Argonne Leadership Computing Facility (ALCF) for technical support. Work at Tulane University was funded by the Louisiana Board of Regents Award # LEQSF(2014-17)-RD-A-10 “Toward Crystal Engineering from First Principles”, by the NSF award # EPS-1003897 “The Louisiana Alliance for Simulation-Guided Materials Applications (LA-SiGMA)”, and by the Tulane Committee on Research Summer Fellowship. Work at the Technical University of Munich was supported by the Solar Technologies Go Hybrid initiative of the State of Bavaria, Germany. Computer time was provided by the Argonne Leadership Computing Facility (ALCF), which is supported by the Office of Science of the U.S. Department of Energy under contract DE-AC02-06CH11357. Submission 13: This work would not have been possible without funding from Khalifa University’s College of Engineering. I would like to acknowledge Prof. Robert Bennell and Prof. Bayan Sharif for supporting me in acquiring the resources needed to carry out this research. Dr. Louise Price is thanked for her guidance on the use of DMACRYS and NEIGHCRYS during the course of this research. She is also thanked for useful discussions and numerous e-mail exchanges concerning the blind test. Prof. Sarah Price is acknowledged for her support and guidance over many years and for providing access to DMACRYS and NEIGHCRYS. Submission 15: The work was supported by the United Kingdom’s Engineering and Physical Sciences Research Council (EPSRC) (EP/J003840/1, EP/J014958/1) and was made possible through access to computational resources and support from the High Performance Computing Cluster at Imperial College London. We are grateful to Professor Sarah L. Price for supplying the DMACRYS code for use within CrystalOptimizer, and to her and her research group for support with DMACRYS and feedback on CrystalPredictor and CrystalOptimizer. Submission 16: R. J. N. acknowledges financial support from the Engineering and Physical Sciences Research Council (EPSRC) of the U.K. [EP/J017639/1]. R. J. N. and C. J. P. acknowledge use of the Archer facilities of the U.K.’s national high-performance computing service (for which access was obtained via the UKCP consortium [EP/K014560/1]). C. J. P. also acknowledges a Leadership Fellowship Grant [EP/K013688/1]. B. M. acknowledges Robinson College, Cambridge, and the Cambridge Philosophical Society for a Henslow Research Fellowship. Submission 17: The work at the University of Delaware was supported by the Army Research Office under Grant W911NF-13-1- 0387 and by the National Science Foundation Grant CHE-1152899. The work at the University of Silesia was supported by the Polish National Science Centre Grant No. DEC-2012/05/B/ST4/00086. Submission 18: We would like to thank Constantinos Pantelides, Claire Adjiman and Isaac Sugden of Imperial College for their support of our use of CrystalPredictor and CrystalOptimizer in this and Submission 19. The CSP work of the group is supported by EPSRC, though grant ESPRC EP/K039229/1, and Eli Lilly. The PhD students support: RKH by a joint UCL Max-Planck Society Magdeburg Impact studentship, REW by a UCL Impact studentship; LI by the Cambridge Crystallographic Data Centre and the M3S Centre for Doctoral Training (EPSRC EP/G036675/1). Submission 19: The potential generation work at the University of Delaware was supported by the Army Research Office under Grant W911NF-13-1-0387 and by the National Science Foundation Grant CHE-1152899. Submission 20: The work at New York University was supported, in part, by the U.S. Army Research Laboratory and the U.S. Army Research Office under contract/grant number W911NF-13-1-0387 (MET and LV) and, in part, by the Materials Research Science and Engineering Center (MRSEC) program of the National Science Foundation under Award Number DMR-1420073 (MET and ES). The work at the University of Delaware was supported by the U.S. Army Research Laboratory and the U.S. Army Research Office under contract/grant number W911NF-13-1- 0387 and by the National Science Foundation Grant CHE-1152899. Submission 21: We thank the National Science Foundation (DMR-1231586), the Government of Russian Federation (Grant No. 14.A12.31.0003), the Foreign Talents Introduction and Academic Exchange Program (No. B08040) and the Russian Science Foundation, project no. 14-43-00052, base organization Photochemistry Center of the Russian Academy of Sciences. Calculations were performed on the Rurik supercomputer at Moscow Institute of Physics and Technology. Submission 22: The computational results presented have been achieved in part using the Vienna Scientific Cluster (VSC). Submission 24: The potential generation work at the University of Delaware was supported by the Army Research Office under Grant W911NF-13-1-0387 and by the National Science Foundation Grant CHE-1152899. Submission 25: J.H. and A.T. acknowledge the support from the Deutsche Forschungsgemeinschaft under the program DFG-SPP 1807. H-Y.K., R.A.D., and R.C. acknowledge support from the Department of Energy (DOE) under Grant Nos. DE-SC0008626. This research used resources of the Argonne Leadership Computing Facility at Argonne National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-06CH11357. This research used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DEAC02-05CH11231. Additional computational resources were provided by the Terascale Infrastructure for Groundbreaking Research in Science and Engineering (TIGRESS) High Performance Computing Center and Visualization Laboratory at Princeton University.This is the final version of the article. It first appeared from Wiley via http://dx.doi.org/10.1107/S2052520616007447

Report on the sixth blind test of organic crystal-structure prediction methods

The sixth blind test of organic crystal-structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal, and a bulky flexible molecule. This blind test has seen substantial growth in the number of submissions, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and "best practices" for performing CSP calculations. All of the targets, apart from a single potentially disordered Z` = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms

Interaction With the Lipid Membrane Influences Fentanyl Pharmacology

Overdose deaths from fentanyl have reached epidemic proportions in the USA and are increasing worldwide. Fentanyl is a potent opioid agonist that is less well reversed by naloxone than morphine. Due to fentanyl’s high lipophilicity and elongated structure we hypothesised that its unusual pharmacology may be explained by its interactions with the lipid membrane on route to binding to the μ-opioid receptor (MOPr). Through coarse-grained molecular dynamics simulations, electrophysiological recordings and cell signalling assays, we determined how fentanyl and morphine access the orthosteric pocket of MOPr. Morphine accesses MOPr via the aqueous pathway; first binding to an extracellular vestibule, then diffusing into the orthosteric pocket. In contrast, fentanyl may take a novel route; first partitioning into the membrane, before accessing the orthosteric site by diffusing through a ligand-induced gap between the transmembrane helices. In electrophysiological recordings fentanyl-induced currents returned after washout, suggesting fentanyl deposits in the lipid membrane. However, mutation of residues forming the potential MOPr transmembrane access site did not alter fentanyl’s pharmacological profile in vitro. A high local concentration of fentanyl in the lipid membrane, possibly in combination with a novel lipophilic binding route, may explain the high potency and lower susceptibility of fentanyl to reversal by naloxone

Murder Ballads

Violence and murder have a strong cultural currency, the implications of which should be pursued by those with an interest in law and society, crime, and justice. Murder ballads are songs about death and killing with a history stretching back to the nineteenth century. Drawing out the major themes of this genre can help scholars gain a handle on how murder has been treated in popular culture, thereupon providing an enhanced understanding of the human condition. As an example of such examination, 2016 marked the twentieth anniversary of Nick Cave and the Bad Seeds’ Murder Ballads, their most famous and, perhaps, defining album. More than any other Bad Seeds album, Murder Ballads captures the essence of a band at its most comfortable in exploring the dark and the taboo: violence, killing, death. In producing a whole album on murder, the band left a calling card by which the wider public could define them. This article will explore the album by considering its key themes and, in so doing, reflect on the need to understand the use of murder in such popular music. The use of murder and death in popular music has not been properly studied, yet it offers potential social insight for several fields of study such as law, criminology, and psychology. In particular, little considered issues around the treatment of murder in popular culture such as humour are identified, while others that require greater attention such as attitudes to women are also given due consideratio

Higher airborne pollen concentrations correlated with increased SARS-CoV-2 infection rates, as evidenced from 31 countries across the globe

Pollen exposure weakens the immunity against certain seasonal respiratory viruses by diminishing the antiviral interferon response. Here we investigate whether the same applies to the pandemic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which is sensitive to antiviral interferons, if infection waves coincide with high airborne pollen concentrations. Our original hypothesis was that more airborne pollen would lead to increases in infection rates. To examine this, we performed a cross-sectional and longitudinal data analysis on SARS-CoV-2 infection, airborne pollen, and meteorological factors. Our dataset is the most comprehensive, largest possible worldwide from 130 stations, across 31 countries and five continents. To explicitly investigate the effects of social contact, we additionally considered population density of each study area, as well as lockdown effects, in all possible combinations: without any lockdown, with mixed lockdown−no lockdown regime, and under complete lockdown. We found that airborne pollen, sometimes in synergy with humidity and temperature, explained, on average, 44% of the infection rate variability. Infection rates increased after higher pollen concentrations most frequently during the four previous days. Without lockdown, an increase of pollen abundance by 100 pollen/m3 resulted in a 4% average increase of infection rates. Lockdown halved infection rates under similar pollen concentrations. As there can be no preventive measures against airborne pollen exposure, we suggest wide dissemination of pollen−virus coexposure dire effect information to encourage high-risk individuals to wear particle filter masks during high springtime pollen concentrations.</p

ILC Reference Design Report Volume 1 - Executive Summary

The International Linear Collider (ILC) is a 200-500 GeV center-of-mass high-luminosity linear electron-positron collider, based on 1.3 GHz superconducting radio-frequency (SCRF) accelerating cavities. The ILC has a total footprint of about 31 km and is designed for a peak luminosity of 2x10^34 cm^-2s^-1. This report is the Executive Summary (Volume I) of the four volume Reference Design Report. It gives an overview of the physics at the ILC, the accelerator design and value estimate, the detector concepts, and the next steps towards project realization.The International Linear Collider (ILC) is a 200-500 GeV center-of-mass high-luminosity linear electron-positron collider, based on 1.3 GHz superconducting radio-frequency (SCRF) accelerating cavities. The ILC has a total footprint of about 31 km and is designed for a peak luminosity of 2x10^34 cm^-2s^-1. This report is the Executive Summary (Volume I) of the four volume Reference Design Report. It gives an overview of the physics at the ILC, the accelerator design and value estimate, the detector concepts, and the next steps towards project realization

ILC Reference Design Report Volume 4 - Detectors

This report, Volume IV of the International Linear Collider Reference Design Report, describes the detectors which will record and measure the charged and neutral particles produced in the ILC's high energy e+e- collisions. The physics of the ILC, and the environment of the machine-detector interface, pose new challenges for detector design. Several conceptual designs for the detector promise the needed performance, and ongoing detector R&D is addressing the outstanding technological issues. Two such detectors, operating in push-pull mode, perfectly instrument the ILC interaction region, and access the full potential of ILC physics.This report, Volume IV of the International Linear Collider Reference Design Report, describes the detectors which will record and measure the charged and neutral particles produced in the ILC's high energy e+e- collisions. The physics of the ILC, and the environment of the machine-detector interface, pose new challenges for detector design. Several conceptual designs for the detector promise the needed performance, and ongoing detector R&D is addressing the outstanding technological issues. Two such detectors, operating in push-pull mode, perfectly instrument the ILC interaction region, and access the full potential of ILC physics