Remote homology detection of integral membrane proteins using conserved sequence features

A Continuum of Evolving De Novo Genes Drives Protein-Coding Novelty in Drosophila

Orphan genes, lacking detectable homologs in outgroup species, typically represent 10-30% of eukaryotic genomes. Efforts to find the source of these young genes indicate that de novo emergence from non-coding DNA may in part explain their prevalence. Here, we investigate the roots of orphan gene emergence in the Drosophila genus. Across the annotated proteomes of twelve species, we find 6297 orphan genes within 4953 taxon-specific clusters of orthologs. By inferring the ancestral DNA as non-coding for between 550 and 2467 (8.7-39.2%) of these genes, we describe for the first time how de novo emergence contributes to the abundance of clade-specific Drosophila genes. In support of them having functional roles, we show that de novo genes have robust expression and translational support. However, the distinct nucleotide sequences of de novo genes, which have characteristics intermediate between intergenic regions and conserved genes, reflect their recent birth from non-coding DNA. We find that de novo genes encode more disordered proteins than both older genes and intergenic regions. Together, our results suggest that gene emergence from non-coding DNA provides an abundant source of material for the evolution of new proteins. Following gene birth, gradual evolution over large evolutionary timescales moulds sequence properties towards those of conserved genes, resulting in a continuum of properties whose starting points depend on the nucleotide sequences of an initial pool of novel genes.

Diversity of the Epsilonproteobacteria Dsb (disulfide bond) systems

The bacterial proteins of the Dsb family-important components of the post-translational protein modification system-catalyze the formation of disulfide bridges, a process that is crucial for protein structure stabilization and activity. Dsb systems play an essential role in the assembly of many virulence factors. Recent rapid advances in global analysis of bacteria have thrown light on the enormous diversity among bacterial Dsb systems. While the Escherichia coli disulfide bond-forming system is quite well understood, the mechanisms of action of Dsb systems in other bacteria, including members of class Epsilonproteobacteria that contain pathogenic and non-pathogenic bacteria colonizing extremely diverse ecological niches, are poorly characterized. Here we present a review of current knowledge on Epsilonproteobacteria Dsb systems. We have focused on the Dsb systems of Campylobacter spp. and Helicobacter spp. because our knowledge about Dsb proteins of Wolinella and Arcobacter spp. is still scarce and comes mainly from bioinformatic studies. Helicobacter pylori is a common human pathogen that colonizes the gastric epithelium of humans with severe consequences. Campylobacter spp. is a leading cause of zoonotic enteric bacterial infections in most developed and developing nations. We focus on various aspects of the diversity of the Dsb systems and their influence on pathogenicity, particularly because Dsb proteins are considered as potential targets for a new class of anti-virulence drugs to treat human infections by Campylobacter or Helicobacter spp.

Alignment of helical membrane protein sequences using AlignMe

Few sequence alignment methods have been designed specifically for integral membrane proteins, even though these important proteins have distinct evolutionary and structural properties that might affect their alignments. Existing approaches typically consider membrane-related information either by using membrane-specific substitution matrices or by assigning distinct penalties for gap creation in transmembrane and non-transmembrane regions. Here, we ask whether favoring matching of predicted transmembrane segments within a standard dynamic programming algorithm can improve the accuracy of pairwise membrane protein sequence alignments. We tested various strategies using a specifically designed program called AlignMe. An updated set of homologous membrane protein structures, called HOMEP2, was used as a reference for optimizing the gap penalties. The best of the membrane-protein optimized approaches were then tested on an independent reference set of membrane protein sequence alignments from the BAliBASE collection. When secondary structure (S) matching was combined with evolutionary information (using a position-specific substitution matrix (P)), in an approach we called AlignMePS, the resultant pairwise alignments were typically among the most accurate over a broad range of sequence similarities when compared to available methods. Matching transmembrane predictions (T), in addition to evolutionary information, and secondary-structure predictions, in an approach called AlignMePST, generally reduces the accuracy of the alignments of closely-related proteins in the BAliBASE set relative to AlignMePS, but may be useful in cases of extremely distantly related proteins for which sequence information is less informative. The open source AlignMe code is available at https://sourceforge.net/projects/alignme/, and at http://www.forrestlab.org, along with an online server and the HOMEP2 data set.

Internal duplications in α-helical membrane protein topologies are common but the nonduplicated forms are rare

