Proteome-Wide Structural Computations Provide Insights into Empirical Amino Acid Substitution Matrices.

The relative contribution of mutation and selection to the amino acid substitution rates observed in empirical matrices is unclear. Herein, we present a neutral continuous fitness-stability model, inspired by the Arrhenius law (qij=aije-ΔΔGij). The model postulates that the rate of amino acid substitution (i→j) is determined by the product of a pre-exponential factor, which is influenced by the genetic code structure, and an exponential term reflecting the relative fitness of the amino acid substitutions. To assess the validity of our model, we computed changes in stability of 14,094 proteins, for which 137,073,638 in silico mutants were analyzed. These site-specific data were summarized into a 20 square matrix, whose entries, ΔΔGij, were obtained after averaging through all the sites in all the proteins. We found a significant positive correlation between these energy values and the disease-causing potential of each substitution, suggesting that the exponential term accurately summarizes the fitness effect. A remarkable observation was that amino acids that were highly destabilizing when acting as the source, tended to have little effect when acting as the destination, and vice versa (source → destination). The Arrhenius model accurately reproduced the pattern of substitution rates collected in the empirical matrices, suggesting a relevant role for the genetic code structure and a tuning role for purifying selection exerted via protein stability.

Pine has two glutamine synthetase paralogs, GS1b.1 and GS1b.2, exhibiting distinct biochemical properties.

The enzyme glutamine synthetase (EC 6.3.1.2) is mainly responsible for the incorporation of inorganic nitrogen into organic molecules in plants. In the present work, a pine (Pinus pinaster) GS1 (PpGS1b.2) gene was identified, showing a high sequence identity with the GS1b.1 gene previously characterized in conifers. Phylogenetic analysis revealed that the presence of PpGS1b.2 is restricted to the genera Pinus and Picea and is not found in other conifers. Gene expression data suggest a putative role of PpGS1b.2 in plant development, similar to other GS1b genes from angiosperms, suggesting evolutionary convergence. The characterization of GS1b.1 and GS1b.2 at the structural, physicochemical, and kinetic levels has shown differences even though they have high sequence homology. GS1b.2 had a lower optimum pH (6 vs. 6.5) and was less thermally stable than GS1b.1. GS1b.2 exhibited positive cooperativity for glutamate and substrate inhibition for ammonium. However, GS1b.1 exhibited substrate inhibition behavior for glutamate and ATP. Alterations in the kinetic characteristics produced by site-directed mutagenesis carried out in this work strongly suggest an implication of amino acids at positions 264 and 267 in the active center of pine GS1b.1 and GS1b.2 being involved in affinity toward ammonium. Therefore, the amino acid differences between GS1b.1 and GS1b.2 would support the functioning of both enzymes to meet distinct plant needs.

A Census of Human Methionine-Rich Prion-like Domain-Containing Proteins.

Methionine-rich prion-like proteins can regulate liquid-liquid phase separation processes in response to stresses. To date, however, very few proteins have been identified as methionine-rich prion-like. Herein, we have performed a computational survey of the human proteome to search for methionine-rich prion-like domains. We present a census of 51 manually curated methionine-rich prion-like proteins. Our results show that these proteins tend to be modular in nature, with molecular sizes significantly greater than those we would expect due to random sampling effects. These proteins also exhibit a remarkably high degree of spatial compaction when compared to average human proteins, even when protein size is accounted for. Computational evidence suggests that such a high degree of compactness might be due to the aggregation of methionine residues, pointing to a potential redox regulation of compactness. Gene ontology and network analyses, performed to shed light on the biological processes in which these proteins might participate, indicate that methionine-rich and non-methionine-rich prion-like proteins share gene ontology terms related to the regulation of transcription and translation but, more interestingly, these analyses also reveal that proteins from the methionine-rich group tend to share more gene ontology terms among them than they do with their non-methionine-rich prion-like counterparts.

Phylogenies from unaligned proteomes using sequence environments of amino acid residues.

