Projected expansion of Trichodesmium's geographical distribution and increase in growth potential in response to climate change.

Estimates of marine N2 fixation range from 52 to 73 Tg N/year, of which we calculate up to 84% is from Trichodesmium based on previous measurements of nifH gene abundance and our new model of Trichodesmium growth. Here, we assess the likely effects of four major climate change-related abiotic factors on the spatiotemporal distribution and growth potential of Trichodesmium for the last glacial maximum (LGM), the present (2006-2015) and the end of this century (2100) by mapping our model of Trichodesmium growth onto inferred global surface ocean fields of pCO2 , temperature, light and Fe. We conclude that growth rate was severely limited by low pCO2 at the LGM, that current pCO2 levels do not significantly limit Trichodesmium growth and thus, the potential for enhanced growth from future increases in CO2 is small. We also found that the area of the ocean where sea surface temperatures (SST) are within Trichodesmium's thermal niche increased by 32% from the LGM to present, but further increases in SST due to continued global warming will reduce this area by 9%. However, the range reduction at the equator is likely to be offset by enhanced growth associated with expansion of regions with optimal or near optimal Fe and light availability. Between now and 2100, the ocean area of optimal SST and irradiance is projected to increase by 7%, and the ocean area of optimal SST, irradiance and iron is projected to increase by 173%. Given the major contribution of this keystone species to annual N2 fixation and thus pelagic ecology, biogeochemistry and CO2 sequestration, the projected increase in the geographical range for optimal growth could provide a negative feedback to increasing atmospheric CO2 concentrations.

One motif to bind them: A small-XXX-small motif affects transmembrane domain 1 oligomerization, function, localization, and cross-talk between two yeast GPCRs.

G protein-coupled receptors (GPCRs) are the largest family of cell-surface receptors in mammals and facilitate a range of physiological responses triggered by a variety of ligands. GPCRs were thought to function as monomers, however it is now accepted that GPCR homo- and hetero-oligomers also exist and influence receptor properties. The Schizosaccharomyces pombe GPCR Mam2 is a pheromone-sensing receptor involved in mating and has previously been shown to form oligomers in vivo. The first transmembrane domain (TMD) of Mam2 contains a small-XXX-small motif, overrepresented in membrane proteins and well-known for promoting helix-helix interactions. An ortholog of Mam2 in Saccharomyces cerevisiae, Ste2, contains an analogous small-XXX-small motif which has been shown to contribute to receptor homo-oligomerization, localization and function. Here we have used experimental and computational techniques to characterize the role of the small-XXX-small motif in function and assembly of Mam2 for the first time. We find that disruption of the motif via mutagenesis leads to reduction of Mam2 TMD1 homo-oligomerization and pheromone-responsive cellular signaling of the full-length protein. It also impairs correct targeting to the plasma membrane. Mutation of the analogous motif in Ste2 yielded similar results, suggesting a conserved mechanism for assembly. Using co-expression of the two fungal receptors in conjunction with computational models, we demonstrate a functional change in G protein specificity and propose that this is brought about through hetero-dimeric interactions of Mam2 with Ste2 via the complementary small-XXX-small motifs. This highlights the potential of these motifs to affect a range of properties that can be investigated in other GPCRs.

Enhancing human islet transplantation by localized release of trophic factors from PLG scaffolds.

Islet transplantation represents a potential cure for type 1 diabetes, yet the clinical approach of intrahepatic delivery is limited by the microenvironment. Microporous scaffolds enable extrahepatic transplantation, and the microenvironment can be designed to enhance islet engraftment and function. We investigated localized trophic factor delivery in a xenogeneic human islet to mouse model of islet transplantation. Double emulsion microspheres containing exendin-4 (Ex4) or insulin-like growth factor-1 (IGF-1) were incorporated into a layered scaffold design consisting of porous outer layers for islet transplantation and a center layer for sustained factor release. Protein encapsulation and release were dependent on both the polymer concentration and the identity of the protein. Proteins retained bioactivity upon release from scaffolds in vitro. A minimal human islet mass transplanted on Ex4-releasing scaffolds demonstrated significant improvement and prolongation of graft function relative to blank scaffolds carrying no protein, and the release profile significantly impacted the duration over which the graft functioned. Ex4-releasing scaffolds enabled better glycemic control in animals subjected to an intraperitoneal glucose tolerance test. Scaffolds releasing IGF-1 lowered blood glucose levels, yet the reduction was insufficient to achieve euglycemia. Ex4-delivering scaffolds provide an extrahepatic transplantation site for modulating the islet microenvironment to enhance islet function posttransplant.

