Evolution of a Signaling Nexus Constrained by Protein Interfaces and Conformational States

Heterotrimeric G proteins act as the physical nexus between numerous receptors that respond to extracellular signals and proteins that drive the cytoplasmic response. The Gα subunit of the G protein, in particular, is highly constrained due to its many interactions with proteins that control or react to its conformational state. Various organisms contain differing sets of Gα-interacting proteins, clearly indicating that shifts in sequence and associated Gα functionality were acquired over time. These numerous interactions constrained much of Gα evolution; yet Gα has diversified, through poorly understood processes, into several functionally specialized classes, each with a unique set of interacting proteins. Applying a synthetic sequence-based approach to mammalian Gα subunits, we established a set of seventy-five evolutionarily important class-distinctive residues, sites where a single Gα class is differentiated from the three other classes. We tested the hypothesis that shifts at these sites are important for class-specific functionality. Importantly, we mapped known and well-studied class-specific functionalities from all four mammalian classes to sixteen of our class-distinctive sites, validating the hypothesis. Our results show how unique functionality can evolve through the recruitment of residues that were ancestrally functional. We also studied acquisition of functionalities by following these evolutionarily important sites in non-mammalian organisms. Our results suggest that many class-distinctive sites were established early on in eukaryotic diversification and were critical for the establishment of new Gα classes, whereas others arose in punctuated bursts throughout metazoan evolution. These Gα class-distinctive residues are rational targets for future structural and functional studies

CodonTest: Modeling Amino Acid Substitution Preferences in Coding Sequences

Codon models of evolution have facilitated the interpretation of selective forces operating on genomes. These models, however, assume a single rate of non-synonymous substitution irrespective of the nature of amino acids being exchanged. Recent developments have shown that models which allow for amino acid pairs to have independent rates of substitution offer improved fit over single rate models. However, these approaches have been limited by the necessity for large alignments in their estimation. An alternative approach is to assume that substitution rates between amino acid pairs can be subdivided into rate classes, dependent on the information content of the alignment. However, given the combinatorially large number of such models, an efficient model search strategy is needed. Here we develop a Genetic Algorithm (GA) method for the estimation of such models. A GA is used to assign amino acid substitution pairs to a series of rate classes, where is estimated from the alignment. Other parameters of the phylogenetic Markov model, including substitution rates, character frequencies and branch lengths are estimated using standard maximum likelihood optimization procedures. We apply the GA to empirical alignments and show improved model fit over existing models of codon evolution. Our results suggest that current models are poor approximations of protein evolution and thus gene and organism specific multi-rate models that incorporate amino acid substitution biases are preferred. We further anticipate that the clustering of amino acid substitution rates into classes will be biologically informative, such that genes with similar functions exhibit similar clustering, and hence this clustering will be useful for the evolutionary fingerprinting of genes

Modeling protein network evolution under genome duplication and domain shuffling

<p>Abstract</p> <p>Background</p> <p>Successive whole genome duplications have recently been firmly established in all major eukaryote kingdoms. Such <it>exponential </it>evolutionary processes must have largely contributed to shape the topology of protein-protein interaction (PPI) networks by outweighing, in particular, all <it>time-linear </it>network growths modeled so far.</p> <p>Results</p> <p>We propose and solve a mathematical model of PPI network evolution under successive genome duplications. This demonstrates, from first principles, that evolutionary conservation and scale-free topology are intrinsically linked properties of PPI networks and emerge from <it>i) </it>prevailing <it>exponential </it>network dynamics under duplication and <it>ii) asymmetric divergence </it>of gene duplicates. While required, we argue that this asymmetric divergence arises, in fact, spontaneously at the level of protein-binding sites. This supports a refined model of PPI network evolution in terms of protein domains under exponential and asymmetric duplication/divergence dynamics, with multidomain proteins underlying the combinatorial formation of protein complexes. Genome duplication then provides a powerful source of PPI network innovation by promoting local rearrangements of multidomain proteins on a genome wide scale. Yet, we show that the overall conservation and topology of PPI networks are robust to extensive domain shuffling of multidomain proteins as well as to finer details of protein interaction and evolution. Finally, large scale features of <it>direct </it>and <it>indirect </it>PPI networks of <it>S. cerevisiae </it>are well reproduced numerically with only two adjusted parameters of clear biological significance (<it>i.e</it>. network effective growth rate and average number of protein-binding domains per protein).</p> <p>Conclusion</p> <p>This study demonstrates the statistical consequences of genome duplication and domain shuffling on the conservation and topology of PPI networks over a broad evolutionary scale across eukaryote kingdoms. In particular, scale-free topologies of PPI networks, which are found to be robust to extensive shuffling of protein domains, appear to be a simple consequence of the conservation of protein-binding domains under asymmetric duplication/divergence dynamics in the course of evolution.</p

