Browse
Submit Data
Databases
API
Help

Dataset Information

26 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

An experimentally determined evolutionary model dramatically improves phylogenetic fit.

ABSTRACT: All modern approaches to molecular phylogenetics require a quantitative model for how genes evolve. Unfortunately, existing evolutionary models do not realistically represent the site-heterogeneous selection that governs actual sequence change. Attempts to remedy this problem have involved augmenting these models with a burgeoning number of free parameters. Here, I demonstrate an alternative: Experimental determination of a parameter-free evolutionary model via mutagenesis, functional selection, and deep sequencing. Using this strategy, I create an evolutionary model for influenza nucleoprotein that describes the gene phylogeny far better than existing models with dozens or even hundreds of free parameters. Emerging high-throughput experimental strategies such as the one employed here provide fundamentally new information that has the potential to transform the sensitivity of phylogenetic and genetic analyses.

SUBMITTER: Bloom JD

PROVIDER: S-EPMC4104320 | biostudies-literature | 2014 Aug

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

An experimentally determined evolutionary model dramatically improves phylogenetic fit.

Bloom Jesse D JD

Molecular biology and evolution 20140524 8

All modern approaches to molecular phylogenetics require a quantitative model for how genes evolve. Unfortunately, existing evolutionary models do not realistically represent the site-heterogeneous selection that governs actual sequence change. Attempts to remedy this problem have involved augmenting these models with a burgeoning number of free parameters. Here, I demonstrate an alternative: Experimental determination of a parameter-free evolutionary model via mutagenesis, functional selection, ...[more]

PMID: 24859245

Similar Datasets

An experimentally informed evolutionary model improves phylogenetic fit to divergent lactamase homologs.

