A machine learning approach for the prediction of protein surface loop flexibility

Proteins often undergo conformational changes when binding to each other. A major fraction of backbone conformational changes involves motion on the protein surface, particularly in loops. Accounting for the motion of protein surface loops represents a challenge for protein-protein docking algorithms. A first step in addressing this challenge is to distinguish protein surface loops that are likely to undergo backbone conformational changes upon protein-protein binding (mobile loops) from those that are not (stationary loops). In this study, we developed a machine learning strategy based on support vector machines (SVMs). Our SVM uses three features of loop residues in the unbound protein structures-Ramachandran angles, crystallographic B-factors, and relative accessible surface area-to distinguish mobile loops from stationary ones. This method yields an average prediction accuracy of 75.3% compared with a random prediction accuracy of 50%, and an average of 0.79 area under the receiver operating characteristic (ROC) curve using cross-validation. Testing the method on an independent dataset, we obtained a prediction accuracy of 70.5%. Finally, we applied the method to 11 complexes that involve members from the Ras superfamily and achieved prediction accuracy of 92.8% for the Ras superfamily proteins and 74.4% for their binding partners.     

Next generation genome-wide association tool: design and coverage of a high-throughput European-optimized SNP array

The success of genome-wide association studies has paralleled the development of efficient genotyping technologies. We describe the development of a next-generation microarray based on the new highly-efficient Affymetrix Axiom genotyping technology that we are using to genotype individuals of European ancestry from the Kaiser Permanente Research Program on Genes, Environment and Health (RPGEH). The array contains 674,517 SNPs, and provides excellent genome-wide as well as gene-based and candidate-SNP coverage. Coverage was calculated using an approach based on imputation and cross validation. Preliminary results for the first 80,301 saliva-derived DNA samples from the RPGEH demonstrate very high quality genotypes, with sample success rates above 94% and over 98% of successful samples having SNP call rates exceeding 98%. At steady state, we have produced 462 million genotypes per week for each Axiom system. The new array provides a valuable addition to the repertoire of tools for large scale genome-wide association studies.     

Gene expression profiling-based identification of CD28 and PI3K as new biomarkers for chronic graft-versus-host disease

Chronic graft-versus-host disease (cGVHD) is a major cause of morbidity and mortality after allogeneic hematopoietic stem cell transplantation. Currently, no reliable biomarkers are available to predict the onset or progression of cGVHD. Therefore, in this study, we collected peripheral blood mononuclear cells from four patients with cGVHD and four ones with non-GVHD after hematopoietic stem cell transplantation and employed Affymetrix GeneChip Human U133 Plus 2.0 microarrays to screen the genes differentially expressed in cGVHD versus non-GVHD groups, with the aim to identify potential clinical biomarkers to predict cGVHD risk or progression. Microarray analysis demonstrated that the expression of 3180 genes changed significantly in cGVHD versus non-GVHD, with 879 genes upregulated and 2301 genes downregulated. Among them we chose CD28 and PI3K as candidates for further verification. Flow cytometry and quantitative real-time polymerase chain reaction analysis confirmed the significant upregulation of CD28 and PI3K in samples from patients with cGVHD compared with patients with non-GVHD, respectively. In conclusion, our study suggested that the upregulation of CD28 and PI3K contributed to the onset and progression of cGVHD and provided evidence that CD28 and PI3K may serve as promising biomarkers for cGVHD.     

Metabolic capabilities and systems fluctuations in Haloarcula marismortui revealed by integrative genomics and proteomics analyses

The 1310 Haloarcula marismortui proteins identified from mid-log and late-log phase soluble and membrane proteomes were analyzed in metabolic and cellular process networks to predict the available systems and systems fluctuations upon environmental stresses. When the connected metabolic reactions with identified proteins were examined, the availability of a number of metabolic pathways and a highly connected amino acid metabolic network were revealed. Quantitative spectral count analyses suggested 300 or more proteins might have expression changes in late-log phase. Among these, integrative network analyses indicated approximately 106 were metabolic proteins that might have growth-phase dependent changes. Interestingly, a large proportion of proteins in affected biomodules had the same trend of changes in spectral counts. Disregard the magnitude of changes, we had successfully predicted and validated the expression changes of nine genes including the rimK, gltCP, rrnAC0132, and argC in lysine biosynthesis pathway which were downregulated in late-log phase. This study had not only revealed the expressed proteins but also the availability of biological systems in two growth phases, systems level changes in response to the stresses in late-log phase, cellular locations of identified proteins, and the likely regulated genes to facilitate further analyses in the postgenomic era.     

Structural variation in two human genomes mapped at single-nucleotide resolution by whole genome de novo assembly

Here we use whole-genome de novo assembly of second-generation sequencing reads to map structural variation (SV) in an Asian genome and an African genome. Our approach identifies small- and intermediate-size homozygous variants (1-50 kb) including insertions, deletions, inversions and their precise breakpoints, and in contrast to other methods, can resolve complex rearrangements. In total, we identified 277,243 SVs ranging in length from 1-23 kb. Validation using computational and experimental methods suggests that we achieve overall <6% false-positive rate and <10% false-negative rate in genomic regions that can be assembled, which outperforms other methods. Analysis of the SVs in the genomes of 106 individuals sequenced as part of the 1000 Genomes Project suggests that SVs account for a greater fraction of the diversity between individuals than do single-nucleotide polymorphisms (SNPs). These findings demonstrate that whole-genome de novo assembly is a feasible approach to deriving more comprehensive maps of genetic variation.     

