Percentile Queries in Multi-Dimensional Markov Decision Processes

Markov decision processes (MDPs) with multi-dimensional weights are useful to analyze systems with multiple objectives that may be conflicting and require the analysis of trade-offs. We study the complexity of percentile queries in such MDPs and give algorithms to synthesize strategies that enforce such constraints. Given a multi-dimensional weighted MDP and a quantitative payoff function $f$, thresholds $v_i$ (one per dimension), and probability thresholds $\alpha_i$, we show how to compute a single strategy to enforce that for all dimensions $i$, the probability of outcomes $\rho$ satisfying $f_i(\rho) \geq v_i$ is at least $\alpha_i$. We consider classical quantitative payoffs from the literature (sup, inf, lim sup, lim inf, mean-payoff, truncated sum, discounted sum). Our work extends to the quantitative case the multi-objective model checking problem studied by Etessami et al. in unweighted MDPs.Comment: Extended version of CAV 2015 pape

Identification of Coevolving Residues and Coevolution Potentials Emphasizing Structure, Bond Formation and Catalytic Coordination in Protein Evolution

The structure and function of a protein is dependent on coordinated interactions between its residues. The selective pressures associated with a mutation at one site should therefore depend on the amino acid identity of interacting sites. Mutual information has previously been applied to multiple sequence alignments as a means of detecting coevolutionary interactions. Here, we introduce a refinement of the mutual information method that: 1) removes a significant, non-coevolutionary bias and 2) accounts for heteroscedasticity. Using a large, non-overlapping database of protein alignments, we demonstrate that predicted coevolving residue-pairs tend to lie in close physical proximity. We introduce coevolution potentials as a novel measure of the propensity for the 20 amino acids to pair amongst predicted coevolutionary interactions. Ionic, hydrogen, and disulfide bond-forming pairs exhibited the highest potentials. Finally, we demonstrate that pairs of catalytic residues have a significantly increased likelihood to be identified as coevolving. These correlations to distinct protein features verify the accuracy of our algorithm and are consistent with a model of coevolution in which selective pressures towards preserving residue interactions act to shape the mutational landscape of a protein by restricting the set of admissible neutral mutations

A review of the phytochemical support for the shifting defence hypothesis

Several theories have been developed to explain why invasive species are very successful and develop into pest species in their new area. The shifting defence hypothesis (SDH) argues that invasive plant species quickly evolve towards new defence levels in the invaded area because they lack their specialist herbivores but are still under attack by local (new) generalist herbivores. The SDH predicts that plants should increase their cheap, toxic defence compounds and lower their expensive digestibility reducing compounds. As a net result resources are saved that can be allocated to growth and reproduction giving these plants a competitive edge over the local plant species. We conducted a literature study to test whether toxic defence compounds in general are increased in the invaded area and if digestibility reducing compounds are lowered. We specifically studied the levels of pyrrolizidine alkaloids, a toxin which is known for its beneficial and detrimental impact against specialists and generalists, respectively. Digestibility reducers did not show a clear trend which might be due to the small number of studies and traits measured. The meta analysis showed that toxic compounds in general and pyrrolizidine alkaloid levels specifically, increased significantly in the invaded area, supporting the predictions of the SDH that a fast evolution takes place in the allocation towards defence

Homopolymer tract length dependent enrichments in functional regions of 27 eukaryotes and their novel dependence on the organism DNA (G+C)% composition

