Search

Protein-protein interactions drive functions in eukaryotes that can be described by short linear motifs (SLiMs). Conservation of SLiMs help illuminate functional SLiMs in eukaryotic protein families. However, the simplicity of eukaryotic SLiMs makes them appear by chance due to mutational processes not only in eukaryotes but also in pathogenic bacteria and viruses. Further, functional eukaryotic SLiMs are often found in disordered regions. Although proteomes from pathogenic bacteria and viruses have less disorder than eukaryotic proteomes, their proteins can successfully mimic eukaryotic SLiMs and disrupt host cellular function. Identifying important SLiMs in pathogens is difficult but essential for understanding potential host-pathogen interactions. We performed a comparative analysis of structural features for experimentally verified SLiMs from the Eukaryotic Linear Motif (ELM) database across viruses, bacteria, and eukaryotes. Our results revealed that many viral SLiMs and specific motifs found across viruses and eukaryotes, such as some glycosylation motifs, have less disorder. Analyzing the disorder and coil properties of equivalent SLiMs from pathogens and eukaryotes revealed that some motifs are more structured in pathogens than their eukaryotic counterparts and vice versa. These results support a varying mechanism of interaction between pathogens and their eukaryotic hosts for some of the same motifs.

Heterotachy-the change in sequence evolutionary rate over time-is a common feature of protein molecular evolution. Decades of studies have shed light on the conditions under which heterotachy occurs, and there is evidence that site-specific evolutionary rate shifts are correlated with changes in protein function. Here, we present a large-scale, computational analysis using thousands of protein sequence alignments from animal and plant proteomes, representing genes related either by orthology (speciation events) or paralogy (gene duplication), to compare sequence divergence patterns in orthologous vs. paralogous sequence alignments. We use sequence-based phylogenetic analyses to infer overall sequence divergence (tree length/number of sequences) and to fit site-specific rates to a discrete gamma distribution with a shape parameter α. This inference method is applied to real protein sequence alignments, as well as alignments simulated under various models of protein sequence evolution. Our simulations indicate that sequence divergence and the α parameter are positively correlated when sequences evolve with heterotachy, meaning that inferred site rate distributions appear more uniform as sequences diverge. Divergence and α are also positively correlated in both orthologous and paralogous genes, but the average increase in α (as a function of divergence) is significantly higher in paralogous protein alignments than in orthologous alignments. This result is consistent with the widely held view that recently duplicated proteins initially evolve under relaxed selective pressure, promoting functional divergence by accumulation of amino acid replacements, and hence experience more evolutionary rate fluctuations than orthologous proteins. We discuss these findings in the context of the ortholog conjecture, a long-standing assumption in molecular evolution, which posits that protein sequences related by orthology tend to be more functionally conserved than paralogous proteins.

Dengue virus (DENV) challenges vaccine design due to antibody-dependent enhancement (ADE) and evidence suggests that Zika virus (ZIKV) experiences ADE with DENV and West Nile virus (WNV) antibodies. Thus, multiple flaviviruses must be considered when developing novel therapies against ZIKV. We analyzed 42 flavivirus polyproteins in their evolutionary context to identify motifs conserved in sequence with low real-time and evolutionary conformational flexibility, thought to be fitness-critical sites. We also analyzed evolutionary rate-shifts between clades for insight on vector specificity. For mosquito-borne flaviviruses, two conserved motifs were identified within the RNA-dependent RNA polymerase (RdRP), critical for flavivirus genome replication. Clade-specific motifs were identified for the ZIKV+DENV and WNV clades, many of which were also in RdRP. Six sites in motifs for WNV experienced significant evolutionary rate-shifts, suggesting their importance for functional divergence. Overall, some of these motifs are prime candidates as broadly neutralizing antiviral drug targets across different mosquito-borne flaviviruses.

Protein evolution often occurs at unequal rates in different sites along an amino acid chain. Site-specific evolutionary rates have been linked to several structural and functional properties of proteins. Previous analyses of this phenomenon have involved relatively small datasets and, in some cases, the interaction among multiple structural factors is not evaluated. Here, we present the results of a large-scale phylogenetic and statistical analysis, testing the effects and interactions of three structural properties on amino acid replacement rates. We used sequence-based computational methods to predict (i) intrinsic disorder propensity, (ii) secondary structure, and (iii) functional domain involvement across millions of amino acid sites in thousands of sequence alignments of metazoan proteins. Our results somewhat corroborate earlier findings that intrinsically disordered sites tend to be more variable than ordered sites, but there is considerable overlap among their rate distributions, and a significant confounding interaction exists between intrinsic disorder and secondary structure. Notably, protein sites that are consistently predicted to be both intrinsically disordered and involved in secondary structures tend to be the most conserved at the amino acid level, suggesting that they are highly constrained and functionally important. In addition, a significant interaction exists between functional domain involvement and secondary structure. These findings suggest that multiple structural drivers of protein evolution should be evaluated simultaneously in order to get a clear picture of their individual effects as well as any confounding interactions among them.

Conformational flexibility conferred though regions of intrinsic structural disorder allows proteins to behave as dynamic molecules. While it is well-known that intrinsically disordered regions can undergo disorder-to-order transitions in real-time as part of their function, we also are beginning to learn more about the dynamics of disorder-to-order transitions along evolutionary time-scales. Intrinsically disordered regions endow proteins with functional promiscuity, which is further enhanced by the ability of some of these regions to undergo real-time disorder-to-order transitions. Disorder content affects gene retention after whole genome duplication, but it is not necessarily conserved. Altered patterns of disorder resulting from evolutionary disorder-to-order transitions indicate that disorder evolves to modify function through refining stability, regulation, and interactions. Here, we review the evolution of intrinsically disordered regions in eukaryotic proteins. We discuss the interplay between secondary structure and disorder on evolutionary time-scales, the importance of disorder for eukaryotic proteome expansion and functional divergence, and the evolutionary dynamics of disorder.