Many α-helical membrane proteins contain internal symmetries, indicating that they might have evolved through a gene duplication and fusion event. Here, we have characterized internal duplications among membrane proteins of known structure and in three complete genomes. We found that the majority of large transmembrane (TM) proteins contain an internal duplication. The duplications found showed a large variability both in the number of TM-segments included and in their orientation. Surprisingly, an approximately equal number of antiparallel duplications and parallel duplications were found. However, of all 11 superfamilies with an internal duplication, only for one, the AcrB Multidrug Efflux Pump, the duplicated unit could be found in its nonduplicated form. An evolutionary analysis of the AcrB homologs indicates that several independent fusions have occurred, including the fusion of the SecD and SecF proteins into the 12-TM-protein SecDF in Brucella and Staphylococcus aureus. In one additional case, the Vitamin B12 transporter-like ABC transporters, the protein had undergone an additional fusion to form protein with 20 TM-helices in several bacterial genomes. Finally, homologs to all human membrane proteins were used to detect the presence of duplicated and nonduplicated proteins. This confirmed that only in rare cases can homologs with different duplication status be found, although internal symmetry is frequent among these proteins. One possible explanation is that it is frequent that duplication and fusion events happen simultaneously and that there is almost always a strong selective advantage for the fused form.

On the mechanism of chloroquine resistance in Plasmodium falciparum

Resistance to chloroquine of malaria strains is known to be associated with a parasite protein named PfCRT, the mutated form of which is able to reduce chloroquine accumulation in the digestive vacuole of the pathogen. Whether the protein mediates extrusion of the drug acting as a channel or as a carrier and which is the protonation state of its chloroquine substrate is the subject of a scientific debate. We present here an analytical approach that explores which combination of hypotheses on the mechanism of transport and the protonation state of chloroquine are consistent with available equilibrium experimental data. We show that the available experimental data are not, by themselves, sufficient to conclude whether the protein acts as a channel or as a transporter, which explains the origin of their different interpretation by different authors. Interestingly, though, each of the two models is only consistent with a subset of hypotheses on the protonation state of the transported molecule. The combination of these results with a sequence and structure analysis of PfCRT, which strongly suggests that the molecule is a carrier, indicates that the transported species is either or both the mono and di-protonated forms of chloroquine. We believe that our results, besides shedding light on the mechanism of chloroquine resistance in P. falciparum, have implications for the development of novel therapies against resistant malaria strains and demonstrate the usefulness of an approach combining systems biology strategies with structural bioinformatics and experimental data.

More than 1,001 problems with protein domain databases: transmembrane regions, signal peptides and the issue of sequence homology

Large-scale genome sequencing gained general importance for life science because functional annotation of otherwise experimentally uncharacterized sequences is made possible by the theory of biomolecular sequence homology. Historically, the paradigm of similarity of protein sequences implying common structure, function and ancestry was generalized based on studies of globular domains. Having the same fold imposes strict conditions over the packing in the hydrophobic core requiring similarity of hydrophobic patterns. The implications of sequence similarity among non-globular protein segments have not been studied to the same extent; nevertheless, homology considerations are silently extended for them. This appears especially detrimental in the case of transmembrane helices (TMs) and signal peptides (SPs) where sequence similarity is necessarily a consequence of physical requirements rather than common ancestry. Thus, matching of SPs/TMs creates the illusion of matching hydrophobic cores. Therefore, inclusion of SPs/TMs into domain models can give rise to wrong annotations. More than 1001 domains among the 10,340 models of Pfam release 23 and 18 domains of SMART version 6 (out of 809) contain SP/TM regions. As expected, fragment-mode HMM searches generate promiscuous hits limited to solely the SP/TM part among clearly unrelated proteins. More worryingly, we show explicit examples that the scores of clearly false-positive hits, even in global-mode searches, can be elevated into the significance range just by matching the hydrophobic runs. In the PIR iProClass database v3.74 using conservative criteria, we find that at least between 2.1% and 13.6% of its annotated Pfam hits appear unjustified for a set of validated domain models. Thus, false-positive domain hits enforced by SP/TM regions can lead to dramatic annotation errors where the hit has nothing in common with the problematic domain model except the SP/TM region itself. We suggest a workflow of flagging problematic hits arising from SP/TM-containing models for critical reconsideration by annotation users.

Sequence homology is a fundamental principle of biology. It implies common phylogenetic ancestry of genes and, subsequently, similarity of their protein products with regard to amino acid sequence, three-dimensional structure and molecular and cellular function. Originally an esoteric concept, homology with the proxy of sequence similarity is used to justify the transfer of functional annotation from well-studied protein examples to new sequences. Yet, functional annotation via sequence similarity seems to have hit a plateau in recent years since relentless annotation transfer led to error propagation across sequence databases; thus, leading experimental follow-up work astray. It must be emphasized that the trinity of sequence, 3D structural and functional similarity has only been proven for globular segments of proteins. For non-globular regions, similarity of sequence is not necessarily a result of divergent evolution from a common ancestor but the consequence of amino acid sequence bias. In our investigation, we found that protein domain databases contain many domain models with transmembrane regions and signal peptides, non-globular segments of proteins having hydrophobic bias. Many proteins have inherited completely wrong function assignments from these domain models. We fear that future function predictions will turn out futile if this issue is not immediately addressed.