Integrated Assessment of Genomic Correlates of Protein Evolutionary Rate

Rates of evolution differ widely among proteins, but the causes and consequences of such differences remain under debate. With the advent of high-throughput functional genomics, it is now possible to rigorously assess the genomic correlates of protein evolutionary rate. However, dissecting the correlations among evolutionary rate and these genomic features remains a major challenge. Here, we use an integrated probabilistic modeling approach to study genomic correlates of protein evolutionary rate in Saccharomyces cerevisiae. We measure and rank degrees of association between (i) an approximate measure of protein evolutionary rate with high genome coverage, and (ii) a diverse list of protein properties (sequence, structural, functional, network, and phenotypic). We observe, among many statistically significant correlations, that slowly evolving proteins tend to be regulated by more transcription factors, deficient in predicted structural disorder, involved in characteristic biological functions (such as translation), biased in amino acid composition, and are generally more abundant, more essential, and enriched for interaction partners. Many of these results are in agreement with recent studies. In addition, we assess information contribution of different subsets of these protein properties in the task of predicting slowly evolving proteins. We employ a logistic regression model on binned data that is able to account for intercorrelation, non-linearity, and heterogeneity within features. Our model considers features both individually and in natural ensembles (“meta-features”) in order to assess joint information contribution and degree of contribution independence. Meta-features based on protein abundance and amino acid composition make strong, partially independent contributions to the task of predicting slowly evolving proteins; other meta-features make additional minor contributions. The combination of all meta-features yields predictions comparable to those based on paired species comparisons, and approaching the predictive limit of optimal lineage-insensitive features. Our integrated assessment framework can be readily extended to other correlational analyses at the genome scale

Attachment and Entry of Chlamydia Have Distinct Requirements for Host Protein Disulfide Isomerase

Chlamydia is an obligate intracellular pathogen that causes a wide range of diseases in humans. Attachment and entry are key processes in infectivity and subsequent pathogenesis of Chlamydia, yet the mechanisms governing these interactions are unknown. It was recently shown that a cell line, CHO6, that is resistant to attachment, and thus infectivity, of multiple Chlamydia species has a defect in protein disulfide isomerase (PDI) N–terminal signal sequence processing. Ectopic expression of PDI in CHO6 cells led to restoration of Chlamydia attachment and infectivity; however, the mechanism leading to this recovery was not ascertained. To advance our understanding of the role of PDI in Chlamydia infection, we used RNA interference to establish that cellular PDI is essential for bacterial attachment to cells, making PDI the only host protein identified as necessary for attachment of multiple species of Chlamydia. Genetic complementation and PDI-specific inhibitors were used to determine that cell surface PDI enzymatic activity is required for bacterial entry into cells, but enzymatic function was not required for bacterial attachment. We further determined that it is a PDI-mediated reduction at the cell surface that triggers bacterial uptake. While PDI is necessary for Chlamydia attachment to cells, the bacteria do not appear to utilize plasma membrane–associated PDI as a receptor, suggesting that Chlamydia binds a cell surface protein that requires structural association with PDI. Our findings demonstrate that PDI has two essential and independent roles in the process of chlamydial infectivity: it is structurally required for chlamydial attachment, and the thiol-mediated oxido-reductive function of PDI is necessary for entry

Stability of domain structures in multi-domain proteins

