Inference algorithms for gene networks: a statistical mechanics analysis

The inference of gene regulatory networks from high throughput gene expression data is one of the major challenges in systems biology. This paper aims at analysing and comparing two different algorithmic approaches. The first approach uses pairwise correlations between regulated and regulating genes; the second one uses message-passing techniques for inferring activating and inhibiting regulatory interactions. The performance of these two algorithms can be analysed theoretically on well-defined test sets, using tools from the statistical physics of disordered systems like the replica method. We find that the second algorithm outperforms the first one since it takes into account collective effects of multiple regulators

A new species of western Atlantic lizardfish (Teleostei: Synodontidae: Synodus) and resurrection of Synodus bondi Fowler, 1939, as a valid species from the Caribbean with redescriptions of S. bondi, S. foetens (Linnaeus, 1766), and S. intermedius (Agassiz, 1829)

Western Atlantic synodontid species were studied as part of an ongoing effort to reanalyze Caribbean shorefish diversity. A neighbor-joining tree constructed from cytochrome c oxidase I (COI) data revealed 2 highly divergent genetic lineages within both Synodus intermedius (Agassiz, 1829) (Sand Diver) and S. foetens (Linnaeus, 1766) (Inshore Lizardfish). A new species, Synodus macrostigmus, is described for one of the S. intermedius lineages. Synodus macrostigmus and S. intermedius differ in number of lateral-line scales, caudal pigmentation, size of the scapular blotch, and shape of the anterior-nostril flap. Synodus macrostigmus and S. intermedius have overlapping geographic and depth distributions, but S. macrostigmus generally inhabits deeper water (>28 m) than does S. intermedius and is known only from coastal waters of the southeastern United States and the Gulf of Mexico, in contrast to those areas and the Caribbean for S. intermedius. Synodus bondi Fowler, 1939, is resurrected from the synonymy of S. foetens for one of the S. foetens genetic lineages. The 2 species differ in length and shape of the snout, number of anal-fin rays, and shape of the anterior-nostril flap. Synodus bondi and S. foetens co-occur in the central Caribbean, but S. bondi otherwise has a more southerly distribution than does S. foetens. Redescriptions are provided for S. intermedius, S. foetens, and S. bondi. Neotypes are designated for S. intermedius and S. foetens. A revised key to Synodus species in the western Atlantic is presented

Complete Genomic Sequence of Bacteriophage Felix O1†

Bacteriophage O1 is a Myoviridae A1 group member used historically for identifying Salmonella. Sequencing revealed a single, linear, 86,155-base-pair genome with 39% average G+C content, 131 open reading frames, and 22 tRNAs. Closest protein homologs occur in Erwinia amylovora phage φEa21-4 and Escherichia coli phage wV8. Proteomic analysis indentified structural proteins: Gp23, Gp36 (major tail protein), Gp49, Gp53, Gp54, Gp55, Gp57, Gp58 (major capsid protein), Gp59, Gp63, Gp64, Gp67, Gp68, Gp69, Gp73, Gp74 and Gp77 (tail fiber). Based on phage-host codon differences, 7 tRNAs could affect translation rate during infection. Introns, holin-lysin cassettes, bacterial toxin homologs and host RNA polymerase-modifying genes were absent

Direct-coupling analysis of residue co-evolution captures native contacts across many protein families

The similarity in the three-dimensional structures of homologous proteins imposes strong constraints on their sequence variability. It has long been suggested that the resulting correlations among amino acid compositions at different sequence positions can be exploited to infer spatial contacts within the tertiary protein structure. Crucial to this inference is the ability to disentangle direct and indirect correlations, as accomplished by the recently introduced Direct Coupling Analysis (DCA) (Weigt et al. (2009) Proc Natl Acad Sci 106:67). Here we develop a computationally efficient implementation of DCA, which allows us to evaluate the accuracy of contact prediction by DCA for a large number of protein domains, based purely on sequence information. DCA is shown to yield a large number of correctly predicted contacts, recapitulating the global structure of the contact map for the majority of the protein domains examined. Furthermore, our analysis captures clear signals beyond intra- domain residue contacts, arising, e.g., from alternative protein conformations, ligand- mediated residue couplings, and inter-domain interactions in protein oligomers. Our findings suggest that contacts predicted by DCA can be used as a reliable guide to facilitate computational predictions of alternative protein conformations, protein complex formation, and even the de novo prediction of protein domain structures, provided the existence of a large number of homologous sequences which are being rapidly made available due to advances in genome sequencing.Comment: 28 pages, 7 figures, to appear in PNA

Integrated Analysis of Residue Coevolution and Protein Structure in ABC Transporters