UK catchment nutrient loads 1993-2003, a new approach using harmonised monitoring scheme data: temporal changes, geographical distribution, limiting nutrients and loads to coastal waters.

The work provides robust estimates of nutrient loads (nitrate and phosphate) from all UK catchments: as required by the Water Framework Directive to monitor catchments' health, and to inform management of these environments. To calculate nutrient loads, data for nutrient concentrations and water flow are combined. In the UK, flow data are typically available at hourly intervals at more than 1300 gauging stations but concentration data are collected less frequently (roughly weekly) and at fewer locations (about 280). The sparseness of the concentration data limits the occasions for which load can be calculated, so a mathematical model was derived which was used to interpolate the concentrations between measurements. The model's parameters provide useful information about the annual nutrient concentration cycles within any catchment, and permitted improved estimates of both the annual loads of N and P, and of the N : P ratios, from mainland UK catchments. Data from 1993-2003 showed nitrate loads from UK catchments were generally constant, while orthophosphate loads generally declined. N : P ratios suggested that most catchments in the north and west of the UK were potentially P-limited although a few were potentially N-limited, while many in central and eastern UK oscillated seasonally between N and P limitation. Knowledge of the nutrient which is potentially limiting to biological productivity is a key factor for management of a catchment's nutrient loads. Calculations of nutrient export loads to coastal regions showed that UK catchments contributed only about 16.5% of total fluvial loads of nitrate to the North Sea, or about 3% of the total N loads when inputs from the Atlantic were included. Orthophosphate loads from the UK catchments into the North Sea were only 1.7% of the total P inputs from rivers and the Atlantic but did not include riverine inputs of P adsorbed to particles.

Assessing the effect of dynamics on the closed-loop protein-folding hypothesis.

The closed-loop (loop-n-lock) hypothesis of protein folding suggests that loops of about 25 residues, closed through interactions between the loop ends (locks), play an important role in protein structure. Coarse-grain elastic network simulations, and examination of loop lengths in a diverse set of proteins, each supports a bias towards loops of close to 25 residues in length between residues of high stability. Previous studies have established a correlation between total contact distance (TCD), a metric of sequence distances between contacting residues (cf. contact order), and the log-folding rate of a protein. In a set of 43 proteins, we identify an improved correlation (r(2) = 0.76), when the metric is restricted to residues contacting the locks, compared to the equivalent result when all residues are considered (r(2) = 0.65). This provides qualified support for the hypothesis, albeit with an increased emphasis upon the importance of a much larger set of residues surrounding the locks. Evidence of a similar-sized protein core/extended nucleus (with significant overlap) was obtained from TCD calculations in which residues were successively eliminated according to their hydrophobicity and connectivity, and from molecular dynamics simulations. Our results suggest that while folding is determined by a subset of residues that can be predicted by application of the closed-loop hypothesis, the original hypothesis is too simplistic; efficient protein folding is dependent on a considerably larger subset of residues than those involved in lock formation.

Do plants contain g protein-coupled receptors?

Whether G protein-coupled receptors (GPCRs) exist in plants is a fundamental biological question. Interest in deorphanizing new GPCRs arises because of their importance in signaling. Within plants, this is controversial, as genome analysis has identified 56 putative GPCRs, including G protein-coupled receptor1 (GCR1), which is reportedly a remote homolog to class A, B, and E GPCRs. Of these, GCR2 is not a GPCR; more recently, it has been proposed that none are, not even GCR1. We have addressed this disparity between genome analysis and biological evidence through a structural bioinformatics study, involving fold recognition methods, from which only GCR1 emerges as a strong candidate. To further probe GCR1, we have developed a novel helix-alignment method, which has been benchmarked against the class A-class B-class F GPCR alignments. In addition, we have presented a mutually consistent set of alignments of GCR1 homologs to class A, class B, and class F GPCRs and shown that GCR1 is closer to class A and/or class B GPCRs than class A, class B, or class F GPCRs are to each other. To further probe GCR1, we have aligned transmembrane helix 3 of GCR1 to each of the six GPCR classes. Variability comparisons provide additional evidence that GCR1 homologs have the GPCR fold. From the alignments and a GCR1 comparative model, we have identified motifs that are common to GCR1, class A, B, and E GPCRs. We discuss the possibilities that emerge from this controversial evidence that GCR1 has a GPCR fold.