Multi-domain proteins have many advantages with respect to stability and folding inside cells. Here we attempt to understand the intricate relationship between the domain-domain interactions and the stability of domains in isolation. We provide quantitative treatment and proof for prevailing intuitive ideas on the strategies employed by nature to stabilize otherwise unstable domains. We find that domains incapable of independent stability are stabilized by favourable interactions with tethered domains in the multi-domain context. Stability of such folds to exist independently is optimized by evolution. Specific residue mutations in the sites equivalent to inter-domain interface enhance the overall solvation, thereby stabilizing these domain folds independently. A few naturally occurring variants at these sites alter communication between domains and affect stability leading to disease manifestation. Our analysis provides safe guidelines for mutagenesis which have attractive applications in obtaining stable fragments and domain constructs essential for structural studies by crystallography and NMR

Inhibition of Competence Development, Horizontal Gene Transfer and Virulence in Streptococcus pneumoniae by a Modified Competence Stimulating Peptide

Competence stimulating peptide (CSP) is a 17-amino acid peptide pheromone secreted by Streptococcus pneumoniae. Upon binding of CSP to its membrane-associated receptor kinase ComD, a cascade of signaling events is initiated, leading to activation of the competence regulon by the response regulator ComE. Genes encoding proteins that are involved in DNA uptake and transformation, as well as virulence, are upregulated. Previous studies have shown that disruption of key components in the competence regulon inhibits DNA transformation and attenuates virulence. Thus, synthetic analogues that competitively inhibit CSPs may serve as attractive drugs to control pneumococcal infection and to reduce horizontal gene transfer during infection. We performed amino acid substitutions on conserved amino acid residues of CSP1 in an effort to disable DNA transformation and to attenuate the virulence of S. pneumoniae. One of the mutated peptides, CSP1-E1A, inhibited development of competence in DNA transformation by outcompeting CSP1 in time and concentration-dependent manners. CSP1-E1A reduced the expression of pneumococcal virulence factors choline binding protein D (CbpD) and autolysin A (LytA) in vitro, and significantly reduced mouse mortality after lung infection. Furthermore, CSP1-E1A attenuated the acquisition of an antibiotic resistance gene and a capsule gene in vivo. Finally, we demonstrated that the strategy of using a peptide inhibitor is applicable to other CSP subtype, including CSP2. CSP1-E1A and CSP2-E1A were able to cross inhibit the induction of competence and DNA transformation in pneumococcal strains with incompatible ComD subtypes. These results demonstrate the applicability of generating competitive analogues of CSPs as drugs to control horizontal transfer of antibiotic resistance and virulence genes, and to attenuate virulence during infection by S. pneumoniae

Transcriptomic profiling of host-parasite interactions in the microsporidian <i>Trachipleistophora hominis</i>

BACKGROUND: Trachipleistophora hominis was isolated from an HIV/AIDS patient and is a member of a highly successful group of obligate intracellular parasites. METHODS: Here we have investigated the evolution of the parasite and the interplay between host and parasite gene expression using transcriptomics of T. hominis-infected rabbit kidney cells. RESULTS: T. hominis has about 30 % more genes than small-genome microsporidians. Highly expressed genes include those involved in growth, replication, defence against oxidative stress, and a large fraction of uncharacterised genes. Chaperones are also highly expressed and may buffer the deleterious effects of the large number of non-synonymous mutations observed in essential T. hominis genes. Host expression suggests a general cellular shutdown upon infection, but ATP, amino sugar and nucleotide sugar production appear enhanced, potentially providing the parasite with substrates it cannot make itself. Expression divergence of duplicated genes, including transporters used to acquire host metabolites, demonstrates ongoing functional diversification during microsporidian evolution. We identified overlapping transcription at more than 100 loci in the sparse T. hominis genome, demonstrating that this feature is not caused by genome compaction. The detection of additional transposons of insect origin strongly suggests that the natural host for T. hominis is an insect. CONCLUSIONS: Our results reveal that the evolution of contemporary microsporidian genomes is highly dynamic and innovative. Moreover, highly expressed T. hominis genes of unknown function include a cohort that are shared among all microsporidians, indicating that some strongly conserved features of the biology of these enormously successful parasites remain uncharacterised. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12864-015-1989-z) contains supplementary material, which is available to authorized users