Reactamole: Functional Reactive Molecular Programming

Chemical reaction networks (CRNs) are an important tool for molecular programming, a field that is rapidly expanding our ability to deploy computer programs into biological systems for a variety of applications. However, CRNs are also difficult to work with due to their massively parallel nature, leading to the need for higher-level languages that allow for easier computation with CRNs. Recently, research has been conducted into a variety of higher-level languages for deterministic CRNs but modeling CRN parallelism, managing error accumulation, and finding natural CRN representations are ongoing challenges. We introduce Reactamole, a higher-level language for deterministic CRNs that utilizes the functional reactive programming (FRP) paradigm to represent CRNs as a reactive dataflow network. Reactamole equates a CRN with a functional reactive program, implementing the key primitives of the FRP paradigm directly as CRNs. The functional nature of Reactamole makes reasoning about molecular programs easier, and its strong static typing allows us to ensure that a CRN is well-formed by virtue of being well-typed. In this paper, we describe the design of Reactamole and how we use CRNs to represent the common datatypes and operations found in FRP. We also demonstrate the potential of this functional reactive approach to molecular programming by giving an extended example where a CRN is constructed using FRP to modulate and demodulate an amplitude modulated signal

ISCB Ebola Award for Important Future Research on the Computational Biology of Ebola Virus

Speed is of the essence in combating Ebola; thus, computational approaches should form a significant component of Ebola research. As for the development of any modern drug, computational biology is uniquely positioned to contribute through comparative analysis of the genome sequences of Ebola strains as well as 3-D protein modeling. Other computational approaches to Ebola may include large-scale docking studies of Ebola proteins with human proteins and with small-molecule libraries, computational modeling of the spread of the virus, computational mining of the Ebola literature, and creation of a curated Ebola database. Taken together, such computational efforts could significantly accelerate traditional scientific approaches. In recognition of the need for important and immediate solutions from the field of computational biology against Ebola, the International Society for Computational Biology (ISCB) announces a prize for an important computational advance in fighting the Ebola virus. ISCB will confer the ISCB Fight against Ebola Award, along with a prize of US$2,000, at its July 2016 annual meeting (ISCB Intelligent Systems for Molecular Biology [ISMB] 2016, Orlando, Florida)

Predicting the effects of climate change on water yield and forest production in the northeastern United States

Rapid and simultaneous changes in temperature, precipitation and the atmospheric concentration of CO2 are predicted to occur over the next century. Simple, well-validated models of ecosystem function are required to predict the effects of these changes. This paper describes an improved version of a forest carbon and water balance model (PnET-II) and the application of the model to predict stand- and regional-level effects of changes in temperature, precipitation and atmospheric CO2 concentration. PnET-II is a simple, generalized, monthly time-step model of water and carbon balances (gross and net) driven by nitrogen availability as expressed through foliar N concentration. Improvements from the original model include a complete carbon balance and improvements in the prediction of canopy phenology, as well as in the computation of canopy structure and photosynthesis. The model was parameterized and run for 4 forest/site combinations and validated against available data for water yield, gross and net carbon exchange and biomass production. The validation exercise suggests that the determination of actual water availability to stands and the occurrence or non-occurrence of soil-based water stress are critical to accurate modeling of forest net primary production (NPP) and net ecosystem production (NEP). The model was then run for the entire NewEngland/New York (USA) region using a 1 km resolution geographic information system. Predicted long-term NEP ranged from -85 to +275 g C m-2 yr-1 for the 4 forest/site combinations, and from -150 to 350 g C m-2 yr-1 for the region, with a regional average of 76 g C m-2 yr-1. A combination of increased temperature (+6*C), decreased precipitation (-15%) and increased water use efficiency (2x, due to doubling of CO2) resulted generally in increases in NPP and decreases in water yield over the region

Comprehensive User Engagement Sites (CUES) in Philadelphia: A Constructive Proposal

This paper is a study about Philadelphia’s comprehensive user engagement sites (CUESs) as the authors address and examine issues related to the upcoming implementation of a CUES while seeking solutions for its disputed questions and plans. Beginning with the federal drug schedules, the authors visit some of the medical and public health issues vis-à-vis safe injection facilities (SIFs). Insite, a successful Canadian SIF, has been thoroughly researched as it represents a paradigm for which a Philadelphia CUES can expand upon. Also, the existing criticisms against SIFs are revisited while critically unpackaged and responded to in favor of the establishment. In the main section, the authors propose the layout and services of the upcoming CUES, much of which would be in congruent to Vancouver’s Insite. On the other hand, the CUES would be distinct from Insite, as the authors emphasize, in that it will offer an information center run by individuals in recovery and place additional emphasis on early education for young healthcare professionals by providing them a platform to work at the site. The paper will also briefly investigate the implementation of a CUES site under an ethical scope of the Harm Reduction Theory. Lastly, the authors recommend some strategic plans that the Philadelphia City government may consider employing at this crucial stage