BACKGROUND: DNA homopolymer tracts, poly(dA).poly(dT) and poly(dG).poly(dC), are the simplest of simple sequence repeats. Homopolymer tracts have been systematically examined in the coding, intron and flanking regions of a limited number of eukaryotes. As the number of DNA sequences publicly available increases, the representation (over and under) of homopolymer tracts of different lengths in these regions of different genomes can be compared. RESULTS: We carried out a survey of the extent of homopolymer tract over-representation (enrichment) and over-proportional length distribution (above expected length) primarily in the single gene documents, but including some whole chromosomes of 27 eukaryotics across the (G+C)% composition range from 20 – 60%. A total of 5.2 × 10(7 )bases from 15,560 cleaned (redundancy removed) sequence documents were analyzed. Calculated frequencies of non-overlapping long homopolymer tracts were found over-represented in non-coding sequences of eukaryotes. Long poly(dA).poly(dT) tracts demonstrated an exponential increase with tract length compared to predicted frequencies. A novel negative slope was observed for all eukaryotes between their (G+C)% composition and the threshold length N where poly(dA).poly(dT) tracts exhibited over-representation and a corresponding positive slope was observed for poly(dG).poly(dC) tracts. Tract size thresholds where over-representation of tracts in different eukaryotes began to occur was between 4 – 11 bp depending upon the organism (G+C)% composition. The higher the GC%, the lower the threshold N value was for poly(dA).poly(dT) tracts, meaning that the over-representation happens at relatively lower tract length in more GC-rich surrounding sequence. We also observed a novel relationship between the highest over-representations, as well as lengths of homopolymer tracts in excess of their random occurrence expected maximum lengths. CONCLUSIONS: We discuss how our novel tract over-representation observations can be accounted for by a few models. A likely model for poly(dA).poly(dT) tract over-representation involves the known insertion into genomes of DNA synthesized from retroviral mRNAs containing 3' polyA tails. A proposed model that can account for a number of our observed results, concerns the origin of the isochore nature of eukaryotic genomes via a non-equilibrium GC% dependent mutation rate mechanism. Our data also suggest that tract lengthening via slip strand replication is not governed by a simple thermodynamic loop energy model

Correlated Mutations: A Hallmark of Phenotypic Amino Acid Substitutions

Point mutations resulting in the substitution of a single amino acid can cause severe functional consequences, but can also be completely harmless. Understanding what determines the phenotypical impact is important both for planning targeted mutation experiments in the laboratory and for analyzing naturally occurring mutations found in patients. Common wisdom suggests using the extent of evolutionary conservation of a residue or a sequence motif as an indicator of its functional importance and thus vulnerability in case of mutation. In this work, we put forward the hypothesis that in addition to conservation, co-evolution of residues in a protein influences the likelihood of a residue to be functionally important and thus associated with disease. While the basic idea of a relation between co-evolution and functional sites has been explored before, we have conducted the first systematic and comprehensive analysis of point mutations causing disease in humans with respect to correlated mutations. We included 14,211 distinct positions with known disease-causing point mutations in 1,153 human proteins in our analysis. Our data show that (1) correlated positions are significantly more likely to be disease-associated than expected by chance, and that (2) this signal cannot be explained by conservation patterns of individual sequence positions. Although correlated residues have primarily been used to predict contact sites, our data are in agreement with previous observations that (3) many such correlations do not relate to physical contacts between amino acid residues. Access to our analysis results are provided at http://webclu.bio.wzw.tum.de/~pagel/supplements/correlated-positions/

A Coevolutionary Residue Network at the Site of a Functionally Important Conformational Change in a Phosphohexomutase Enzyme Family

Coevolution analyses identify residues that co-vary with each other during evolution, revealing sequence relationships unobservable from traditional multiple sequence alignments. Here we describe a coevolutionary analysis of phosphomannomutase/phosphoglucomutase (PMM/PGM), a widespread and diverse enzyme family involved in carbohydrate biosynthesis. Mutual information and graph theory were utilized to identify a network of highly connected residues with high significance. An examination of the most tightly connected regions of the coevolutionary network reveals that most of the involved residues are localized near an interdomain interface of this enzyme, known to be the site of a functionally important conformational change. The roles of four interface residues found in this network were examined via site-directed mutagenesis and kinetic characterization. For three of these residues, mutation to alanine reduces enzyme specificity to ∼10% or less of wild-type, while the other has ∼45% activity of wild-type enzyme. An additional mutant of an interface residue that is not densely connected in the coevolutionary network was also characterized, and shows no change in activity relative to wild-type enzyme. The results of these studies are interpreted in the context of structural and functional data on PMM/PGM. Together, they demonstrate that a network of coevolving residues links the highly conserved active site with the interdomain conformational change necessary for the multi-step catalytic reaction. This work adds to our understanding of the functional roles of coevolving residue networks, and has implications for the definition of catalytically important residues