Similarity between class A and class B G-protein-coupled receptors exemplified through calcitonin gene-related peptide receptor modelling and mutagenesis studies.

Modelling class B G-protein-coupled receptors (GPCRs) using class A GPCR structural templates is difficult due to lack of homology. The plant GPCR, GCR1, has homology to both class A and class B GPCRs. We have used this to generate a class A-class B alignment, and by incorporating maximum lagged correlation of entropy and hydrophobicity into a consensus score, we have been able to align receptor transmembrane regions. We have applied this analysis to generate active and inactive homology models of the class B calcitonin gene-related peptide (CGRP) receptor, and have supported it with site-directed mutagenesis data using 122 CGRP receptor residues and 144 published mutagenesis results on other class B GPCRs. The variation of sequence variability with structure, the analysis of polarity violations, the alignment of group-conserved residues and the mutagenesis results at 27 key positions were particularly informative in distinguishing between the proposed and plausible alternative alignments. Furthermore, we have been able to associate the key molecular features of the class B GPCR signalling machinery with their class A counterparts for the first time. These include the [K/R]KLH motif in intracellular loop 1, [I/L]xxxL and KxxK at the intracellular end of TM5 and TM6, the NPXXY/VAVLY motif on TM7 and small group-conserved residues in TM1, TM2, TM3 and TM7. The equivalent of the class A DRY motif is proposed to involve Arg(2.39), His(2.43) and Glu(3.46), which makes a polar lock with T(6.37). These alignments and models provide useful tools for understanding class B GPCR function.

Motif effects in Affymetrix GeneChips seriously affect probe intensities.

An Affymetrix GeneChip consists of an array of hundreds of thousands of probes (each a sequence of 25 bases) with the probe values being used to infer the extent to which genes are expressed in the biological material under investigation. In this article, we demonstrate that these probe values are also strongly influenced by their precise base sequence. We use data from >28 000 CEL files relating to 10 different Affymetrix GeneChip platforms and involving nearly 1000 experiments. Our results confirm known effects (those due to the T7-primer and the formation of G-quadruplexes) but reveal other effects. We show that there can be huge variations from one experiment to another, and that there may also be sizeable disparities between batches within an experiment and between CEL files within a batch.

G-protein-coupled receptor dynamics: dimerization and activation models compared with experiment.

Our previously derived models of the active state of the ß2-adrenergic receptor are compared with recently published X-ray crystallographic structures of activated GPCRs (G-protein-coupled receptors). These molecular dynamics-based models using experimental data derived from biophysical experiments on activation were used to restrain the receptor to an active state that gave high enrichment for agonists in virtual screening. The ß2-adrenergic receptor active model and X-ray structures are in good agreement over both the transmembrane region and the orthosteric binding site, although in some regions the active model is more similar to the active rhodopsin X-ray structures. The general features of the microswitches were well reproduced, but with minor differences, partly because of the unexpected X-ray results for the rotamer toggle switch. In addition, most of the interacting residues between the receptor and the G-protein were identified. This analysis of the modelling has also given important additional insight into GPCR dimerization: re-analysis of results on photoaffinity analogues of rhodopsin provided additional evidence that TM4 (transmembrane helix 4) resides at the dimer interface and that ligands such as bivalent ligands may pass between the mobile helices. A comparison, and discussion, is also carried out between the use of implicit and explicit solvent for active-state modelling.

Normalized Affymetrix expression data are biased by G-quadruplex formation.

Probes with runs of four or more guanines (G-stacks) in their sequences can exhibit a level of hybridization that is unrelated to the expression levels of the mRNA that they are intended to measure. This is most likely caused by the formation of G-quadruplexes, where inter-probe guanines form Hoogsteen hydrogen bonds, which probes with G-stacks are capable of forming. We demonstrate that for a specific microarray data set using the Human HG_U133A Affymetrix GeneChip and RMA normalization there is significant bias in the expression levels, the fold change and the correlations between expression levels. These effects grow more pronounced as the number of G-stack probes in a probe set increases. Approximately 14% of the probe sets are directly affected. The analysis was repeated for a number of other normalization pipelines and two, FARMS and PLIER, minimized the bias to some extent. We estimate that ∼15% of the data sets deposited in the GEO database are susceptible to the effect. The inclusion of G-stack probes in the affected data sets can bias key parameters used in the selection and clustering of genes. The elimination of these probes from any analysis in such affected data sets outweighs the increase of noise in the signal.

Antisense suppression of the small chloroplast protein CP12 in tobacco: a transcriptional viewpoint.