New genetic loci implicated in fasting glucose homeostasis and their impact on type 2 diabetes risk.

Levels of circulating glucose are tightly regulated. To identify new loci influencing glycemic traits, we performed meta-analyses of 21 genome-wide association studies informative for fasting glucose, fasting insulin and indices of beta-cell function (HOMA-B) and insulin resistance (HOMA-IR) in up to 46,186 nondiabetic participants. Follow-up of 25 loci in up to 76,558 additional subjects identified 16 loci associated with fasting glucose and HOMA-B and two loci associated with fasting insulin and HOMA-IR. These include nine loci newly associated with fasting glucose (in or near ADCY5, MADD, ADRA2A, CRY2, FADS1, GLIS3, SLC2A2, PROX1 and C2CD4B) and one influencing fasting insulin and HOMA-IR (near IGF1). We also demonstrated association of ADCY5, PROX1, GCK, GCKR and DGKB-TMEM195 with type 2 diabetes. Within these loci, likely biological candidate genes influence signal transduction, cell proliferation, development, glucose-sensing and circadian regulation. Our results demonstrate that genetic studies of glycemic traits can identify type 2 diabetes risk loci, as well as loci containing gene variants that are associated with a modest elevation in glucose levels but are not associated with overt diabetes

Hundreds of variants clustered in genomic loci and biological pathways affect human height

Most common human traits and diseases have a polygenic pattern of inheritance: DNA sequence variants at many genetic loci influence the phenotype. Genome-wide association (GWA) studies have identified more than 600 variants associated with human traits, but these typically explain small fractions of phenotypic variation, raising questions about the use of further studies. Here, using 183,727 individuals, we show that hundreds of genetic variants, in at least 180 loci, influence adult height, a highly heritable and classic polygenic trait. The large number of loci reveals patterns with important implications for genetic studies of common human diseases and traits. First, the 180 loci are not random, but instead are enriched for genes that are connected in biological pathways (P = 0.016) and that underlie skeletal growth defects (P < 0.001). Second, the likely causal gene is often located near the most strongly associated variant: in 13 of 21 loci containing a known skeletal growth gene, that gene was closest to the associated variant. Third, at least 19 loci have multiple independently associated variants, suggesting that allelic heterogeneity is a frequent feature of polygenic traits, that comprehensive explorations of already-discovered loci should discover additional variants and that an appreciable fraction of associated loci may have been identified. Fourth, associated variants are enriched for likely functional effects on genes, being over-represented among variants that alter amino-acid structure of proteins and expression levels of nearby genes. Our data explain approximately 10% of the phenotypic variation in height, and we estimate that unidentified common variants of similar effect sizes would increase this figure to approximately 16% of phenotypic variation (approximately 20% of heritable variation). Although additional approaches are needed to dissect the genetic architecture of polygenic human traits fully, our findings indicate that GWA studies can identify large numbers of loci that implicate biologically relevant genes and pathways.

Genetic risk and a primary role for cell-mediated immune mechanisms in multiple sclerosis.

Multiple sclerosis is a common disease of the central nervous system in which the interplay between inflammatory and neurodegenerative processes typically results in intermittent neurological disturbance followed by progressive accumulation of disability. Epidemiological studies have shown that genetic factors are primarily responsible for the substantially increased frequency of the disease seen in the relatives of affected individuals, and systematic attempts to identify linkage in multiplex families have confirmed that variation within the major histocompatibility complex (MHC) exerts the greatest individual effect on risk. Modestly powered genome-wide association studies (GWAS) have enabled more than 20 additional risk loci to be identified and have shown that multiple variants exerting modest individual effects have a key role in disease susceptibility. Most of the genetic architecture underlying susceptibility to the disease remains to be defined and is anticipated to require the analysis of sample sizes that are beyond the numbers currently available to individual research groups. In a collaborative GWAS involving 9,772 cases of European descent collected by 23 research groups working in 15 different countries, we have replicated almost all of the previously suggested associations and identified at least a further 29 novel susceptibility loci. Within the MHC we have refined the identity of the HLA-DRB1 risk alleles and confirmed that variation in the HLA-A gene underlies the independent protective effect attributable to the class I region. Immunologically relevant genes are significantly overrepresented among those mapping close to the identified loci and particularly implicate T-helper-cell differentiation in the pathogenesis of multiple sclerosis