Alignment-free methods for sequence comparison and phylogeny inference have attracted a great deal of attention in recent years. Several algorithms have been implemented in diverse software packages. Despite the great number of existing methods, most of them are based on word statistics. Although they propose different filtering and weighting strategies and explore different metrics, their performance may be limited by the phylogenetic signal preserved in these words. Herein, we present a different approach based on the species-specific amino acid neighborhood preferences. These differential preferences can be assessed in the context of vector spaces. In this way, a distance-based method to build phylogenies has been developed and implemented into an easy-to-use R package. Tests run on real-world datasets show that this method can reconstruct phylogenetic relationships with high accuracy, and often outperforms other alignment-free approaches. Furthermore, we present evidence that the new method can perform reliably on datasets formed by non-orthologous protein sequences, that is, the method not only does not require the identification of orthologous proteins, but also does not require their presence in the analyzed dataset. These results suggest that the neighborhood preference of amino acids conveys a phylogenetic signal that may be of great utility in phylogenomics.

renz: An R package for the analysis of enzyme kinetic data.

BACKGROUND: Complex enzymatic models are required for analyzing kinetic data derived under conditions that may not satisfy the assumptions associated with Michaelis-Menten kinetics. To analyze these data, several software packages have been developed. However, the complexity introduced by these programs is often dispensable when analyzing data conforming to the canonical Michaelis-Menten model. In these cases, the sophisticated routines of these packages become inefficient and unnecessarily intricated for the intended purpose, reason for which most users resort to general-purpose graphing programs. However, this approach, in addition of being time-consuming, is prone to human error, and can lead to misleading estimates of kinetic parameters, particularly when unweighted regression analyses of transformed kinetic data are performed. RESULTS: To fill the existing gap between highly specialized and general-purpose software, we have developed an easy-to-use R package, renz, designed for accurate and efficient estimation of enzyme kinetic parameters. The package provides different methods that can be clustered into four categories, depending on whether they are based on data fitting to a single progress curve (evolution of substrate concentration over time) or, alternatively, based on the dependency of initial rates on substrate concentration (differential rate equation). A second criterion to be considered is whether the experimental data need to be manipulated to obtain linear functions or, alternatively, data are directly fitted using non-linear regression analysis. The current program is a cross-platform, free and open-source software that can be obtained from the CRAN repository. The package is accompanied by five vignettes, which are intended to guide users to choose the appropriate method in each case, as well as providing the basic theoretical foundations of each method. These vignettes use real experimental data to illustrate the use of the package utilities. CONCLUSIONS: renz is a rigorous and yet easy-to-use software devoted to the analysis of kinetic data. This application has been designed to meet the needs of users who are not practicing enzymologists, but who need to accurately estimate the kinetic parameters of enzymes. The current software saves time and minimizes the risk of making mistakes or introducing biases due to uncorrected error propagation effects.

The Role of Methionine Residues in the Regulation of Liquid-Liquid Phase Separation.

Membraneless organelles are non-stoichiometric supramolecular structures in the micron scale. These structures can be quickly assembled/disassembled in a regulated fashion in response to specific stimuli. Membraneless organelles contribute to the spatiotemporal compartmentalization of the cell, and they are involved in diverse cellular processes often, but not exclusively, related to RNA metabolism. Liquid-liquid phase separation, a reversible event involving demixing into two distinct liquid phases, provides a physical framework to gain insights concerning the molecular forces underlying the process and how they can be tuned according to the cellular needs. Proteins able to undergo phase separation usually present a modular architecture, which favors a multivalency-driven demixing. We discuss the role of low complexity regions in establishing networks of intra- and intermolecular interactions that collectively control the phase regime. Post-translational modifications of the residues present in these domains provide a convenient strategy to reshape the residue-residue interaction networks that determine the dynamics of phase separation. Focus will be placed on those proteins with low complexity domains exhibiting a biased composition towards the amino acid methionine and the prominent role that reversible methionine sulfoxidation plays in the assembly/disassembly of biomolecular condensates.

ptm: an R package for the study of methionine sulfoxidation and other posttranslational modifications.

