Enabling high-throughput data management for systems biology: the Bioinformatics Resource Manager

The Bioinformatics Resource Manager (BRM) is a software environment that provides the user with data management, retrieval and integration capabilities. Designed in collaboration with biologists, BRM simplifies mundane analysis tasks of merging microarray and proteomic data across platforms, facilitates integration of users' data with functional annotation and interaction data from public sources and provides connectivity to visual analytic tools through reformatting of the data for easy import or dynamic launching capability. BRM is developed using Java and other open-source technologies for free distribution. AVAILABILITY: BRM, sample data sets and a user manual can be downloaded from http://www.sysbio.org/dataresources/brm.stm.     

Bayesian spatial modeling of moose count data: increasing estimator efficiency and exploring ecological hypotheses

Moose management throughout much of Alaska and Canada relies on aerial count data, which are commonly collected using the geospatial population estimator (GSPE) protocol. The GSPE uses a model-based analytical approach and finite-population block kriging to estimate abundance from a collection of sampled survey units. Widespread implementation and well-documented analytical software have resulted in reliable estimates of moose abundance, density, and composition across a large proportion of their range. Analysis is conducted almost exclusively using the GSPE software, which fits a fixed model structure to data collected within a single year. The downside of this approach to analysis is that the fixed model structure is inefficient for estimation, leading to more field effort than would otherwise be necessary to achieve a desired level of estimator precision. We developed a more easily modified and flexible Bayesian spatial general additive model approach (BSG) that accommodates spatial and temporal covariates (e.g., habitat characteristics, trend), multiple survey events, prior information, and incomplete detection. Using a series of 6 GSPE surveys conducted in Yukon-Charley Rivers National Preserve, Alaska, USA, from 2003–2019, we established the equivalence of the 2 approaches under similar model structures. We then extended the BSG to demonstrate how a more comprehensive approach to analysis can affect estimator precision and be used to assess ecological relationships. The precision of annual abundance estimators from the BSG were improved by an average of 43% over those based on the standard GSPE analysis, highlighting the very real costs of assuming a fixed (i.e., suboptimal) model structure. The population increased at a rate of 2.3%/year (95% CrI = 0.8–3.8%), and the increase was largely explained by a parallel increase in wildfire extent (i.e., high quality moose habitat). These results suggest that our approach could be used to increase estimator efficiency or decrease future survey costs without any modifications to the basic protocol. While modification of the GSPE software is possible, practitioners may find the BSG approach more convenient for quickly developing model structures for a particular application, thereby allowing them to extract more information from existing and future datasets.     

Enhancing the Detection of Giardia duodenalis Cysts in Foods by Inertial Microfluidic Separation

The sensitivity and specificity of current Giardia cyst detection methods for foods are largely determined by the effectiveness of the elution, separation, and concentration methods used. The aim of these methods is to produce a final suspension with an adequate concentration of Giardia cysts for detection and a low concentration of interfering food debris. In the present study, a microfluidic device, which makes use of inertial separation, was designed and fabricated for the separation of Giardia cysts. A cyclical pumping platform and protocol was developed to concentrate 10-ml suspensions down to less than 1 ml. Tests involving Giardia duodenalis cysts and 1.90-μm microbeads in pure suspensions demonstrated the specificity of the microfluidic chip for cysts over smaller nonspecific particles. As the suspension cycled through the chip, a large number of beads were removed (70%) and the majority of the cysts were concentrated (82%). Subsequently, the microfluidic inertial separation chip was integrated into a method for the detection of G. duodenalis cysts from lettuce samples. The method greatly reduced the concentration of background debris in the final suspensions (10-fold reduction) in comparison to that obtained by a conventional method. The method also recovered an average of 68.4% of cysts from 25-g lettuce samples and had a limit of detection (LOD) of 38 cysts. While the recovery of cysts by inertial separation was slightly lower, and the LOD slightly higher, than with the conventional method, the sample analysis time was greatly reduced, as there were far fewer background food particles interfering with the detection of cysts by immunofluorescence microscopy.     

Survey of phosphorylation near drug binding sites in the Protein Data Bank (PDB) and their effects