Spatial Variation among Lakes within Landscapes: Ecological Organization along Lake Chains

Although limnologists have long been interested in regional patterns in lake attributes, only recently have they considered lakes connected and organized across the landscape, rather than as spatially independent entities. Here we explore the spatial organization of lake districts through the concept of landscape position, a concept that considers lakes longitudinally along gradients of geomorphology and hydrology. We analyzed long-term chemical and biological data from nine lake chains (lakes in a series connected through surface or groundwater flow) from seven lake districts of diverse hydrologic and geomorphic settings across North America. Spatial patterns in lake variables driven by landscape position were surprisingly common across lake districts and across a wide range of variables. On the other hand, temporal patterns of lake variables, quantified using synchrony, the degree to which pairs of lakes exhibit similar dynamics through time, related to landscape position only for lake chains with lake water residence times that spanned a wide range and were generally long (close to or greater than 1 year). Highest synchrony of lakes within a lake chain occurred when lakes had short water residence times. Our results from both the spatial and temporal analyses suggest that certain features of the landscape position concept are robust enough to span a wide range of seemingly disparate lake types. The strong spatial patterns observed in this analysis, and some unexplained patterns, suggest the need to further study these scales and to continue to view lake ecosystems spatially, longitudinally, and broadly across the landscape.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/42434/1/10021-2-5-395_2n5p395.pd

All-codon scanning identifies p53 cancer rescue mutations

In vitro scanning mutagenesis strategies are valuable tools to identify critical residues in proteins and to generate proteins with modified properties. We describe the fast and simple All-Codon Scanning (ACS) strategy that creates a defined gene library wherein each individual codon within a specific target region is changed into all possible codons with only a single codon change per mutagenesis product. ACS is based on a multiplexed overlapping mutagenesis primer design that saturates only the targeted gene region with single codon changes. We have used ACS to produce single amino-acid changes in small and large regions of the human tumor suppressor protein p53 to identify single amino-acid substitutions that can restore activity to inactive p53 found in human cancers. Single-tube reactions were used to saturate defined 30-nt regions with all possible codon changes. The same technique was used in 20 parallel reactions to scan the 600-bp fragment encoding the entire p53 core domain. Identification of several novel p53 cancer rescue mutations demonstrated the utility of the ACS approach. ACS is a fast, simple and versatile method, which is useful for protein structure–function analyses and protein design or evolution problems

Predicting Positive p53 Cancer Rescue Regions Using Most Informative Positive (MIP) Active Learning

Many protein engineering problems involve finding mutations that produce proteins with a particular function. Computational active learning is an attractive approach to discover desired biological activities. Traditional active learning techniques have been optimized to iteratively improve classifier accuracy, not to quickly discover biologically significant results. We report here a novel active learning technique, Most Informative Positive (MIP), which is tailored to biological problems because it seeks novel and informative positive results. MIP active learning differs from traditional active learning methods in two ways: (1) it preferentially seeks Positive (functionally active) examples; and (2) it may be effectively extended to select gene regions suitable for high throughput combinatorial mutagenesis. We applied MIP to discover mutations in the tumor suppressor protein p53 that reactivate mutated p53 found in human cancers. This is an important biomedical goal because p53 mutants have been implicated in half of all human cancers, and restoring active p53 in tumors leads to tumor regression. MIP found Positive (cancer rescue) p53 mutants in silico using 33% fewer experiments than traditional non-MIP active learning, with only a minor decrease in classifier accuracy. Applying MIP to in vivo experimentation yielded immediate Positive results. Ten different p53 mutations found in human cancers were paired in silico with all possible single amino acid rescue mutations, from which MIP was used to select a Positive Region predicted to be enriched for p53 cancer rescue mutants. In vivo assays showed that the predicted Positive Region: (1) had significantly more (p<0.01) new strong cancer rescue mutants than control regions (Negative, and non-MIP active learning); (2) had slightly more new strong cancer rescue mutants than an Expert region selected for purely biological considerations; and (3) rescued for the first time the previously unrescuable p53 cancer mutant P152L