SUMMARY: Methionine sulfoxidation is a posttranslational modification (PTM) playing important roles in cell signaling. Herein, we present ptm, an R package for the study of this modification. However, since many of the analyses applied to methionine modification can be extended to other modifications, the package can be useful to thoroughly analyze PTMs in general. Thus, within a single software environment, ptm can integrate information from up to 11 databases covering nine modifications. Different functions can work coordinately to form pipelines allowing the programmatic analysis of thousands of proteins. Alternatively, the user can simultaneously perform different analyses on the same protein of interest, combining the results into a single output. The flexibility of ptm makes it a suitable tool to address site- and protein-centric hypotheses related to PTMs. Accompanying the package, we maintain a web page containing extended documentation and examples of the tasks that can be performed with ptm. AVAILABILITY AND IMPLEMENTATION: ptm is implemented in R. Release versions are available via CRAN and work on all major operating systems. The development version is maintained at https://bitbucket.org/jcaledo/ptm. Extended documentation can be found at https://metositeptm.com. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MetOSite: an integrated resource for the study of methionine residues sulfoxidation.

MOTIVATION: The oxidation of protein-bound methionine to form methionine sulfoxide has traditionally been regarded as an oxidative damage. However, growing evidences support the view of this reversible reaction also as a regulatory post-translational modification. Thus, the oxidation of methionine residues has been reported to have multiple and varied implications for protein function. However, despite the importance of this modification and the abundance of reports, all these data are scattered in the literature. No database/resource on methionine sulfoxidation exists currently. Since this information is useful to gain further insights into the redox regulation of cellular proteins, we have created a primary database of experimentally confirmed sulfoxidation sites. RESULTS: MetOSite currently contains 7242 methionine sulfoxide sites found in 3562 different proteins from 23 species, with Homo sapiens, Arabidopsis thaliana and Bacillus cereus as the main contributors. Each collected site has been classified according to the effect of its sulfoxidation on the biological properties of the modified protein. Thus, MetOSite documents cases where the sulfoxidation of methionine leads to (i) gain of activity, (ii) loss of activity, (iii) increased protein-protein interaction susceptibility, (iv) decreased protein-protein interaction susceptibility, (v) changes in protein stability and (vi) changes in subcellular location. AVAILABILITY AND IMPLEMENTATION: MetOSite is available at https://metosite.uma.es.

Multisite phosphorylation provides a reliable mechanism for making decisions in noisy environments.

The ability to make decisions at the cellular level is absolutely critical for the survival of organisms. Eukaryotic cells are constantly making binary decisions in response to internal and environmental signals. Among the most notable transducers of information are protein kinases. The regulation of these signaling proteins often relies on the activity of other protein kinases located upstream in the signaling cascade. However, these signaling systems are by their own nature an important source of molecular noise. Herein, we have assessed the role of multisite phosphorylation in detecting signals in the face of molecular noise. To address this issue, we have conceptually envisioned the biochemical transduction machinery as a classifier model that can lead to four possible outputs: true positives and negatives, and false positives and negatives. In this probabilistic framework, we show that multisite phosphorylation represents a mechanism to filter noise during the decision-making process. We present results showing that nonessential phosphorylation sites contribute to increase the rate of true positives while, at the same time, they can lessen the rate of false positives. This simultaneous increase in sensitivity and specificity, makes multisite phosphorylation a valuable and easily implemented mechanism to reliably transduce information in noisy contexts.

Inferring Methionine Sulfoxidation and serine Phosphorylation crosstalk from Phylogenetic analyses.

