Search

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.

Introduction. The emergence of SARS-CoV-2 has taken humanity off guard. Following an outbreak of SARS-CoV in 2002, and MERS-CoV about 10 years later, SARS-CoV-2 is the third coronavirus in less than 20 years to cross the species barrier and start spreading by human-to-human transmission. It is the most infectious of the three, currently causing the COVID-19 pandemic. No treatment has been approved for COVID-19. We previously proposed targets that can serve as binding sites for antiviral drugs for multiple coronaviruses, and here we set out to find current drugs that can be repurposed as COVID-19 therapeutics. Aim. To identify drugs against COVID-19, we performed an in silico virtual screen with the US Food and Drug Administration (FDA)-approved drugs targeting the RNA-dependent RNA polymerase (RdRP), a critical enzyme for coronavirus replication. Methodology. Initially, no RdRP structure of SARS-CoV-2 was available. We performed basic sequence and structural analysis to determine if RdRP from SARS-CoV was a suitable replacement. We performed molecular dynamics simulations to generate multiple starting conformations that were used for the in silico virtual screen. During this work, a structure of RdRP from SARS-CoV-2 became available and was also included in the in silico virtual screen. Results. The virtual screen identified several drugs predicted to bind in the conserved RNA tunnel of RdRP, where many of the proposed targets were located. Among these candidates, quinupristin is particularly interesting because it is expected to bind across the RNA tunnel, blocking access from both sides and suggesting that it has the potential to arrest viral replication by preventing viral RNA synthesis. Quinupristin is an antibiotic that has been in clinical use for two decades and is known to cause relatively minor side effects. Conclusion. Quinupristin represents a potential anti-SARS-CoV-2 therapeutic. At present, we have no evidence that this drug is effective against SARS-CoV-2 but expect that the biomedical community will expeditiously follow up on our in silico findings.

Prion diseases are infectious and fatal neurodegenerative diseases affecting humans and animals. Transmission is possible within and between species with zoonotic potential. Currently, no prophylaxis or treatment exists. Prions are composed of the misfolded isoform PrPSc of the cellular prion protein PrPC. Expression of PrPC is a prerequisite for prion infection, and conformational conversion of PrPC is induced upon its direct interaction with PrPSc. Inhibition of this interaction can abrogate prion propagation, and we have previously established peptide aptamers (PAs) binding to PrPC as new anti-prion compounds. Here, we mapped the interaction site of PA8 in PrP and modeled the complex in silico to design targeted mutations in PA8 which presumably enhance binding properties. Using these PA8 variants, we could improve PA-mediated inhibition of PrPSc replication and de novo infection of neuronal cells. Furthermore, we demonstrate that binding of PA8 and its variants increases PrPC α-cleavage and interferes with its internalization. This gives rise to high levels of the membrane-anchored PrP-C1 fragment, a transdominant negative inhibitor of prion replication. PA8 and its variants interact with PrPC at its central and most highly conserved domain, a region which is crucial for prion conversion and facilitates toxic signaling of Aβ oligomers characteristic for Alzheimer’s disease. Our strategy allows for the first time to induce α-cleavage, which occurs within this central domain, independent of targeting the responsible protease. Therefore, interaction of PAs with PrPC and enhancement of α-cleavage represent mechanisms that can be beneficial for the treatment of prion and other neurodegenerative diseases. Electronic supplementary material The online version of this article (10.1007/s12035-018-0944-9) contains supplementary material, which is available to authorized users.

Within the last 15 years, two related coronaviruses (Severe Acute Respiratory Syndrome [SARS]-CoV and Middle East Respiratory Syndrome [MERS]-CoV) expanded their host range to include humans, with increased virulence in their new host. Coronaviruses were recently found to have little intrinsic disorder compared with many other virus families. Because intrinsically disordered regions have been proposed to be important for rewiring interactions between virus and host, we investigated the conservation of intrinsic disorder and secondary structure in coronaviruses in an evolutionary context. We found that regions of intrinsic disorder are rarely conserved among different coronavirus protein families, with the primary exception of the nucleocapsid. Also, secondary structure predictions are only conserved across 50–80% of sites for most protein families, with the implication that 20–50% of sites do not have conserved secondary structure prediction. Furthermore, nonconserved structure sites are significantly less constrained in sequence divergence than either sites conserved in the secondary structure or sites conserved in loop. Avoiding regions symptomatic of conformational flexibility such as disordered sites and sites with nonconserved secondary structure to identify potential broad-specificity antiviral targets, only one sequence motif (five residues or longer) remains from the >10,000 starting sites across all coronaviruses in this study. The identified sequence motif is found within the nonstructural protein (NSP) 12 and constitutes an antiviral target potentially effective against the present day and future coronaviruses. On shorter evolutionary timescales, the SARS and MERS clades have more sequence motifs fulfilling the criteria applied. Interestingly, many motifs map to NSP12 making this a prime target for coronavirus antivirals.

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.

