The variation among sites of protein structure divergence is shaped by mutation and scaled by selection.
Ontology highlight
ABSTRACT:
SUBMITTER: Marcos ML
PROVIDER: S-EPMC8244499 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:We present a novel method to identify sites under selection in protein-coding genes. Our method combines the traditional Goldman-Yang model of coding-sequence evolution with the information obtained from the 3D structure of the evolving protein, specifically the relative solvent accessibility (RSA) of individual residues. We develop a random-effects likelihood sites model in which rate classes are RSA dependent. The RSA dependence is modeled with linear functions. We demonstrate that our RSA-dependent model provides a significantly better fit to molecular sequence data than does a traditional, RSA-independent model. We further show that our model provides a natural, RSA-dependent neutral baseline for the evolutionary rate ratio ω = dN/dS Sites that deviate from this neutral baseline likely experience selection pressure for function. We apply our method to the influenza proteins hemagglutinin and neuraminidase. For hemagglutinin, our method recovers positively selected sites near the sialic acid-binding site and negatively selected sites that may be important for trimerization. For neuraminidase, our method recovers the oseltamivir resistance site and otherwise suggests that few sites deviate from the neutral baseline. Our method is broadly applicable to any protein sequences for which structural data are available or can be obtained via homology modeling or threading.
| S-EPMC3525147 | biostudies-literature
Project description:Human mtDNA shows striking regional variation, traditionally attributed to genetic drift. However, it is not easy to account for the fact that only two mtDNA lineages (M and N) left Africa to colonize Eurasia and that lineages A, C, D, and G show a 5-fold enrichment from central Asia to Siberia. As an alternative to drift, natural selection might have enriched for certain mtDNA lineages as people migrated north into colder climates. To test this hypothesis we analyzed 104 complete mtDNA sequences from all global regions and lineages. African mtDNA variation did not significantly deviate from the standard neutral model, but European, Asian, and Siberian plus Native American variations did. Analysis of amino acid substitution mutations (nonsynonymous, Ka) versus neutral mutations (synonymous, Ks) (kaks) for all 13 mtDNA protein-coding genes revealed that the ATP6 gene had the highest amino acid sequence variation of any human mtDNA gene, even though ATP6 is one of the more conserved mtDNA proteins. Comparison of the kaks ratios for each mtDNA gene from the tropical, temperate, and arctic zones revealed that ATP6 was highly variable in the mtDNAs from the arctic zone, cytochrome b was particularly variable in the temperate zone, and cytochrome oxidase I was notably more variable in the tropics. Moreover, multiple amino acid changes found in ATP6, cytochrome b, and cytochrome oxidase I appeared to be functionally significant. From these analyses we conclude that selection may have played a role in shaping human regional mtDNA variation and that one of the selective influences was climate.
| S-EPMC140917 | biostudies-literature
Project description:To develop new strategies for searching for genetic associations with complex human diseases, we analyzed 2784 single-nucleotide polymorphisms (SNPs) in 396 protein-coding genes involved in biological processes relevant to cancer and other complex diseases, with respect to gene diversity within samples of individuals representing the three major historic human populations (African, European, and Asian) and with respect to interpopulation genetic distance. Reduced levels of both intrapopulation gene diversity and interpopulation genetic distance were seen in the case of SNPs located within the 5'-UTR and at nonsynonymous SNPs, causing radical changes to protein structure. Reduction of gene diversity at SNP loci in these categories was evidence of purifying selection acting at these sites, which in turn causes a reduction in interpopulation divergence. By contrast, a small number of SNP sites in these categories revealed unusually high genetic distances between the two most diverged populations (African and Asian); these loci may have historically been subject to divergent selection pressures.
| S-EPMC1242184 | biostudies-literature
Project description:Nuclear magnetic resonance paramagnetic relaxation enhancement (PRE) measures long-range distances to isotopically labeled residues, providing useful constraints for protein structure prediction. The method usually requires labor-intensive conjugation of nitroxide labels to multiple locations on the protein, one at a time. Here a computational procedure, based on protein sequence and simple secondary structure models, is presented to facilitate optimal placement of a minimum number of labels needed to determine the correct topology of a helical transmembrane protein. Tests on DsbB (four helices) using just one label lead to correct topology predictions in four of five cases, with the predicted structures <6 Å to the native structure. Benchmark results using simulated PRE data show that we can generally predict the correct topology for five and six to seven helices using two and three labels, respectively, with an average success rate of 76% and structures of similar precision. The results show promise in facilitating experimentally constrained structure prediction of membrane proteins.
| S-EPMC3099474 | biostudies-literature
Project description:Numerous approaches have been developed to infer natural selection based on the comparison of polymorphism within species and divergence between species. These methods are especially powerful for the detection of uniform selection operating across a gene. However, empirical analyses have demonstrated that regions of protein-coding genes exhibiting clusters of amino acid substitutions are subject to different levels of selection relative to other regions of the same gene. To quantify this heterogeneity of selection within coding sequences, we developed Model Averaged Site Selection via Poisson Random Field (MASS-PRF). MASS-PRF identifies an ensemble of intragenic clustering models for polymorphic and divergent sites. This ensemble of models is used within the Poisson Random Field framework to estimate selection intensity on a site-by-site basis. Using simulations, we demonstrate that MASS-PRF has high power to detect clusters of amino acid variants in small genic regions, can reliably estimate the probability of a variant occurring at each nucleotide site in sequence data and is robust to historical demographic trends and recombination. We applied MASS-PRF to human gene polymorphism derived from the 1,000 Genomes Project and divergence data from the common chimpanzee. On the basis of this analysis, we discovered striking regional variation in selection intensity, indicative of positive or negative selection, in well-defined domains of genes that have previously been associated with neurological processing, immunity, and reproduction. We suggest that amino acid-altering substitutions within these regions likely are or have been selectively advantageous in the human lineage, playing important roles in protein function.
| S-EPMC5850860 | biostudies-literature
Project description:The oligomerization of amyloid beta (A?) peptides into soluble non-fibrillar species plays a critical role in the pathogenesis of Alzheimer's disease. However, it has been challenging to characterize the tertiary and quaternary structures of A? peptides due to their disordered nature and high aggregation propensity. In this work, replica exchange molecular dynamics simulations were used to explore the conformational space of A?42 monomer. Among the most populated transient states, we identified a particular conformation which was able to generate ring-shaped pentamers and hexamers, when docked onto itself. The structures of these aggregates were stable during microsecond all-atom MD simulations in explicit solvent. In addition to high resolution models of these oligomers, this study provides support for the conformational selection mechanism of A? peptide self-assembly.
| S-EPMC4751476 | biostudies-literature
Project description:DNA diversity varies across the genome of many species. Variation in diversity across a genome might arise from regional variation in the mutation rate, variation in the intensity and mode of natural selection, and regional variation in the recombination rate. We show that both noncoding and nonsynonymous diversity are positively correlated to a measure of the mutation rate and the recombination rate and negatively correlated to the density of conserved sequences in 50 kb windows across the genomes of humans and nonhuman homininae. Interestingly, we find that although noncoding diversity is equally affected by these three genomic variables, nonsynonymous diversity is mostly dominated by the density of conserved sequences. The positive correlation between diversity and our measure of the mutation rate seems to be largely a direct consequence of regions with higher mutation rates having more diversity. However, the positive correlation with recombination rate and the negative correlation with the density of conserved sequences suggest that selection at linked sites also affect levels of diversity. This is supported by the observation that the ratio of the number of nonsynonymous to noncoding polymorphisms is negatively correlated to a measure of the effective population size across the genome. We show these patterns persist even when we restrict our analysis to GC-conservative mutations, demonstrating that the patterns are not driven by GC biased gene conversion. In conclusion, our comparative analyses describe how recombination rate, gene density, and mutation rate interact to produce the patterns of DNA diversity that we observe along the hominine genomes.
| S-EPMC6944223 | biostudies-literature
Project description:Recent work has demonstrated that evolutionary processes shape ecological dynamics on relatively short timescales (eco-evolutionary dynamics), but demonstrating these effects at large spatial scales in natural landscapes has proven difficult. We used empirical studies and modeling to investigate how selective pressures from fire and predispersal seed predation affect the evolution of serotiny, an ecologically important trait. Serotiny is a highly heritable key reproductive trait in Rocky Mountain lodgepole pine (Pinus contorta subsp. latifolia), a conifer that dominates millions of hectares in western North America. In these forests, the frequency of serotiny determines postfire seedling density with corresponding community- and ecosystem-level effects. We found that serotinous individuals have a selective advantage at high fire frequencies and low predation pressure; however, very high seed predation shifted the selective advantage to nonserotinous individuals even at high fire frequencies. Simulation modeling suggests that spatial variation in the frequency of serotiny results from heterogeneity in these two selective agents. These results, combined with previous findings showing a negative association between the density of seed predators and the frequency of serotiny at both landscape and continental scales, demonstrate that contemporary patterns in serotiny reflect an evolutionary response to conflicting selection pressures from fire and seed predation. Thus, we show that variation in the frequency of a heritable polygenic trait depends on spatial variation in two dominant selective agents, and, importantly, the effects of the local trait variation propagate with profound consequences to the structure and function of communities and ecosystems across a large landscape.
| S-EPMC4084486 | biostudies-literature
Project description:The distance-dependent structure-derived potentials developed so far all employed a reference state that can be characterized as a residue (atom)-averaged state. Here, we establish a new reference state called the distance-scaled, finite ideal-gas reference (DFIRE) state. The reference state is used to construct a residue-specific all-atom potential of mean force from a database of 1011 nonhomologous (less than 30% homology) protein structures with resolution less than 2 A. The new all-atom potential recognizes more native proteins from 32 multiple decoy sets, and raises an average Z-score by 1.4 units more than two previously developed, residue-specific, all-atom knowledge-based potentials. When only backbone and C(beta) atoms are used in scoring, the performance of the DFIRE-based potential, although is worse than that of the all-atom version, is comparable to those of the previously developed potentials on the all-atom level. In addition, the DFIRE-based all-atom potential provides the most accurate prediction of the stabilities of 895 mutants among three knowledge-based all-atom potentials. Comparison with several physical-based potentials is made.
| S-EPMC2373736 | biostudies-literature
Project description:Mutation generates the heritable variation that genetic drift and natural selection shape. In classical quantitative genetic models, drift is a function of the effective population size and acts uniformly across traits, whereas mutation and selection act trait-specifically. We identified thousands of quantitative trait loci (QTLs) influencing transcript abundance traits in a cross of two Caenorhabditis elegans strains; although trait-specific mutation and selection explained some of the observed pattern of QTL distribution, the pattern was better explained by trait-independent variation in the intensity of selection on linked sites. Our results suggest that traits in C. elegans exhibit different levels of variation less because of their own attributes than because of differences in the effective population sizes of the genomic regions harboring their underlying loci.
| S-EPMC3138179 | biostudies-literature