BACKGROUND: The sulfoxidation of methionine residues within the phosphorylation motif of protein kinase substrates, may provide a mechanism to couple oxidative signals to changes in protein phosphorylation. Herein, we hypothesize that if the residues within a pair of phosphorylatable-sulfoxidable sites are functionally linked, then they might have been coevolving. To test this hypothesis a number of site pairs previously detected on human stress-related proteins has been subjected to analysis using eukaryote ortholog sequences and a phylogenetic approach. RESULTS: Overall, the results support the conclusion that in the eIF2α protein, serine phosphorylation at position 218 and methionine oxidation at position 222, belong to the same functional network. First, the observed data were much better fitted by Markovian models that assumed coevolution of both sites, with respect to their counterparts assuming independent evolution (p-value = 0.003). Second, this conclusion was robust with respect to the methods used to reconstruct the phylogenetic relationship between the 233 eukaryotic species analyzed. Third, the co-distribution of phosphorylatable and sulfoxidable residues at these positions showed multiple origins throughout the evolution of eukaryotes, which further supports the view of an adaptive value for this co-occurrence. Fourth, the possibility that the coevolution of these two sites might be due to structure-driven compensatory mutations was evaluated. The results suggested that factors other than those merely structural were behind the observed coevolution. Finally, the relationship detected between other modifiable site pairs from ataxin-2 (S814-M815), ataxin-2-like (S211-M215) and Pumilio homolog 1 (S124-M125), reinforce the view of a role for phosphorylation-sulfoxidation crosstalk. CONCLUSIONS: For the four stress-related proteins analyzed herein, their respective pairs of PTM sites (phosphorylatable serine and sulfoxidable methionine) were found to be evolving in a correlated fashion, which suggests a relevant role for methionine sulfoxidation and serine phosphorylation crosstalk in the control of protein translation under stress conditions.

Aging is neither a failure nor an achievement of natural selection.

In contraposition to the view of aging as a stochastic time-dependent accumulation of damage, phenoptotic theories of aging postulate that senescence may provide supra-individual advantages, and therefore it might have been promoted by natural selection. We reason that although programmed aging theories are subjectively appealing because they convey a cure for aging, they also raise a number of objections that need to be dealt with, before we may be entitled to contemplate aging as an adaptive function evolved through natural selection. As an alternative view, we present metabolism as an endless source of by-products and errors causing cellular damage. Although this damage accumulation event is a spontaneous entropy-driven process, its kinetics can be genetically and environmentally modulated, giving place to the wide range of lifespans we observe. Mild forms of damage may be accumulating during a long enough period of time to allow reproduction before the fatal failure happens. Hence, aging would be a stochastic process out of the reach of natural selection. However, those genetic pathways influencing the rate of aging and consequently determining longevity may be targets of natural selection and may contribute to shaping the optimal strategy according to the ecological context. In this sense, short- and long-lived organisms represent two extreme strategies that, in terms of biological fitness, can perform equally well, each within its own niche.

Aging is neither a Failure nor an Achievement of Natural Selection.

In contraposition to the view of aging as a stochastic time-dependent accumulation of damage, phenoptotic theories of aging postulate that senescence may provide supra-individual advantages, and therefore it might have been promoted by natural selection. We reason that although programmed aging theories are subjectively appealing because they convey a cure for aging, they also raise a number of objections that need to be dealt with, before we may be entitled to contemplate aging as an adaptive function evolved through natural selection. As an alternative view, we present metabolism as an endless source of by-products and errors causing cellular damage. Although this damage accumulation event is a spontaneous entropy-driven process, its kinetics can be genetically and environmentally modulated, giving place to the wide range of lifespans we observe. Mild forms of damage may be accumulating during a long enough period of time to allow reproduction before the fatal failure happens. Hence, aging would be a stochastic process out of the reach of natural selection. However, those genetic pathways influencing the rate of aging and consequently determining longevity may be targets of natural selection and may contribute to shaping the optimal strategy according to the ecological context. In this sense, short- and long-lived organisms represent two extreme strategies that, in terms of biological fitness, can perform equally well, each within its own niche.

Protein-protein interfaces from cytochrome c oxidase I evolve faster than nonbinding surfaces, yet negative selection is the driving force.

