Phagosomal transporters (Phts), required for intracellular growth of Legionella pneumophila, comprise a novel family of multispanning a-helical proteins within the major facilitator superfamily (MFS). The members of this family derive exclusively from bacteria. Multiple paralogues are present in a restricted group of Alpha-and Gammaproteobacteria, but single members were also found in Chlamydia and Cyanobacteria. Their protein sequences were aligned, yielding a phylogenetic tree showing the relations of the proteins to each other. Topological analyses revealed a probable 12 a-helical transmembrane segment (TMS) topology. Motif identification and statistical analyses provided convincing evidence that these proteins arose from a six TMS precursor by intragenic duplication. The phylogenetic tree revealed some potential orthologous relationships, suggestive of common function. However, several probable examples of lateral transfer of the encoding genetic material between bacteria were identified and analysed. The Pht family most closely resembles a smaller MFS family (the UMF9 family) with no functionally characterized members. However, the UMF9 family occurs in a broader range of prokaryotic organism types, including Archaea. These two families differ in that organisms bearing members of the Pht family often have numerous paralogues, whereas organisms bearing members of the UMF9 family never have more than two. This work serves to characterize two novel families within the MFS and provides compelling evidence for horizontal transfer of some of the family members.

Directed evolution is a technique that enables the identification of mutants of a particular protein that carry a desired property by successive rounds of random mutagenesis, screening, and selection. This technique has many applications, including the development of G proteincoupled receptor-based biosensors and designer drugs for personalized medicine. Although effective, directed evolution is not without challenges and can greatly benefit from the development of computational techniques to predict the functional outcome of single-point amino acid substitutions. In this article, we describe a molecular dynamics-based approach to predict the effects of single amino acid substitutions on agonist binding (salicin) to a human bitter taste receptor (hT2R16). An experimentally determined functional map of single-point amino acid substitutions was used to validate the whole-protein molecular dynamics-based predictive functions. Molecular docking was used to construct a wild-type agonist-receptor complex, providing a starting structure for single-point substitution simulations. The effects of each single amino acid substitution in the functional response of the receptor to its agonist were estimated using three binding energy schemes with increasing inclusion of solvation effects. We show that molecular docking combined with molecular mechanics simulations of single-point mutants of the agonist-receptor complex accurately predicts the functional outcome of single amino acid substitutions in a human bitter taste receptor.