Intraprotein side chain contacts can couple the evolutionary process of amino acid substitution at one position to that at another. This coupling, known as residue coevolution, may vary in strength. Conserved contacts thus not only define 3-dimensional protein structure, but also indicate which residue-residue interactions are crucial to a protein’s function. Therefore, prediction of strongly coevolving residue-pairs helps clarify molecular mechanisms underlying function. Previously, various coevolution detectors have been employed separately to predict these pairs purely from multiple sequence alignments, while disregarding available structural information. This study introduces an integrative framework that improves the accuracy of such predictions, relative to previous approaches, by combining multiple coevolution detectors and incorporating structural contact information. This framework is applied to the ABC-B and ABC-C transporter families, which include the drug exporter P-glycoprotein involved in multidrug resistance of cancer cells, as well as the CFTR chloride channel linked to cystic fibrosis disease. The predicted coevolving pairs are further analyzed based on conformational changes inferred from outward- and inward-facing transporter structures. The analysis suggests that some pairs coevolved to directly regulate conformational changes of the alternating-access transport mechanism, while others to stabilize rigid-body-like components of the protein structure. Moreover, some identified pairs correspond to residues previously implicated in cystic fibrosis

Structural and Functional Characterization of Mature Forms of Metalloprotease E495 from Arctic Sea-Ice Bacterium Pseudoalteromonas sp. SM495

E495 is the most abundant protease secreted by the Arctic sea-ice bacterium Pseudoalteromonas sp. SM495. As a thermolysin family metalloprotease, E495 was found to have multiple active forms in the culture of strain SM495. E495-M (containing only the catalytic domain) and E495-M-C1 (containing the catalytic domain and one PPC domain) were two stable mature forms, and E495-M-C1-C2 (containing the catalytic domain and two PPC domains) might be an intermediate. Compared to E495-M, E495-M-C1 had similar affinity and catalytic efficiency to oligopeptides, but higher affinity and catalytic efficiency to proteins. The PPC domains from E495 were expressed as GST-fused proteins. Both of the recombinant PPC domains were shown to have binding ability to proteins C-phycocyanin and casein, and domain PPC1 had higher affinity to C-phycocyanin than domain PPC2. These results indicated that the domain PPC1 in E495-M-C1 could be helpful in binding protein substrate, and therefore, improving the catalytic efficiency. Site-directed mutagenesis on the PPC domains showed that the conserved polar and aromatic residues, D26, D28, Y30, Y/W65, in the PPC domains played key roles in protein binding. Our study may shed light on the mechanism of organic nitrogen degradation in the Arctic sea ice

Helicobacter pylori infection and circulating ghrelin levels - A systematic review

BACKGROUND: The nature of the association between ghrelin, an orexigenic hormone produced mainly in the stomach, and Helicobacter pylori (H pylori), a bacterium that colonises the stomach, is still controversial. We examined available evidence to determine whether an association exists between the two; and if one exists, in what direction. METHODS: We reviewed original English language studies on humans reporting circulating ghrelin levels in H pylori infected and un-infected participants; and circulating ghrelin levels before and after H pylori eradication. Meta-analyses were conducted for eligible studies by combining study specific estimates using the inverse variance method with weighted average for continuous outcomes in a random effects model. RESULTS: Seventeen out of 27 papers that reported ghrelin levels in H pylori positive and negative subjects found lower circulating ghrelin levels in H pylori positive subjects; while 10 found no difference. A meta-analysis of 19 studies with a total of 1801 participants showed a significantly higher circulating ghrelin concentration in H pylori negative participants than in H pylori positive participants (Effect estimate (95%CI) = -0.48 (-0.60, -0.36)). However, eradicating H pylori did not have any significant effect on circulating ghrelin levels (Effect estimate (95% CI) = 0.08 (-0.33, 0.16); Test for overall effect: Z = 0.67 (P = 0.5)). CONCLUSIONS: We conclude that circulating ghrelin levels are lower in H pylori infected people compared to those not infected; but the relationship between circulating ghrelin and eradication of H pylori is more complex

Structural Constraints Identified with Covariation Analysis in Ribosomal RNA

Covariation analysis is used to identify those positions with similar patterns of sequence variation in an alignment of RNA sequences. These constraints on the evolution of two positions are usually associated with a base pair in a helix. While mutual information (MI) has been used to accurately predict an RNA secondary structure and a few of its tertiary interactions, early studies revealed that phylogenetic event counting methods are more sensitive and provide extra confidence in the prediction of base pairs. We developed a novel and powerful phylogenetic events counting method (PEC) for quantifying positional covariation with the Gutell lab’s new RNA Comparative Analysis Database (rCAD). The PEC and MI-based methods each identify unique base pairs, and jointly identify many other base pairs. In total, both methods in combination with an N-best and helix-extension strategy identify the maximal number of base pairs. While covariation methods have effectively and accurately predicted RNAs secondary structure, only a few tertiary structure base pairs have been identified. Analysis presented herein and at the Gutell lab’s Comparative RNA Web (CRW) Site reveal that the majority of these latter base pairs do not covary with one another. However, covariation analysis does reveal a weaker although significant covariation between sets of nucleotides that are in proximity in the three-dimensional RNA structure. This reveals that covariation analysis identifies other types of structural constraints beyond the two nucleotides that form a base pair

A review of bronchiolitis obliterans syndrome and therapeutic strategies

Lung transplantation is an important treatment option for patients with advanced lung disease. Survival rates for lung transplant recipients have improved; however, the major obstacle limiting better survival is bronchiolitis obliterans syndrome (BOS). In the last decade, survival after lung retransplantation has improved for transplant recipients with BOS. This manuscript reviews BOS along with the current therapeutic strategies, including recent outcomes for lung retransplantation