The chloroplast protein CP12 forms a multi-enzyme complex with the Calvin-Benson cycle enzymes phosphoribulokinase (PRK) and NADP-glyceraldehyde-3-phosphate dehydrogenase (GAPDH). PRK and GAPDH are inactivated when present in this complex, a process shown in vitro to be dependent upon oxidized CP12. Recently we reported on the importance of CP12 in vivo to higher plant metabolism using antisense suppression of CP12 in tobacco (Nicotiana tabacum). Our results indicated that while only minor changes in photosynthetic carbon fixation and in PRK and GAPDH activities were observed, striking changes in growth rates and morphology were seen. In this article we present data on the transcriptional changes observed in one of the antisense lines and we discuss the major findings in light of the metabolic phenotype described.

Advancing islet transplantation: from engraftment to the immune response.

The promise and progress of islet transplantation for treating type 1 diabetes has been challenged by obstacles to patient accessibility and long-term graft function that may be overcome by integrating emerging technologies in biomaterials, drug delivery and immunomodulation. The hepatic microenvironment and traditional systemic immunosuppression stress the vulnerable islets and contribute to the limited success of transplantation. Locally delivering extracellular matrix proteins and trophic factors can enhance transplantation at extrahepatic sites by promoting islet engraftment, revascularisation and long-term function while avoiding unintended systemic effects. Cell- and cytokine-based therapies for immune cell recruitment and reprogramming can inhibit local and systemic immune system activation that normally attacks transplanted islets. Combined with antigen-specific immunotherapies, states of operational tolerance may be achievable, reducing or eliminating the long-term pharmaceutical burden. Integration of these technologies to enhance engraftment and combat rejection may help to advance the therapeutic efficacy and availability of islet transplantation.

The detection of blur in affymetrix GeneChips.

High correlations were obtained between probes in seemingly unrelated probe sets, following an examination of the data from thousands of Affymetrix GeneChips. Investigation revealed that these unexpected correlations were between probes that were adjacent to high-valued probes. Using carefully selected probes, together with simple linear models, the extent of blur has been measured for each CEL file. The cause is shown to be attributable to poorly performing scanners. Blur can result in the doubling of the values of thousands of probes. This in turn can lead to the doubling of the expression level for hundreds of probe sets.

Closed loop folding units from structural alignments: experimental foldons revisited.

Nonoverlapping closed loops of around 25-35 amino acids formed via nonlocal interactions at the loop ends have been proposed as an important unit of protein structure. This hypothesis is significant as such short loops can fold quickly and so would not be bound by the Leventhal paradox, giving insight into the possible nature of the funnel in protein folding. Previously, these closed loops have been identified either by sequence analysis (conservation and autocorrelation) or studies of the geometry of individual proteins. Given the potential significance of the closed loop hypothesis, we have explored a new strategy for determining closed loops from the insertions identified by the structural alignment of proteins sharing the same overall fold. We determined the locations of the closed loops in 37 pairs of proteins and obtained excellent agreement with previously published closed loops. The relevance of NMR structures to closed loop determination is briefly discussed. For cytochrome c, cytochrome b(562) and triosephophate isomerase, independent folding units have been determined on the basis of hydrogen exchange experiments and misincorporation proton-alkyl exchange experiments. The correspondence between these experimentally derived foldons and the theoretically derived closed loops indicates that the closed loop hypothesis may provide a useful framework for analyzing such experimental data.

A Comparative Study of the Impact of G-Stack Probes on Various Affymetrix GeneChips of Mammalia.

We have previously discovered that probes containing runs of four or more contiguous guanines are not reliable for measuring gene expression in the Human HG_U133A Affymetrix GeneChip data. These probes are not correlated with other members of their probe set, but they are correlated with each other. We now extend our analysis to different 3' GeneChip designs of mouse, rat, and human. We find that, in all these chip designs, the G-stack probes (probes with a run of exactly four consecutive guanines) are correlated highly with each other, indicating that such probes are not reliable measures of gene expression in mammalian studies. Furthermore, there is no specific position of G-stack where the correlation is highest in all the chips. We also find that the latest designs of rat and mouse chips have significantly fewer G-stack probes compared to their predecessors, whereas there has not been a similar reduction in G-stack density across the changes in human chips. Moreover, we find significant changes in RMA values (after removing G-stack probes) as the number of G-stack probes increases.

Identifying the impact of G-quadruplexes on Affymetrix 3' arrays using cloud computing.