Project description:Phylogenetic analyses of molecular data require a quantitative model for how sequences evolve. Traditionally, the details of the site-specific selection that governs sequence evolution are not known a priori, making it challenging to create evolutionary models that adequately capture the heterogeneity of selection at different sites. However, recent advances in high-throughput experiments have made it possible to quantify the effects of all single mutations on gene function. I have previously shown that such high-throughput experiments can be combined with knowledge of underlying mutation rates to create a parameter-free evolutionary model that describes the phylogeny of influenza nucleoprotein far better than commonly used existing models. Here, I extend this work by showing that published experimental data on TEM-1 beta-lactamase (Firnberg E, Labonte JW, Gray JJ, Ostermeier M. 2014. A comprehensive, high-resolution map of a gene's fitness landscape. Mol Biol Evol. 31:1581-1592) can be combined with a few mutation rate parameters to create an evolutionary model that describes beta-lactamase phylogenies much better than most common existing models. This experimentally informed evolutionary model is superior even for homologs that are substantially diverged (about 35% divergence at the protein level) from the TEM-1 parent that was the subject of the experimental study. These results suggest that experimental measurements can inform phylogenetic evolutionary models that are applicable to homologs that span a substantial range of sequence divergence.

| S-EPMC4166927 | biostudies-literature

An experimentally determined protein folding energy landscape.

Project description:Energy landscapes have been used to conceptually describe and model protein folding but have been difficult to measure experimentally, in large part because of the myriad of partly folded protein conformations that cannot be isolated and thermodynamically characterized. Here we experimentally determine a detailed energy landscape for protein folding. We generated a series of overlapping constructs containing subsets of the seven ankyrin repeats of the Drosophila Notch receptor, a protein domain whose linear arrangement of modular structural units can be fragmented without disrupting structure. To a good approximation, stabilities of each construct can be described as a sum of energy terms associated with each repeat. The magnitude of each energy term indicates that each repeat is intrinsically unstable but is strongly stabilized by interactions with its nearest neighbors. These linear energy terms define an equilibrium free energy landscape, which shows an early free energy barrier and suggests preferred low-energy routes for folding.

| S-EPMC521126 | biostudies-literature

Inverse Potts model improves accuracy of phylogenetic profiling.

Project description:Phylogenetic profiling is a powerful computational method for revealing the functions of function-unknown genes. Although conventional similarity metrics in phylogenetic profiling achieved high prediction accuracy, they have two estimation biases: an evolutionary bias and a spurious correlation bias. While previous studies reduced the evolutionary bias by considering a phylogenetic tree, few studies have analyzed the spurious correlation bias. To reduce the spurious correlation bias, we developed metrics based on the inverse Potts model (IPM) for phylogenetic profiling. We also developed a metric based on both the IPM and a phylogenetic tree. In an empirical dataset analysis, we demonstrated that these IPM-based metrics improved the prediction performance of phylogenetic profiling. In addition, we found that the integration of several metrics, including the IPM-based metrics, had superior performance to a single metric. The source code is freely available at https://github.com/fukunagatsu/Ipm. Supplementary data are available at Bioinformatics online.

| S-EPMC8963296 | biostudies-literature

Rigidification Dramatically Improves Inhibitor Selectivity for RAF Kinases.

Project description:One effective means to achieve inhibitor specificity for RAF kinases, an important family of cancer drug targets, has been to target the monomeric inactive state conformation of the kinase domain, which, unlike most other kinases, can accommodate sulfonamide-containing drugs such as vemurafenib and dabrafenib because of the presence of a unique pocket specific to inactive RAF kinases. We previously reported an alternate strategy whereby rigidification of a nonselective pyrazolo[3,4-d]pyrimidine-based inhibitor through ring closure afforded moderate but appreciable increases in selectivity for RAF kinases. Here, we show that a further application of the rigidification strategy to a different pyrazolopyrimidine-based scaffold dramatically improved selectivity for RAF kinases. Crystal structure analysis confirmed our inhibitor design hypothesis revealing that 2l engages an active-like state conformation of BRAF normally associated with poorly discriminating inhibitors. When screened against a panel of distinct cancer cell lines, the optimized inhibitor 2l primarily inhibited the proliferation of the expected BRAFV600E-harboring cell lines consistent with its kinome selectivity profile. These results suggest that rigidification could be a general and powerful strategy for enhancing inhibitor selectivity against protein kinases, which may open up therapeutic opportunities not afforded by other approaches.

| S-EPMC6627727 | biostudies-literature

Antivirals for allosteric inhibition of Zika virus using a homology model and experimentally determined structure of envelope protein.

Project description:An approach to inhibiting enveloped flaviviruses is to deter the ability of the envelope protein(s) binding onto glycoproteins. In previous work, using a small ~100-amino acid homology model of Zika virus envelope protein (ZVEP), we proved the susceptibility of Zika virus to inhibition. In this work, we verify the efficacy of the homology model based antiviral search method using a larger protein (>400 amino acids) and comparing the results with the experimentally determined one (PDB ID:5IRE).By examining how glycan molecules, small-molecule probes and screened ligands that have a high affinity to ZVEP, we report the mechanics of ZVEP to inhibition via allosteric blockage of the glycan-binding domain while proposing even more possibly potent inhibitors. The small molecular probes based study using the homology model and subsequently verified using actual experimental structure, 5IRE, revealed that ZVEP is druggable. A pharmacophore analysis followed by screening showed at least four ligands that allosterically binds to the glycan binding domain constituted by residues VAL 153 and ASN 154 in 5IRE. Based on further selection criteria ZINC40621658 was identified to have high potential to be a strong antiviral candidate for Zika virus inhibition.

| S-EPMC5534091 | biostudies-literature

Model-driven analysis of experimentally determined growth phenotypes for 465 yeast gene deletion mutants under 16 different conditions.

Project description:Understanding the response of complex biochemical networks to genetic perturbations and environmental variability is a fundamental challenge in biology. Integration of high-throughput experimental assays and genome-scale computational methods is likely to produce insight otherwise unreachable, but specific examples of such integration have only begun to be explored.In this study, we measured growth phenotypes of 465 Saccharomyces cerevisiae gene deletion mutants under 16 metabolically relevant conditions and integrated them with the corresponding flux balance model predictions. We first used discordance between experimental results and model predictions to guide a stage of experimental refinement, which resulted in a significant improvement in the quality of the experimental data. Next, we used discordance still present in the refined experimental data to assess the reliability of yeast metabolism models under different conditions. In addition to estimating predictive capacity based on growth phenotypes, we sought to explain these discordances by examining predicted flux distributions visualized through a new, freely available platform. This analysis led to insight into the glycerol utilization pathway and the potential effects of metabolic shortcuts on model results. Finally, we used model predictions and experimental data to discriminate between alternative raffinose catabolism routes.Our study demonstrates how a new level of integration between high throughput measurements and flux balance model predictions can improve understanding of both experimental and computational results. The added value of a joint analysis is a more reliable platform for specific testing of biological hypotheses, such as the catabolic routes of different carbon sources.

| S-EPMC2592718 | biostudies-literature

Experimentally Determined Hansen Solubility Parameters of Biobased and Biodegradable Polyesters.

Project description:Hansen solubility parameters (HSP) of 15 commercially relevant biobased and biodegradable polyesters were experimentally determined by applying a novel approach to the classic solubility study method. In this approach, the extent of swelling in polymer films was determined using a simple equation based on the mass difference between swollen and nonswollen film samples to obtain normalized solvent uptake (N). Using N and HSPiP software, highly accurate HSP values were obtained for all 15 polyesters. Qualitative evaluation of the HSP values was conducted by predicting the miscibility of poly(3-hydroxybutyrate-co-3-hydroxyhexanoate) (PHB-co-HHx, 7 mol % HHx) and poly(lactic acid) (PLA) with a novel lignin-based plasticizer (ethyl 3-(4-ethoxy-3-methoxyphenyl)propanoate, EP) with a relative energy difference (RED) less than 0.4. Additionally, an HSP-predicted plasticizer (di(2-ethylhexyl) adipate, DA) with a larger RED (>0.7) was used to demonstrate the effects of less-miscible additives. Plasticized samples were analyzed by differential scanning calorimetry and polarized optical microscopy (POM) to determine the Tg depression, with EP showing linear Tg depression up to 50% plasticizer loading, whereas DA shows minimal Tg depression past 10% loading. Further analysis by POM reveals that the DA phase separates from both polymers at loadings as low as 2.5% (PHB-co-HHx, 7 mol % HHx) and 5% (PLA), while the EP phase separates at a much higher loading of 50% (PHB-co-HHx, 7 mol% HHx) and 30% (PLA).

| S-EPMC10865435 | biostudies-literature

Modeling muscle function using experimentally determined subject-specific muscle properties.

Project description:Muscle models are commonly based on intrinsic properties pooled across a number of individuals, often from a different species, and rarely validated against directly measured muscle forces. Here we use a rich data set of rat medial gastrocnemius muscle forces recorded during in-situ and in-vivo isometric, isotonic, and cyclic contractions to test the accuracy of forces predicted using Hill-type muscle models. We identified force-length and force-velocity parameters for each individual, and used either these subject-specific intrinsic properties, or population-averaged properties within the models. The modeled forces for cyclic in-vivo and in-situ contractions matched with measured muscle-tendon forces with r2 between 0.70 and 0.86, and root-mean square errors (RMSE) of 0.10 to 0.13 (values normalized to the maximum isometric force). The modeled forces were least accurate at the highest movement and cycle frequencies and did not show an improvement in r2 when subject-specific intrinsic properties were used; however, there was a reduction in the RMSE with fewer predictions having higher errors. We additionally recorded and tested muscle models specific to proximal and distal regions of the muscle and compared them to measures and models from the whole muscle belly: there was no improvement in model performance when using data from specific anatomical regions. These results show that Hill-type muscle models can yield very good performance for cyclic contractions typical of locomotion, with small reductions in errors when subject-specific intrinsic properties are used.

| S-EPMC7927351 | biostudies-literature

A mixed branch length model of heterotachy improves phylogenetic accuracy.

Project description:Evolutionary relationships are typically inferred from molecular sequence data using a statistical model of the evolutionary process. When the model accurately reflects the underlying process, probabilistic phylogenetic methods recover the correct relationships with high accuracy. There is ample evidence, however, that models commonly used today do not adequately reflect real-world evolutionary dynamics. Virtually all contemporary models assume that relatively fast-evolving sites are fast across the entire tree, whereas slower sites always evolve at relatively slower rates. Many molecular sequences, however, exhibit site-specific changes in evolutionary rates, called "heterotachy." Here we examine the accuracy of 2 phylogenetic methods for incorporating heterotachy, the mixed branch length model--which incorporates site-specific rate changes by summing likelihoods over multiple sets of branch lengths on the same tree--and the covarion model, which uses a hidden Markov process to allow sites to switch between variable and invariable as they evolve. Under a variety of simple heterogeneous simulation conditions, the mixed model was dramatically more accurate than homotachous models, which were subject to topological biases as well as biases in branch length estimates. When data were simulated with strong versions of the types of heterotachy observed in real molecular sequences, the mixed branch length model was more accurate than homotachous techniques. Analyses of empirical data sets confirmed that the mixed branch length model can improve phylogenetic accuracy under conditions that cause homotachous models to fail. In contrast, the covarion model did not improve phylogenetic accuracy compared with homotachous models and was sometimes substantially less accurate. We conclude that a mixed branch length approach, although not the solution to all phylogenetic errors, is a valuable strategy for improving the accuracy of inferred trees.

| S-EPMC3299401 | biostudies-literature

Phylogenetic tree reconstruction accuracy and model fit when proportions of variable sites change across the tree.

Project description:Commonly used phylogenetic models assume a homogeneous process through time in all parts of the tree. However, it is known that these models can be too simplistic as they do not account for nonhomogeneous lineage-specific properties. In particular, it is now widely recognized that as constraints on sequences evolve, the proportion and positions of variable sites can vary between lineages causing heterotachy. The extent to which this model misspecification affects tree reconstruction is still unknown. Here, we evaluate the effect of changes in the proportions and positions of variable sites on model fit and tree estimation. We consider 5 current models of nucleotide sequence evolution in a Bayesian Markov chain Monte Carlo framework as well as maximum parsimony (MP). We show that for a tree with 4 lineages where 2 nonsister taxa undergo a change in the proportion of variable sites tree reconstruction under the best-fitting model, which is chosen using a relative test, often results in the wrong tree. In this case, we found that an absolute test of model fit is a better predictor of tree estimation accuracy. We also found further evidence that MP is not immune to heterotachy. In addition, we show that increased sampling of taxa that have undergone a change in proportion and positions of variable sites is critical for accurate tree reconstruction.

| S-EPMC2850392 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data