HacA-Independent Functions of the ER Stress Sensor IreA Synergize with the Canonical UPR to Influence Virulence Traits in Aspergillus fumigatus

Endoplasmic reticulum (ER) stress is a condition in which the protein folding capacity of the ER becomes overwhelmed by an increased demand for secretion or by exposure to compounds that disrupt ER homeostasis. In yeast and other fungi, the accumulation of unfolded proteins is detected by the ER-transmembrane sensor IreA/Ire1, which responds by cleaving an intron from the downstream cytoplasmic mRNA HacA/Hac1, allowing for the translation of a transcription factor that coordinates a series of adaptive responses that are collectively known as the unfolded protein response (UPR). Here, we examined the contribution of IreA to growth and virulence in the human fungal pathogen Aspergillus fumigatus. Gene expression profiling revealed that A. fumigatus IreA signals predominantly through the canonical IreA-HacA pathway under conditions of severe ER stress. However, in the absence of ER stress IreA controls dual signaling circuits that are both HacA-dependent and HacA-independent. We found that a ΔireA mutant was avirulent in a mouse model of invasive aspergillosis, which contrasts the partial virulence of a ΔhacA mutant, suggesting that IreA contributes to pathogenesis independently of HacA. In support of this conclusion, we found that the ΔireA mutant had more severe defects in the expression of multiple virulence-related traits relative to ΔhacA, including reduced thermotolerance, decreased nutritional versatility, impaired growth under hypoxia, altered cell wall and membrane composition, and increased susceptibility to azole antifungals. In addition, full or partial virulence could be restored to the ΔireA mutant by complementation with either the induced form of the hacA mRNA, hacAi, or an ireA deletion mutant that was incapable of processing the hacA mRNA, ireAΔ10. Together, these findings demonstrate that IreA has both HacA-dependent and HacA-independent functions that contribute to the expression of traits that are essential for virulence in A. fumigatus

A Role for the Unfolded Protein Response (UPR) in Virulence and Antifungal Susceptibility in Aspergillus fumigatus

Filamentous fungi rely heavily on the secretory pathway, both for the delivery of cell wall components to the hyphal tip and the production and secretion of extracellular hydrolytic enzymes needed to support growth on polymeric substrates. Increased demand on the secretory system exerts stress on the endoplasmic reticulum (ER), which is countered by the activation of a coordinated stress response pathway termed the unfolded protein response (UPR). To determine the contribution of the UPR to the growth and virulence of the filamentous fungal pathogen Aspergillus fumigatus, we disrupted the hacA gene, encoding the major transcriptional regulator of the UPR. The ΔhacA mutant was unable to activate the UPR in response to ER stress and was hypersensitive to agents that disrupt ER homeostasis or the cell wall. Failure to induce the UPR did not affect radial growth on rich medium at 37°C, but cell wall integrity was disrupted at 45°C, resulting in a dramatic loss in viability. The ΔhacA mutant displayed a reduced capacity for protease secretion and was growth-impaired when challenged to assimilate nutrients from complex substrates. In addition, the ΔhacA mutant exhibited increased susceptibility to current antifungal agents that disrupt the membrane or cell wall and had attenuated virulence in multiple mouse models of invasive aspergillosis. These results demonstrate the importance of ER homeostasis to the growth and virulence of A. fumigatus and suggest that targeting the UPR, either alone or in combination with other antifungal drugs, would be an effective antifungal strategy