While it is currently estimated that 40 to 50% of eukaryotic proteins are phosphorylated, little is known about the frequency and local effects of phosphorylation near pharmaceutical inhibitor binding sites. In this study, we investigated how frequently phosphorylation may affect the binding of drug inhibitors to target proteins. We examined the 453 non‐redundant structures of soluble mammalian drug target proteins bound to inhibitors currently available in the Protein Data Bank (PDB). We cross‐referenced these structures with phosphorylation data available from the PhosphoSitePlus database. Three hundred twenty‐two of 453 (71%) of drug targets have evidence of phosphorylation that has been validated by multiple methods or labs. For 132 of 453 (29%) of those, the phosphorylation site is within 12 Å of the small molecule‐binding site, where it would likely alter small molecule binding affinity. We propose a framework for distinguishing between drug‐phosphorylation site interactions that are likely to alter the efficacy of drugs versus those that are not. In addition we highlight examples of well‐established drug targets, such as estrogen receptor alpha, for which phosphorylation may affect drug affinity and clinical efficacy. Our data suggest that phosphorylation may affect drug binding and efficacy for a significant fraction of drug target proteins. Proteins 2015; 83:25–36. © 2014 Wiley Periodicals, Inc.     

Glucosylceramides are critical for cell‐type differentiation and organogenesis,but not for cell viability in Arabidopsis

Glucosylceramides (GlcCer), glucose‐conjugated sphingolipids, are major components of the endomembrane system and plasma membrane in most eukaryotic cells. Yet the quantitative significance and cellular functions of GlcCer are not well characterized in plants and other multi‐organ eukaryotes. To address this, we examined Arabidopsis lines that were lacking or deficient in GlcCer by insertional disruption or by RNA interference (RNAi) suppression of the single gene for GlcCer synthase (GCS, At2g19880), the enzyme that catalyzes GlcCer synthesis. Null mutants for GCS (designated ‘gcs‐1’) were viable as seedlings, albeit strongly reduced in size, and failed to develop beyond the seedling stage. Heterozygous plants harboring the insertion allele exhibited reduced transmission through the male gametophyte. Undifferentiated calli generated from gcs‐1 seedlings and lacking GlcCer proliferated in a manner similar to calli from wild‐type plants. However, gcs‐1 calli, in contrast to wild‐type calli, were unable to develop organs on differentiation media. Consistent with a role for GlcCer in organ‐specific cell differentiation, calli from gcs‐1 mutants formed roots and leaves on media supplemented with the glucosylated sphingosine glucopsychosine, which was readily converted to GlcCer independent of GCS. Underlying these phenotypes, gcs‐1 cells had altered Golgi morphology and fewer cisternae per Golgi apparatus relative to wild‐type cells, indicative of protein trafficking defects. Despite seedling lethality in the null mutant, GCS RNAi suppression lines with ≤2% of wild‐type GlcCer levels were viable and fertile. Collectively, these results indicate that GlcCer are essential for cell‐type differentiation and organogenesis, and plant cells produce amounts of GlcCer in excess of that required for normal development.     

Sequential unfolding of the hemolysin two‐partner secretion domain from Proteus mirabilis

Protein secretion is a major contributor to Gram‐negative bacterial virulence. Type Vb or two‐partner secretion (TPS) pathways utilize a membrane bound β‐barrel B component (TpsB) to translocate large and predominantly virulent exoproteins (TpsA) through a nucleotide independent mechanism. We focused our studies on a truncated TpsA member termed hemolysin A (HpmA265), a structurally and functionally characterized TPS domain from Proteus mirabilis. Contrary to the expectation that the TPS domain of HpmA265 would denature in a single cooperative transition, we found that the unfolding follows a sequential model with three distinct transitions linking four states. The solvent inaccessible core of HpmA265 can be divided into two different regions. The C‐proximal region contains nonpolar residues and forms a prototypical hydrophobic core as found in globular proteins. The N‐proximal region of the solvent inaccessible core, however, contains polar residues. To understand the contributions of the hydrophobic and polar interiors to overall TPS domain stability, we conducted unfolding studies on HpmA265 and site‐specific mutants of HpmA265. By correlating the effect of individual site‐specific mutations with the sequential unfolding results we were able to divide the HpmA265 TPS domain into polar core, nonpolar core, and C‐terminal subdomains. Moreover, the unfolding studies provide quantitative evidence that the folding free energy for the polar core subdomain is more favorable than for the nonpolar core and C‐terminal subdomains. This study implicates the hydrogen bonds shared among these conserved internal residues as a primary means for stabilizing the N‐proximal polar core subdomain.