Integration of Evolutionary Features for the Identification of Functionally Important Residues in Major Facilitator Superfamily Transporters

The identification of functionally important residues is an important challenge for understanding the molecular mechanisms of proteins. Membrane protein transporters operate two-state allosteric conformational changes using functionally important cooperative residues that mediate long-range communication from the substrate binding site to the translocation pathway. In this study, we identified functionally important cooperative residues of membrane protein transporters by integrating sequence conservation and co-evolutionary information. A newly derived evolutionary feature, the co-evolutionary coupling number, was introduced to measure the connectivity of co-evolving residue pairs and was integrated with the sequence conservation score. We tested this method on three Major Facilitator Superfamily (MFS) transporters, LacY, GlpT, and EmrD. MFS transporters are an important family of membrane protein transporters, which utilize diverse substrates, catalyze different modes of transport using unique combinations of functional residues, and have enough characterized functional residues to validate the performance of our method. We found that the conserved cores of evolutionarily coupled residues are involved in specific substrate recognition and translocation of MFS transporters. Furthermore, a subset of the residues forms an interaction network connecting functional sites in the protein structure. We also confirmed that our method is effective on other membrane protein transporters. Our results provide insight into the location of functional residues important for the molecular mechanisms of membrane protein transporters

Two polymorphisms facilitate differences in plasticity between two chicken major histocompatibility complex class I proteins

Major histocompatibility complex class I molecules (MHC I) present peptides to cytotoxic T-cells at the surface of almost all nucleated cells. The function of MHC I molecules is to select high affinity peptides from a large intracellular pool and they are assisted in this process by co-factor molecules, notably tapasin. In contrast to mammals, MHC homozygous chickens express a single MHC I gene locus, termed BF2, which is hypothesised to have co-evolved with the highly polymorphic tapasin within stable haplotypes. The BF2 molecules of the B15 and B19 haplotypes have recently been shown to differ in their interactions with tapasin and in their peptide selection properties. This study investigated whether these observations might be explained by differences in the protein plasticity that is encoded into the MHC I structure by primary sequence polymorphisms. Furthermore, we aimed to demonstrate the utility of a complimentary modelling approach to the understanding of complex experimental data. Combining mechanistic molecular dynamics simulations and the primary sequence based technique of statistical coupling analysis, we show how two of the eight polymorphisms between BF2*15:01 and BF2*19:01 facilitate differences in plasticity. We show that BF2*15:01 is intrinsically more plastic than BF2*19:01, exploring more conformations in the absence of peptide. We identify a protein sector of contiguous residues connecting the membrane bound ?3 domain and the heavy chain peptide binding site. This sector contains two of the eight polymorphic residues. One is residue 22 in the peptide binding domain and the other 220 is in the ?3 domain, a putative tapasin binding site. These observations are in correspondence with the experimentally observed functional differences of these molecules and suggest a mechanism for how modulation of MHC I plasticity by tapasin catalyses peptide selection allosterically

The tomato Prf complex is a molecular trap for bacterial effectors based on Pto transphosphorylation

The bacteria Pseudomonas syringae is a pathogen of many crop species and one of the model pathogens for studying plant and bacterial arms race coevolution. In the current model, plants perceive bacteria pathogens via plasma membrane receptors, and recognition leads to the activation of general defenses. In turn, bacteria inject proteins called effectors into the plant cell to prevent the activation of immune responses. AvrPto and AvrPtoB are two such proteins that inhibit multiple plant kinases. The tomato plant has reacted to these effectors by the evolution of a cytoplasmic resistance complex. This complex is compromised of two proteins, Prf and Pto kinase, and is capable of recognizing the effector proteins. How the Pto kinase is able to avoid inhibition by the effector proteins is currently unknown. Our data shows how the tomato plant utilizes dimerization of resistance proteins to gain advantage over the faster evolving bacterial pathogen. Here we illustrate that oligomerisation of Prf brings into proximity two Pto kinases allowing them to avoid inhibition by the effectors by transphosphorylation and to activate immune responses

Inferring stabilizing mutations from protein phylogenies : application to influenza hemagglutinin

One selection pressure shaping sequence evolution is the requirement that a protein fold with sufficient stability to perform its biological functions. We present a conceptual framework that explains how this requirement causes the probability that a particular amino acid mutation is fixed during evolution to depend on its effect on protein stability. We mathematically formalize this framework to develop a Bayesian approach for inferring the stability effects of individual mutations from homologous protein sequences of known phylogeny. This approach is able to predict published experimentally measured mutational stability effects (ΔΔG values) with an accuracy that exceeds both a state-of-the-art physicochemical modeling program and the sequence-based consensus approach. As a further test, we use our phylogenetic inference approach to predict stabilizing mutations to influenza hemagglutinin. We introduce these mutations into a temperature-sensitive influenza virus with a defect in its hemagglutinin gene and experimentally demonstrate that some of the mutations allow the virus to grow at higher temperatures. Our work therefore describes a powerful new approach for predicting stabilizing mutations that can be successfully applied even to large, complex proteins such as hemagglutinin. This approach also makes a mathematical link between phylogenetics and experimentally measurable protein properties, potentially paving the way for more accurate analyses of molecular evolution

Evolution and diversity of secretome genes in the apicomplexan parasite Theileria annulata