DroPhEA: Drosophila phenotype enrichment analysis for insect functional genomics

SUMMARY: DroPhEA is a core module of a web application that facilitates research in insect functional genomics through enrichment analysis on mutant phenotypes of fruit fly (Drosophila melanogaster). The phenotypes investigated in the analyses can be predefined by FlyBase or customized by users. DroPhEA allows users to specify mutation or ortholog types, displays enriched term results in a hierarchical structure and supports analyses on gene sets of all insect species with a fully sequenced genome.     

Graphene–Cellulose Paper Flexible Supercapacitors

A simple and scalable method to fabricate graphene‐cellulose paper (GCP) membranes is reported; these membranes exhibit great advantages as freestanding and binder‐free electrodes for flexible supercapacitors. The GCP electrode consists of a unique three‐dimensional interwoven structure of graphene nanosheets and cellulose fibers and has excellent mechanical flexibility, good specific capacitance and power performance, and excellent cyclic stability. The electrical conductivity of the GCP membrane shows high stability with a decrease of only 6% after being bent 1000 times. This flexible GCP electrode has a high capacitance per geometric area of 81 mF cm^?2, which is equivalent to a gravimetric capacitance of 120 F g^?1 of graphene, and retains >99% capacitance over 5000 cycles. Several types of flexible GCP‐based polymer supercapacitors with various architectures are assembled to meet the power‐energy requirements of typical flexible or printable electronics. Under highly flexible conditions, the supercapacitors show a high capacitance per geometric area of 46 mF cm^?2 for the complete devices. All the results demonstrate that polymer supercapacitors made using GCP membranes are versatile and may be used for flexible and portable micropower devices.     

Molecular characterization of acquired tolerance of tumor cells to picropodophyllin (PPP)

Background

Picropodophyllin (PPP) is a promising novel anti-neoplastic agent that efficiently kills tumor cells in vitro and causes tumor regression and increased survival in vivo. We have previously reported that PPP treatment induced moderate tolerance in two out of 10 cell lines only, and here report the acquired genomic and expression alterations associated with PPP selection over 1.5 years of treatment.

Methodology/Principal Findings

Copy number alterations monitored using metaphase and array-based comparative genomic hybridization analyses revealed largely overlapping alterations in parental and maximally tolerant cells. Gain/ amplification of the MYC and PVT1 loci in 8q24.21 were verified on the chromosome level. Abnormalities observed in connection to PPP treatment included regular gains and losses, as well as homozygous losses in 10q24.1-q24.2 and 12p12.3-p13.2 in one of the lines and amplification at 5q11.2 in the other. Abnormalities observed in both tolerant derivatives include amplification/gain of 5q11.2, gain of 11q12.1-q14.3 and gain of 13q33.3-qter. Using Nexus software analysis we combined the array-CGH data with data from gene expression profilings and identified genes that were altered in both inputs. A subset of genes identified as downregulated (ALDH1A3, ANXA1, TLR4 and RAB5A) or upregulated (COX6A1, NFIX, ME1, MAPK and TAP2) were validated by siRNA in the tolerant or parental cells to alter sensitivity to PPP and confirmed to alter sensitivity to PPP in further cell lines.

Conclusions

Long-term PPP selection lead to altered gene expression in PPP tolerant cells with increase as well as decrease of genes involved in cell death such as PTEN and BCL2. In addition, acquired genomic copy number alterations were observed that were often reflected by altered mRNA expression levels for genes in the same regions.     

Identification of a stable quantitative trait locus for percentage grains with white chalkiness in rice (Oryza sativa)

High chalkiness is a major problem in many rice-producing areas of the world, especially in hybrid rice (Oryza sativa L.) in China. We previously showed a major quantitative trait locus for the percentage of grains with white chalkiness (QTLqPGWC-8) in the interval G1149-R727 on chromosome 8 using a chromosome segment substitution line (CSSL). Here, we selected the line-CSSL50 harboring the QTLqPGWC-8 allele from the CSSLs derived from a cross between Asominori (as a recurrent parent) and IR24 (as a donor parent), which had higher percentage chalkiness, markedly different from that of Asominori. There were also significant differences in starch granules, appearance of amylose content (AAC) and milling qualities between Asominori and CSSL50, but not in grain size or thousand grain weight (TGW). The BC(4) F(2) and BC(4) F(3) populations from a cross between CSSL50 and Asominori were used for fine mapping of qPGWC-8. We narrowed down the location of this QTL to a 142 kb region between Indel markers 8G-7 and 8G-9. QTLqPGWC-8 accounted for 50.9% of the difference in PGWC between the parents. The markers tightly linked to qPGWC-8 should facilitate cloning of the gene underlying this QTL and will be of value for marker-assisted selection in breeding rice varieties with better grain quality.