Respiratory complexes are encoded by two genomes (mitochondrial DNA [mtDNA] and nuclear DNA [nDNA]). Although the importance of intergenomic coadaptation is acknowledged, the forces and constraints shaping such coevolution are largely unknown. Previous works using cytochrome c oxidase (COX) as a model enzyme have led to the so-called "optimizing interaction" hypothesis. According to this view, mtDNA-encoded residues close to nDNA-encoded residues evolve faster than the rest of positions, favoring the optimization of protein-protein interfaces. Herein, using evolutionary data in combination with structural information of COX, we show that failing to discern the effects of interaction from other structural and functional effects can lead to deceptive conclusions such as the "optimizing hypothesis." Once spurious factors have been accounted for, data analysis shows that mtDNA-encoded residues engaged in contacts are, in general, more constrained than their noncontact counterparts. Nevertheless, noncontact residues from the surface of COX I subunit are a remarkable exception, being subjected to an exceptionally high purifying selection that may be related to the maintenance of a suitable heme environment. We also report that mtDNA-encoded residues involved in contacts with other mtDNA-encoded subunits are more constrained than mtDNA-encoded residues interacting with nDNA-encoded polypeptides. This differential behavior cannot be explained on the basis of predicted thermodynamic stability, as interactions between mtDNA-encoded subunits contribute more weakly to the complex stability than those interactions between subunits encoded by different genomes. Therefore, the higher conservation observed among mtDNA-encoded residues involved in intragenome interactions is likely due to factors other than structural stability.

Life-history Constraints on the Mechanisms that Control the Rate of ROS Production.

The quest to understand why and how we age has led to numerous lines of investigation that have gradually converged to consider mitochondrial metabolism as a major player. During mitochondrial respiration a small and variable amount of the consumed oxygen is converted to reactive species of oxygen (ROS). For many years, these ROS have been perceived as harmful by-products of respiration. However, evidence from recent years indicates that ROS fulfill important roles as cellular messengers. Results obtained using model organisms suggest that ROS-dependent signalling may even activate beneficial cellular stress responses, which eventually may lead to increased lifespan. Nevertheless, when an overload of ROS cannot be properly disposed of, its accumulation generates oxidative stress, which plays a major part in the ageing process. Comparative studies about the rates of ROS production and oxidative damage accumulation, have led to the idea that the lower rate of mitochondrial oxygen radical generation of long-lived animals with respect to that of their short-lived counterpart, could be a primary cause of their slow ageing rate. A hitherto largely under-appreciated alternative view is that such lower rate of ROS production, rather than a cause may be a consequence of the metabolic constraints imposed for the large body sizes that accompany high lifespans. To help understanding the logical underpinning of this rather heterodox view, herein I review the current literature regarding the mechanisms of ROS formation, with particular emphasis on evolutionary aspects.

On the relationships of postcanine tooth size with dietary quality and brain volume in primates: implications for hominin evolution.

Brain volume and cheek-tooth size have traditionally been considered as two traits that show opposite evolutionary trends during the evolution of Homo. As a result, differences in encephalization and molarization among hominins tend to be interpreted in paleobiological grounds, because both traits were presumably linked to the dietary quality of extinct species. Here we show that there is an essential difference between the genus Homo and the living primate species, because postcanine tooth size and brain volume are related to negative allometry in primates and show an inverse relationship in Homo. However, when size effects are removed, the negative relationship between encephalization and molarization holds only for platyrrhines and the genus Homo. In addition, there is no general trend for the relationship between postcanine tooth size and dietary quality among the living primates. If size and phylogeny effects are both removed, this relationship vanishes in many taxonomic groups. As a result, the suggestion that the presence of well-developed postcanine teeth in extinct hominins should be indicative of a poor-quality diet cannot be generalized to all extant and extinct primates.

Plasma long-chain free fatty acids predict mammalian longevity.

Membrane lipid composition is an important correlate of the rate of aging of animals and, therefore, the determination of their longevity. In the present work, the use of high-throughput technologies allowed us to determine the plasma lipidomic profile of 11 mammalian species ranging in maximum longevity from 3.5 to 120 years. The non-targeted approach revealed a specie-specific lipidomic profile that accurately predicts the animal longevity. The regression analysis between lipid species and longevity demonstrated that the longer the longevity of a species, the lower is its plasma long-chain free fatty acid (LC-FFA) concentrations, peroxidizability index, and lipid peroxidation-derived products content. The inverse association between longevity and LC-FFA persisted after correction for body mass and phylogenetic interdependence. These results indicate that the lipidomic signature is an optimized feature associated with animal longevity, emerging LC-FFA as a potential biomarker of longevity.

Mutational bias plays an important role in shaping longevity-related amino acid content in mammalian mtDNA-encoded proteins.

