Structure and age jointly influence rates of protein evolution.
Ontology highlight
ABSTRACT:
SUBMITTER: Toll-Riera M
PROVIDER: S-EPMC3364943 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A general understanding of the complex phenomenon of protein evolution requires the accurate description of the constraints that define the sub-space of proteins with mutations that do not appreciably reduce the fitness of the organism. Such constraints can have multiple origins, in this work we present a model for constrained evolutionary trajectories represented by a markovian process throughout a set of protein-like structures artificially constructed to be topological intermediates between the structure of two natural occurring proteins. The number and type of intermediate steps defines how constrained the total evolutionary process is. By using a coarse-grained representation for the protein structures, we derive an analytic formulation of the transition rates between each of the intermediate structures. The results indicate that compact structures with a high number of hydrogen bonds are more probable and have a higher likelihood to arise during evolution. Knowledge of the transition rates allows for the study of complex evolutionary pathways represented by trajectories through a set of intermediate structures.
| S-EPMC3367531 | biostudies-literature
Project description:Avian influenza viruses (AIVs) have been pivotal to the origination of human pandemic strains. Despite their scientific and public health significance, however, there remains much to be understood about the ecology and evolution of AIVs in wild birds, where major pools of genetic diversity are generated and maintained. Here, we present comparative phylodynamic analyses of human and AIVs in North America, demonstrating (i) significantly higher standing genetic diversity and (ii) phylogenetic trees with a weaker signature of immune escape in AIVs than in human viruses. To explain these differences, we performed statistical analyses to quantify the relative contribution of several potential explanations. We found that HA genetic diversity in avian viruses is determined by a combination of factors, predominantly subtype-specific differences in host immune selective pressure and the ecology of transmission (in particular, the durability of subtypes in aquatic environments). Extending this analysis using a computational model demonstrated that virus durability may lead to long-term, indirect chains of transmission that, when coupled with a short host lifespan, can generate and maintain the observed high levels of genetic diversity. Further evidence in support of this novel finding was found by demonstrating an association between subtype-specific environmental durability and predicted phylogenetic signatures: genetic diversity, variation in phylogenetic tree branch lengths, and tree height. The conclusion that environmental transmission plays an important role in the evolutionary biology of avian influenza viruses-a manifestation of the "storage effect"-highlights the potentially unpredictable impact of wildlife reservoirs for future human pandemics and the need for improved understanding of the natural ecology of these viruses.
| S-EPMC4130664 | biostudies-literature
Project description:The number of amino acids that occupy a given protein site during evolution reflects the selective constraints operating on the site. This evolutionary variability is strongly influenced by the structural properties of the site in the native structure, and it is quantified either through sequence entropy or through substitution rates. However, while the sequence entropy only depends on the equilibrium frequencies of the amino acids, the substitution rate also depends on the exchangeability matrix that describes mutations in the mathematical model of the substitution process. Here we apply two variants of a mathematical model of protein evolution with selection for protein stability, both against unfolding and against misfolding. Exploiting the approximation of independent sites, these models allow computing site-specific substitution processes that satisfy global constraints on folding stability. We find that site-specific substitution rates do not depend only on the selective constraints acting on the site, quantified through its sequence entropy. In fact, polar sites evolve faster than hydrophobic sites even for equal sequence entropy, as a consequence of the fact that polar amino acids are characterized by higher mutational exchangeability than hydrophobic ones. Accordingly, the model predicts that more polar proteins tend to evolve faster. Nevertheless, these results change if we compare proteins that evolve under different mutation biases, such as orthologous proteins in different bacterial genomes. In this case, the substitution rates are faster in genomes that evolve under mutational bias that favor hydrophobic amino acids by preferentially incorporating the nucleotide Thymine that is more frequent in hydrophobic codons. This appearingly contradictory result arises because buried sites occupied by hydrophobic amino acids are characterized by larger selective factors that largely amplify the substitution rate between hydrophobic amino acids, while the selective factors of exposed sites have a weaker effect. Thus, changes in the mutational bias produce deep effects on the biophysical properties of the protein (hydrophobicity) and on its evolutionary properties (sequence entropy and substitution rate) at the same time. The program Prot_evol that implements the two site-specific substitution processes is freely available at https://ub.cbm.uam.es/prot_fold_evol/prot_fold_evol_soft_main.php#Prot_Evol.
| S-EPMC6174885 | biostudies-literature
Project description:In the human genome, most genes undergo splicing, and patterns of codon usage are splicing dependent: guanine and cytosine (GC) content is the highest within single-exon genes and within first exons of multi-exon genes. However, the effects of codon usage on gene expression are typically characterized in unspliced model genes. Here, we measured the effects of splicing on expression in a panel of synonymous reporter genes that varied in nucleotide composition. We found that high GC content increased protein yield, mRNA yield, cytoplasmic mRNA localization, and translation of unspliced reporters. Splicing did not affect the expression of GC-rich variants. However, splicing promoted the expression of AT-rich variants by increasing their steady-state protein and mRNA levels, in part through promoting cytoplasmic localization of mRNA. We propose that splicing promotes the nuclear export of AU-rich mRNAs and that codon- and splicing-dependent effects on expression are under evolutionary pressure in the human genome.
| S-EPMC7181179 | biostudies-literature
Project description:The evolutionary rates of proteins vary over several orders of magnitude. Recent work suggests that analysis of large data sets of evolutionary rates in conjunction with the results from high-throughput functional genomic experiments can identify the factors that cause proteins to evolve at such dramatically different rates. To this end, we estimated the evolutionary rates of >3,000 proteins in four species of the yeast genus Saccharomyces and investigated their relationship with levels of expression and protein dispensability. Each protein's dispensability was estimated by the growth rate of mutants deficient for the protein. Our analyses of these improved evolutionary and functional genomic data sets yield three main results. First, dispensability and expression have independent, significant effects on the rate of protein evolution. Second, measurements of expression levels in the laboratory can be used to filter data sets of dispensability estimates, removing variates that are unlikely to reflect real biological effects. Third, structural equation models show that although we may reasonably infer that dispensability and expression have significant effects on protein evolutionary rate, we cannot yet accurately estimate the relative strengths of these effects.
| S-EPMC555735 | biostudies-literature
Project description:Although the folding rates of proteins have been studied extensively, both experimentally and theoretically, and many native state topological parameters have been proposed to correlate with or predict these rates, unfolding rates have received much less attention. Moreover, unfolding rates have generally been thought either to not relate to native topology in the same manner as folding rates, perhaps depending on different topological parameters, or to be more difficult to predict. Using a dataset of 108 proteins including two-state and multistate folders, we find that both unfolding and folding rates correlate strongly, and comparably well, with well-established measures of native topology, the absolute contact order and the long range order, with correlation coefficient values of 0.75 or higher. In addition, compared to folding rates, the absolute values of unfolding rates vary more strongly with native topology, have a larger range of values, and correlate better with thermodynamic stability. Similar trends are observed for subsets of different protein structural classes. Taken together, these results suggest that choosing a scaffold for protein engineering may require a compromise between a simple topology that will fold sufficiently quickly but also unfold quickly, and a complex topology that will unfold slowly and hence have kinetic stability, but fold slowly. These observations, together with the established role of kinetic stability in determining resistance to thermal and chemical denaturation as well as proteases, have important implications for understanding fundamental aspects of protein unfolding and folding and for protein engineering and design.
| S-EPMC4380988 | biostudies-literature
Project description:The hepatitis C virus (HCV) infects close to 200 million people globally, resulting in a significant need for effective HCV therapies. The HCV polymerase (gene product NS5B) is a valuable target for therapeutics because of its role in replicating the viral genome. Various studies have identified inhibitors for this enzyme, including non-nucleoside inhibitors (NNIs) that bind distal to the enzyme active site. Recently, it has been shown that simultaneously challenging the enzyme with two NNIs results in enhanced inhibition relative to that observed after challenge with individual inhibitors, suggesting that employing multiple NNIs might be the basis of more effective therapeutics. Nevertheless, the molecular mechanisms responsible for this enhanced inhibition remain unclear. We employ molecular dynamics simulations to determine the origin of enhanced inhibition when two NNIs bind to NS5B. Our results suggest that nonoverlapping NNI sites are compatible with simultaneous binding of dual NNIs. We observe that both inhibitors act in concert to induce novel enzyme conformations and dynamics, allowing us to identify molecular mechanisms underlying enhanced inhibition of NS5B. This knowledge will be useful in optimizing combinations of NNIs to target NS5B, helping to prevent the acquisition of viral resistance that remains a significant barrier to the development of HCV therapeutics.
| S-EPMC4918089 | biostudies-literature
Project description:Coevolution with phages drive the evolution of high bacterial mutation rates in vitro, but the relevance of this finding to natural populations is unclear. Here, we investigated how coevolution affects mutation rate evolution in soil, in the presence and absence of the rest of the natural microbial community. Although mutation rate on average increased threefold, neither coevolving phages nor the rest of natural community significantly affected mutation rates. Our results suggest that features of the soil over and above directly interacting organisms constrain the evolution of strong mutators, helping to explain their relatively low frequency compared with some laboratory and clinical settings.
| S-EPMC3806259 | biostudies-literature
Project description:BackgroundMost studies of molecular evolution are focused on individual genes and proteins. However, understanding the design principles and evolutionary properties of molecular networks requires a system-wide perspective. In the present work we connect molecular evolution on the gene level with system properties of a cellular metabolic network. In contrast to protein interaction networks, where several previous studies investigated the molecular evolution of proteins, metabolic networks have a relatively well-defined global function. The ability to consider fluxes in a metabolic network allows us to relate the functional role of each enzyme in a network to its rate of evolution.ResultsOur results, based on the yeast metabolic network, demonstrate that important evolutionary processes, such as the fixation of single nucleotide mutations, gene duplications, and gene deletions, are influenced by the structure and function of the network. Specifically, central and highly connected enzymes evolve more slowly than less connected enzymes. Also, enzymes carrying high metabolic fluxes under natural biological conditions experience higher evolutionary constraints. Genes encoding enzymes with high connectivity and high metabolic flux have higher chances to retain duplicates in evolution. In contrast to protein interaction networks, highly connected enzymes are no more likely to be essential compared to less connected enzymes.ConclusionThe presented analysis of evolutionary constraints, gene duplication, and essentiality demonstrates that the structure and function of a metabolic network shapes the evolution of its enzymes. Our results underscore the need for systems-based approaches in studies of molecular evolution.
| S-EPMC1779518 | biostudies-literature
Project description:Batrachochytrium salamandrivorans (Bsal) is an emerging invasive pathogen that is highly pathogenic to salamander species. Modeling infection dynamics in this system can facilitate proactive efforts to mitigate this pathogen's impact on North American species. Given its widespread distribution and high abundance, the eastern newt (Notophthalmus viridescens) has the potential to significantly influence Bsal epidemiology. We designed experiments to 1) estimate contact rates given different host densities and habitat structure and 2) estimate the probability of transmission from infected to susceptible individuals. Using parameter estimates from data generated during these experiments, we modeled infection and disease outcomes for a population of newts using a system of differential equations. We found that host contact rates were density-dependent, and that adding habitat structure reduced contacts. The probability of Bsal transmission given contact between newts was very high (>90%) even at early stages of infection. Our simulations show rapid transmission of Bsal among individuals following pathogen introduction, with infection prevalence exceeding 90% within one month and >80% mortality of newts in three months. Estimates of basic reproductive rate (R0) of Bsal for eastern newts were 1.9 and 3.2 for complex and simple habitats, respectively. Although reducing host density and increasing habitat complexity might decrease transmission, these management strategies may be ineffective at stopping Bsal invasion in eastern newt populations due to this species' hyper-susceptibility.
| S-EPMC7101388 | biostudies-literature