Despite high sequence identity among mammalian prion proteins, mammals have varying rates of susceptibility to prion disease resulting in a so-called species barrier. The species barrier follows no clear pattern, with closely related species or similar sequences being no more likely to infect each other, and remains an unresolved enigma. Variation of the conformationally flexible regions may alter the thermodynamics of the conformational change, commonly referred to as the conformational conversion, that occurs in the pathogenic process of the mammalian prion protein. A conformational ensemble scenario is supported by the species barrier in prion disease and evidence that there are strains of pathogenic prion with different conformations within species. To study how conformational flexibility has evolved in the prion protein, an investigation was undertaken of the evolutionary dynamics of structurally disordered regions in the mammalian prion protein, non-mammalian prion protein that is not vulnerable to prion disease, and remote homologs Doppel and Shadoo. Structural disorder prediction analyzed in an evolutionary context revealed that the occurrence of increased or altered conformational flexibility in mammalian prion proteins coincides with key events among PrP, Doppel, and Shadoo. Comparatively rapid evolutionary dynamics of conformational flexibility in the prion protein suggest that the species barrier is not a static phenomenon. A small number of amino acid substitutions can repopulate the conformational ensemble and have a disproportionately large effect on pathogenesis.

Conformational and functional flexibility promote protein evolvability. High evolvability allows related proteins to functionally diverge and perhaps to neostructuralize. p53 is a multifunctional protein frequently referred to as the Guardian of the Genome–a hub for e.g. incoming and outgoing signals in apoptosis and DNA repair. p53 has been found to be structurally disordered, an extreme form of conformational flexibility. Here, p53, and its paralogs p63 and p73, were studied for further insights into the evolutionary dynamics of structural disorder, secondary structure, and phosphorylation. This study is focused on the post gene duplication phase for the p53 family in vertebrates, but also visits the origin of the protein family and the early domain loss and gain events. Functional divergence, measured by rapid evolutionary dynamics of protein domains, structural properties, and phosphorylation propensity, is inferred across vertebrate p53 proteins, in p63 and p73 from fish, and between the three paralogs. In particular, structurally disordered regions are redistributed among paralogs, but within clades redistribution of structural disorder also appears to be an ongoing process. Despite its deemed importance as the Guardian of the Genome, p53 is indeed a protein with high evolvability as seen not only in rearranged structural disorder, but also in fluctuating domain sequence signatures among lineages.

Cyclic nucleotides, cAMP and cGMP, are ubiquitous second messengers that regulate metabolic and behavioral responses in diverse organisms. We describe purification, engineering, and characterization of photoactivated nucleotidyl cyclases that can be used to manipulate cAMP and cGMP levels in vivo. We identified the blaC gene encoding a putative photoactivated adenylyl cyclase in the Beggiatoa sp. PS genome. BlaC contains a BLUF domain involved in blue-light sensing using FAD and a nucleotidyl cyclase domain. The blaC gene was overexpressed in Escherichia coli, and its product was purified. Irradiation of BlaC in vitro resulted in a small red shift in flavin absorbance, typical of BLUF photoreceptors. BlaC had adenylyl cyclase activity that was negligible in the dark and up-regulated by light by 2 orders of magnitude. To convert BlaC into a guanylyl cyclase, we constructed a model of the nucleotidyl cyclase domain and mutagenized several residues predicted to be involved in substrate binding. One triple mutant, designated BlgC, was found to have photoactivated guanylyl cyclase in vitro. Irradiation with blue light of the E. coli cya mutant expressing BlaC or BlgC resulted in the significant increases in cAMP or cGMP synthesis, respectively. BlaC, but not BlgC, restored cAMP-dependent growth of the mutant in the presence of light. Small protein sizes, negligible activities in the dark, high light-to-dark activation ratios, functionality at broad temperature range and physiological pH, as well as utilization of the naturally occurring flavins as chromophores make BlaC and BlgC attractive for optogenetic applications in various animal and microbial models.

Protein structure is generally more conserved than sequence, but for regions that can adopt different structures in different environments, does this hold true? Understanding how structurally disordered regions evolve altered secondary structure element propensities as well as conformational flexibility among paralogs are fundamental questions for our understanding of protein structural evolution. We have investigated the evolutionary dynamics of structural disorder in protein families containing both orthologs and paralogs using phylogenetic tree reconstruction, protein structure disorder prediction, and secondary structure prediction in order to shed light upon these questions. Our results indicate that the extent and location of structurally disordered regions are not universally conserved. As structurally disordered regions often have high conformational flexibility, this is likely to have an effect on how protein structure evolves as spatially altered conformational flexibility can also change the secondary structure propensities for homologous regions in a protein family.

A fundamental assumption in biology is that protein structure is more conserved than protein sequence. This opinion stems from observations of the fold distribution of protein structures present in the protein data bank (PDB), where homologous proteins are found to display the same fold. However, the set of proteins in PDB is not representative of protein sequence space. The predominant experimental method used for the structural characterization of the proteins in PDB is X-ray crystallography. Thus, conformational flexibility is lost and a static snapshot of a protein’s dynamic structure remains. Using structural disorder as a proxy for protein dynamism, it seems that at least 30% of eukaryotic proteins have dynamic regions, that is, regions with conformational flexibility. Little is known about the conservation of structurally disordered regions, except that they often evolve at an elevated rate compared to the structured regions. Our preliminary results indicate that there are regions that show fast ev...