How to choose templates for modeling of protein complexes: Insights from benchmarking template-based docking

Comparative docking is based on experimentally determined structures of protein-protein complexes (templates), following the paradigm that proteins with similar sequences and/or structures form similar complexes. Modeling utilizing structure similarity of target monomers to template complexes significantly expands structural coverage of the interactome. Template-based docking by structure alignment can be performed for the entire structures or by aligning targets to the bound interfaces of the experimentally determined complexes. Systematic benchmarking of docking protocols based on full and interface structure alignment showed that both protocols perform similarly, with top 1 docking success rate 26%. However, in terms of the models' quality, the interface-based docking performed marginally better. The interface-based docking is preferable when one would suspect a significant conformational change in the full protein structure upon binding, for example, a rearrangement of the domains in multidomain proteins. Importantly, if the same structure is selected as the top template by both full and interface alignment, the docking success rate increases 2-fold for both top 1 and top 10 predictions. Matching structural annotations of the target and template proteins for template detection, as a computationally less expensive alternative to structural alignment, did not improve the docking performance. Sophisticated remote sequence homology detection added templates to the pool of those identified by structure-based alignment, suggesting that for practical docking, the combination of the structure alignment protocols and the remote sequence homology detection may be useful in order to avoid potential flaws in generation of the structural templates library.     

SubcloneSeeker: a computational framework for reconstructing tumor clone structure for cancer variant interpretation and prioritization

Many tumors are composed of genetically divergent cell subpopulations. We report SubcloneSeeker, a package capable of exhaustive identification of subclone structures and evolutionary histories with bulk somatic variant allele frequency measurements from tumor biopsies. We present a statistical framework to elucidate whether specific sets of mutations are present within the same subclones, and the order in which they occur. We demonstrate how subclone reconstruction provides crucial information about tumorigenesis and relapse mechanisms; guides functional study by variant prioritization, and has the potential as a rational basis for informed therapeutic strategies for the patient. SubcloneSeeker is available at: https://github.com/yiq/SubcloneSeeker.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0443-x) contains supplementary material, which is available to authorized users.     

Escherichia coli O157 Somatic Antigen Is Present in an Isolate of E. fergusonii

A bacterium that tested positive with antibodies specific for Escherichia coli O157 was isolated from beef during routine screening procedures. The bacterium was identified as E. fergusonii by biochemical testing and partial sequencing of 16S rRNA. The isolate was tested for the presence of genes encoding Shiga toxins, the E. coli attaching and effacing factor, enterohemolysin, and the O157 O antigen. The isolate tested negative for Shiga toxins and other E. coli O157 virulence markers but was found to harbor the genes encoding the O157 antigen. These results suggest genetic transfer of the O antigen gene cluster between E. coli O157:H7 and E. fergusonii.     

Evidence for over-dispersion in the distribution of clinical malaria episodes in children

Background

It may be assumed that patterns of clinical malaria in children of similar age under the same level of exposure would follow a Poisson distribution with no over-dispersion. Longitudinal studies that have been conducted over many years suggest that some children may experience more episodes of clinical malaria than would be expected. The aim of this study was to identify this group of children and investigate possible causes for this increased susceptibility.

Methodology and Principal Findings

Using Poisson regression, we chose a group of children whom we designated as ‘more susceptible’ to malaria from 373 children under 10 years of age who were followed up for between 3 to 5 years from 1998–2003. About 21% of the children were categorized as ‘more susceptible’ and although they contributed only 23% of the person-time of follow-up, they experienced 55% of total clinical malaria episodes. Children that were parasite negative at all cross-sectional survey were less likely to belong to this group [AOR = 0.09, (95% CI: 0.14–0.61), p = 0.001].

Conclusions and Significance

The pattern of clinical malaria episodes follows a negative binomial distribution. Use of lack of a clinical malaria episode in a certain time period as endpoints for intervention or immunological studies may not adequately distinguish groups who are more or less immune. It may be useful in such studies, in addition to the usual endpoint of the time to first episode, to include end points which take into account the total number of clinical episodes experienced per child.     

HeatMapper: powerful combined visualization of gene expression profile correlations, genotypes, phenotypes and sample characteristics

Background  

Accurate interpretation of data obtained by unsupervised analysis of large scale expression profiling studies is currently frequently performed by visually combining sample-gene heatmaps and sample characteristics. This method is not optimal for comparing individual samples or groups of samples. Here, we describe an approach to visually integrate the results of unsupervised and supervised cluster analysis using a correlation plot and additional sample metadata.     

Shiga toxin-producing Escherichia coli in sheep and pre-slaughter lambs in eastern Australia

Sheep and lambs from 14 farms in southern Queensland and one from central New South Wales were surveyed to determine the prevalence of Shiga toxin-producing Escherichia coli (STEC). STEC, isolated from 45% of 144 sheep faeces collected on the farms and 36% of 72 lamb faeces from abattoir yards, were tested for the presence of genes encoding virulence factors. Most (64%) of the 117 STEC isolates contained Shiga toxin 1 and 2 genes, 22% contained those encoding Shiga toxin 1, and 14% contained genes encoding Shiga toxin 2. The genes encoding the E. coli attaching and effacing factor were present in 2.6% of STEC and 26% contained the enterohaemolysin gene. The isolates that contained the E. coli attaching and effacing gene were serotype O157:H. This study has shown that STEC are widely distributed in eastern Australian sheep and lambs and are shed in their faeces prior to slaughter. Thus, there is potential for contamination of carcasses and entry of STEC into the human food chain.     

Nucleotide sequence of the nifH gene coding for nitrogen reductase in the acetic acid bacterium Acetobacter diazotrophicus

The nifH gene sequence of the nitrogen-fixing bacterium Acetobacter diazotrophicus was determined with the use of the polymerase chain reaction and universal degenerate oligonucleotide primers. The gene shows highest pair-wise similarity to the nifH gene of Azospirillum brasilense . The phylogenetic relationships of the nifH gene sequences were compared with those inferred from 16S rRNA gene sequences. Knowledge of the sequence of the nifH gene contributes to the growing database of nifH gene sequences, and will allow the detection of Acet. diazotrophicus from environmental samples with nifH gene-based primers.