A tetramer quadruplex structure is formed by four parallel strands of DNA/ RNA containing runs of guanine. These quadruplexes are able to form because guanine can Hoogsteen hydrogen bond to other guanines, and a tetrad of guanines can form a stable arrangement. Recently we have discovered that probes on Affymetrix GeneChips that contain runs of guanine do not measure gene expression reliably. We associate this finding with the likelihood that quadruplexes are forming on the surface of GeneChips. In order to cope with the rapidly expanding size of GeneChip array datasets in the public domain, we are exploring the use of cloud computing to replicate our experiments on 3' arrays to look at the effect of the location of G-spots (runs of guanines). Cloud computing is a recently introduced high-performance solution that takes advantage of the computational infrastructure of large organisations such as Amazon and Google. We expect that cloud computing will become widely adopted because it enables bioinformaticians to avoid capital expenditure on expensive computing resources and to only pay a cloud computing provider for what is used. Moreover, as well as financial efficiency, cloud computing is an ecologically-friendly technology, it enables efficient data-sharing and we expect it to be faster for development purposes. Here we propose the advantageous use of cloud computing to perform a large data-mining analysis of public domain 3' arrays.

Using surveys of Affymetrix GeneChips to study antisense expression.

We have used large surveys of Affymetrix GeneChip data in the public domain to conduct a study of antisense expression across diverse conditions. We derive correlations between groups of probes which map uniquely to the same exon in the antisense direction. When there are no probes assigned to an exon in the sense direction we find that many of the antisense groups fail to detect a coherent block of transcription. We find that only a minority of these groups contain coherent blocks of antisense expression suggesting transcription. We also derive correlations between groups of probes which map uniquely to the same exon in both sense and antisense direction. In some of these cases the locations of sense probes overlap with the antisense probes, and the sense and antisense probe intensities are correlated with each other. This configuration suggests the existence of a Natural Antisense Transcript (NAT) pair. We find the majority of such NAT pairs detected by GeneChips are formed by a transcript of an established gene and either an EST or an mRNA. In order to determine the exact antisense regulatory mechanism indicated by the correlation of sense probes with antisense probes, a further investigation is necessary for every particular case of interest. However, the analysis of microarray data has proved to be a good method to reconfirm known NATs, discover new ones, as well as to notice possible problems in the annotation of antisense transcripts.

On the causes of outliers in Affymetrix GeneChip data.

We describe various types of outliers seen in Affymetrix GeneChip data. We have been able to utilise the data in the Gene Expression Omnibus to screen GeneChips across a range of scales, from single probes, to spatially adjacent fractions of arrays, to whole arrays, to whole experiments. In this review we describe a number of causes for why some reported intensities might be misleading on GeneChips.

Probes containing runs of guanines provide insights into the biophysics and bioinformatics of Affymetrix GeneChips.

The reliable interpretation of Affymetrix GeneChip data is a multi-faceted problem. The interplay between biophysics, bioinformatics and mining of GeneChip surveys is leading to new insights into how best to analyse the data. Many of the molecular processes occurring on the surfaces of GeneChips result from the high surface density of probes. Interactions between neighbouring adjacent probes affect their rate and strength of hybridization to targets. Competing targets may hybridize to the same probe, and targets may partially bind to more than one probe. The formation of these partial hybrids results in a number of probes not reaching thermodynamic equilibrium during hybridization. Moreover, some targets fold up, or cross-hybridize to other targets. Furthermore, probes may fold and can undergo chemical saturation. There are also sequence-dependent differences in the rates of target desorption during the washing stage. Improvements in the mappings between probe sequence and biological databases are leading to more accurate gene expression profiles. Moreover, algorithms that combine the intensities of multiple probes into single measures of expression are increasingly dependent upon models of the hybridization processes occurring on GeneChips. The large repositories of GeneChip data can be searched for systematic effects across many experiments. This data mining has led to the discovery of a family of thousands of probes, which show correlated expression across thousands of GeneChip experiments. These probes contain runs of guanines, suggesting that G-quadruplexes are able to form on GeneChips. We discuss the impact of these structures on the interpretation of data from GeneChip experiments.

Reducing spatial flaws in oligonucleotide arrays by using neighborhood information.

We address the problem of detection and correction of spatial flaws in oligonucleotide microarrays. We present two similar procedures, of which one is intended solely for use with replicates and the other has wider applicability. By constructing a set of replicates, with one realistically flawed, we are able to examine the extent to which our procedures are capable of repairing the flaw. We find that, for this purpose, our procedures are superior to the existing 'Harshlight' procedure.