During the course of evolution, amino acid shifts might have resulted in mitochondrial proteomes better endowed to resist oxidative stress. However, owing to the problem of distinguishing between functional constraints/adaptations in protein sequences and mutation-driven biases in the composition of these sequences, the adaptive value of such amino acid shifts remains under discussion. We have analyzed the coding sequences of mtDNA from 173 mammalian species, dissecting the effect of nucleotide composition on amino acid usages. We found remarkable cysteine avoidance in mtDNA-encoded proteins. However, no effect of longevity on cysteine content could be detected. On the other hand, nucleotide compositional shifts fully accounted for threonine usages. In spite of a strong effect of mutational bias on methionine abundances, our results suggest a role of selection in determining the composition of methionine. Whether this selective effect is linked or not to protection against oxidative stress is still a subject of debate.

Thermodynamic stability explains the differential evolutionary dynamics of cytochrome b and COX I in mammals.

By using a combination of evolutionary and structural data from 231 species, we have addressed the relationship between evolution and structural features of cytochrome b and COX I, two mtDNA-encoded proteins. The interior of cytochrome b, in contrast to that of COX I, exhibits a remarkable tolerance to changes. The higher evolvability of cytochrome b contrasts with the lower rate of synonymous substitutions of its gene when compared to that of COX I, suggesting that the latter is subjected to a stronger purifying selection. We present evidences that the stability effect of mutations (ΔΔG) may be behind these differential behaviour.

SNAT2 transceptor signalling via mTOR: a role in cell growth and proliferation?

We have investigated the effect of chronic competitive inhibition of SNAT2 (System A) amino acid (AA) transport, induced by incubation with a saturating dose of a non-metabolisable System A amino acid analogue (Me-AIB), on growth and proliferation of MCF-7 human breast cancer cells in complete culture medium. These cells express Na+- and pH-dependent SNAT2 AA transport and a saturating concentration of Me-AIB (10 mM) competitively inhibits (>90%) AA uptake via SNAT2. Incubation with Me-AIB for up to 5 days progressively reduced cell proliferation (~2-fold) and depleted intracellular concentrations of not only SNAT2 AA substrates but of essential branched chain AAs (e.g. leucine). Surprisingly, total cellular protein was maintained and cells subjected to chronic Me-AIB incubation exhibited a detectable increase in cell size. Analysis of mTOR signalling revealed that, despite a substantial reduction in size of the intracellular AA pool, Me-AIB elevated mTOR-dependent p70S6K1 phosphorylation. Proteomic analysis of TAP-tag purified SNAT2 fusion proteins identified two novel SNAT2-interacting proteins that may potentially function in conjunction with the SNAT2 transceptor to regulate signalling pathways influencing protein turnover and cell growth.

Mitochondrially encoded methionine is inversely related to longevity in mammals.

Methionine residues in proteins react readily with reactive oxygen species making them particularly sensitive to oxidation. However, because oxidized methionine can be reduced back in a catalyzed reaction, it has been suggested that methionine residues act as oxidant scavengers, protecting not only the proteins where they are located but also the surrounding macromolecules. To investigate whether methionine residues may be selected for or against animal longevity, we carried out a meta-examination of mitochondrial genomes from mammalian species. Our analyses unveiled a hitherto unnoticed observation: mitochondrially encoded polypeptides from short-lived species are enriched in methionine when compared with their long-lived counterparts. We show evidence suggesting that methionine addition to proteins in short-lived species, rather than methionine loss from proteins in long-lived species, is behind the reported difference in methionine usage. The inverse association between longevity and methionine, which persisted after correction for body mass and phylogenetic interdependence, was paralleled by the methionine codon AUA, but not by the codon AUG. Although nuclear encoded mitochondrial polypeptides exhibited higher methionine usage than nonmitochondrial proteins, correlation with longevity was only found within the group of those polypeptides located in the inner mitochondrial membrane. Based on these results, we propose that short-lived animals subjected to higher oxidative stress selectively accumulate methionine in their mitochondrially encoded proteins, which supports the role of oxidative damage in aging.