<b>BACKGROUND</b>: Little is known about how apicomplexan parasites have evolved to infect different host species and cell types. Theileria annulata and Theileria parva invade and transform bovine leukocytes but each species favours a different host cell lineage. Parasite-encoded proteins secreted from the intracellular macroschizont stage within the leukocyte represent a critical interface between host and pathogen systems. Genome sequencing has revealed that several Theileria-specific gene families encoding secreted proteins are positively selected at the inter-species level, indicating diversification between the species. We extend this analysis to the intra-species level, focusing on allelic diversity of two major secretome families. These families represent a well-characterised group of genes implicated in control of the host cell phenotype and a gene family of unknown function. To gain further insight into their evolution and function, this study investigates whether representative genes of these two families are diversifying or constrained within the T. annulata population. <b>RESULTS</b>: Strong evidence is provided that the sub-telomerically encoded SVSP family and the host-nucleus targeted TashAT family have evolved under contrasting pressures within natural T. annulata populations. SVSP genes were found to possess atypical codon usage and be evolving neutrally, with high levels of nucleotide substitutions and multiple indels. No evidence of geographical sub-structuring of allelic sequences was found. In contrast, TashAT family genes, implicated in control of host cell gene expression, are strongly conserved at the protein level and geographically sub-structured allelic sequences were identified among Tunisian and Turkish isolates. Although different copy numbers of DNA binding motifs were identified in alleles of TashAT proteins, motif periodicity was strongly maintained, implying conserved functional activity of these sites. <b>CONCLUSIONS</b>: This analysis provides evidence that two distinct secretome genes families have evolved under contrasting selective pressures. The data supports current hypotheses regarding the biological role of TashAT family proteins in the management of host cell phenotype that may have evolved to allow adaptation of T. annulata to a specific host cell lineage. We provide new evidence of extensive allelic diversity in representative members of the enigmatic SVSP gene family, which supports a putative role for the encoded products in subversion of the host immune response

Structural insights into Clostridium perfringens delta toxin pore formation

Clostridium perfringens Delta toxin is one of the three hemolysin-like proteins produced by C. perfringens type C and possibly type B strains. One of the others, NetB, has been shown to be the major cause of Avian Nectrotic Enteritis, which following the reduction in use of antibiotics as growth promoters, has become an emerging disease of industrial poultry. Delta toxin itself is cytotoxic to the wide range of human and animal macrophages and platelets that present GM2 ganglioside on their membranes. It has sequence similarity with Staphylococcus aureus β-pore forming toxins and is expected to heptamerize and form pores in the lipid bilayer of host cell membranes. Nevertheless, its exact mode of action remains undetermined. Here we report the 2.4 Å crystal structure of monomeric Delta toxin. The superposition of this structure with the structure of the phospholipid-bound F component of S. aureus leucocidin (LukF) revealed that the glycerol molecules bound to Delta toxin and the phospholipids in LukF are accommodated in the same hydrophobic clefts, corresponding to where the toxin is expected to latch onto the membrane, though the binding sites show significant differences. From structure-based sequence alignment with the known structure of staphylococcal α-hemolysin, a model of the Delta toxin pore form has been built. Using electron microscopy, we have validated our model and characterized the Delta toxin pore on liposomes. These results highlight both similarities and differences in the mechanism of Delta toxin (and by extension NetB) cytotoxicity from that of the staphylococcal pore-forming toxins

The hydrogen peroxide-sensitive proteome of the chloroplast in vitro and in vivo

Muthuramalingam M, Matros A, Scheibe R, Mock H-P, Dietz K-J. The hydrogen peroxide-sensitive proteome of the chloroplast in vitro and in vivo. Frontiers in Plant Science. 2013;4:54-1-54-14.Hydrogen peroxide (H2O2) evolves during cellular metabolism and accumulates under various stresses causing serious redox imbalances. Many proteomics studies aiming to identify proteins sensitive to H2O2 used concentrations that were above the physiological range. Here the chloroplast proteins were subjected to partial oxidation by exogenous addition of H2O2 equivalent to 10% of available protein thiols which allowed for the identification of the primary targets of oxidation. The chosen redox proteomic approach employed differential labeling of non-oxidized and oxidized thiols using sequential alkylation with N-ethylmaleimide and biotin maleimide. The in vitro identified proteins are involved in carbohydrate metabolism, photosynthesis, redox homeostasis, and nitrogen assimilation. By using methyl viologen that induces oxidative stress in vivo, mostly the same primary targets of oxidation were identified and several oxidation sites were annotated. Ribulose-1,5-bisphosphate (RubisCO) was a primary oxidation target. Due to its high abundance, RubisCO is suggested to act as a chloroplast redox buffer to maintain a suitable redox state, even in the presence of increased reactive oxygen species release. 2-cysteine peroxiredoxins (2-Cys Prx) undergo redox-dependent modifications and play important roles in antioxidant defense and signaling. The identification of 2-Cys Prx was expected based on its high affinity to H2O2 and is considered as a proof of concept for the approach. Targets of Trx, such as phosphoribulokinase, glyceraldehyde-3-phosphate dehydrogenase, transketolase, and sedoheptulose-1,7-bisphosphatase have at least one regulatory disulfide bridge which supports the conclusion that the identified proteins undergo reversible thiol oxidation. In conclusion, the presented approach enabled the identification of early targets of H2O2 oxidation within the cellular proteome under physiological experimental conditions