Various structural and functional constraints govern the evolution of protein sequences. As a result, the relative rates of amino acid replacement among sites within a protein can vary significantly. Previous large-scale work on Metazoan (Animal) protein sequence alignments indicated that amino acid replacement rates are partially driven by a complex interaction among three factors: intrinsic disorder propensity, secondary structure, and functional domain involvement. Here, we use sequence-based predictors to evaluate the effects of these factors on site-specific sequence evolutionary rates within four eukaryotic lineages: Metazoans, Plants, Saccharomycete Fungi, and Alveolate Protists. Our results show broad, consistent trends across all four Eukaryote groups. In all four lineages, there is a significant increase in amino acid replacement rates when comparing: (i) disordered vs. ordered sites, (ii) random coil sites vs. sites in secondary structures, and (iii) inter-domain linker sites vs. sites in functional domains. Additionally, within Metazoans, Plants, and Saccharomycetes, there is a strong confounding interaction between intrinsic disorder and secondary structure&mdash, alignment sites exhibiting both high disorder propensity and involvement in secondary structures have very low average rates of sequence evolution. Analysis of gene ontology (GO) terms revealed that in all four lineages, a high fraction of sequences containing these conserved, disordered-structured sites are involved in nucleic acid binding. We also observe notable differences in the statistical trends of Alveolates, where intrinsically disordered sites are more variable than in other Eukaryotes and the statistical interactions between disorder and other factors are less pronounced.

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.

One of the largest multigene families in Metazoa are the tyrosine kinases (TKs). These are important multifunctional proteins that have evolved as dynamic switches that perform tyrosine phosphorylation and other noncatalytic activities regulated by various allosteric mechanisms. TKs interact with each other and with other molecules, ultimately activating and inhibiting different signaling pathways. TKs are implicated in cancer and almost 30 FDA-approved TK inhibitors are available. However, specific binding is a challenge when targeting an active site that has been conserved in multiple protein paralogs for millions of years. A cassette domain (CD) containing SH3–SH2–Tyrosine Kinase domains reoccurs in vertebrate nonreceptor TKs. Although part of the CD function is shared between TKs, it also presents TK specific features. Here, the evolutionary dynamics of sequence, structure, and phosphorylation across the CD in 17 TK paralogs have been investigated in a large-scale study. We establish that TKs often have ortholog-specific structural disorder and phosphorylation patterns, while secondary structure elements, as expected, are highly conserved. Further, domain-specific differences are at play. Notably, we found the catalytic domain to fluctuate more in certain secondary structure elements than the regulatory domains. By elucidating how different properties evolve after gene duplications and which properties are specifically conserved within orthologs, the mechanistic understanding of protein evolution is enriched and regions supposedly critical for functional divergence across paralogs are highlighted.

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.

Protein sequence, structure, and function are inherently linked through evolution and population genetics. Our knowledge of protein structure comes from solved structures in the Protein Data Bank (PDB), our knowledge of sequence through sequences found in the NCBI sequence databases (http://www.ncbi.nlm.nih.gov/), and our knowledge of function through a limited set of in-vitro biochemical studies. How these intersect through evolution is described in the first part of the review. In the second part, our understanding of a series of questions is addressed. This includes how sequences evolve within structures, how evolutionary processes enable structural transitions, how the folding process can change through evolution and what the fitness impacts of this might be. Moving beyond static structures, the evolution of protein kinetics (including normal modes) is discussed, as is the evolution of conformational ensembles and structurally disordered proteins. This ties back to a question of the role of neostructuralization and how it relates to selection on sequences for functions. The relationship between metastability, the fitness landscape, sequence divergence, and organismal effective population size is explored. Lastly, a brief discussion of modeling the evolution of sequences of ordered and disordered proteins is entertained.

Protein structure is commonly regarded to be conserved and to dictate function. Most proteins rely on conformational flexibility to some degree. Are regions that convey conformational flexibility conserved over evolutionary time? Can changes in conformational flexibility alter protein function? Here, the evolutionary dynamics of structurally ordered and disordered (flexible) regions are investigated genome-wide in flaviviruses, revealing that the amount and location of structural disorder fluctuates highly among related proteins. Some regions are prone to shift between structured and flexible states. Increased evolutionary dynamics of structural disorder is observed for some lineages but not in others. Lineage-specific transitions of this kind could alter the conformational ensemble accessible to the same protein in different species, causing a functional change, even if the predominant function remains conserved. Thus, rapid evolutionary dynamics of structural disorder is a potential driving force for phenotypic divergence among flaviviruses.