Next Article in Journal
Cellular Pre-Adaptation to the High O2 Concentration Used in Standard Cell Culture Confers Resistance to Subsequent H2O2-Induced Cell Death
Previous Article in Journal
Crosstalk between Oxidative Stress and Inflammation Caused by Noise and Air Pollution—Implications for Neurodegenerative Diseases
Previous Article in Special Issue
Mechanism of Antiradical Activity of Coumarin-Trihydroxybenzohydrazide Derivatives: A Comprehensive Kinetic DFT Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cysteine Is the Only Universally Affected and Disfavored Proteomic Amino Acid under Oxidative Conditions in Animals

by
Mario Schindeldecker
1,2 and
Bernd Moosmann
1,3,*
1
Evolutionary Biochemistry and Redox Medicine, Institute for Pathobiochemistry, University Medical Center of the Johannes Gutenberg University, 55128 Mainz, Germany
2
Institute for Pathology, University Medical Center of the Johannes Gutenberg University, 55131 Mainz, Germany
3
Institute for Quantitative and Computational Biosciences, Johannes Gutenberg University, 55128 Mainz, Germany
*
Author to whom correspondence should be addressed.
Antioxidants 2024, 13(3), 267; https://doi.org/10.3390/antiox13030267
Submission received: 31 December 2023 / Revised: 7 February 2024 / Accepted: 19 February 2024 / Published: 22 February 2024
(This article belongs to the Special Issue Oxidative Stress and Antioxidants in Computational Chemistry)

Abstract

:
Oxidative modifications of amino acid side chains in proteins are a hallmark of oxidative stress, and they are usually regarded as structural damage. However, amino acid oxidation may also have a protective effect and may serve regulatory or structural purposes. Here, we have attempted to characterize the global redox role of the 20 proteinogenic amino acids in animals by analyzing their usage frequency in 5 plausible evolutionary paradigms of increased oxidative burden: (i) peroxisomal proteins versus all proteins, (ii) mitochondrial proteins versus all proteins, (iii) mitochondrially encoded respiratory chain proteins versus all mitochondrial proteins, (iv) proteins from long-lived animals versus those from short-lived animals, and (v) proteins from aerobic, free-living animals versus those from facultatively anaerobic animals. We have found that avoidance of cysteine in the oxidative condition was the most pronounced and significant variation in the majority of comparisons. Beyond this preeminent pattern, only local signals were observed, primarily increases in methionine and glutamine as well as decreases in serine and proline. Hence, certain types of cysteine oxidation appear to enforce its proteome-wide evolutionary avoidance despite its essential role in disulfide bond formation and metal ligation. The susceptibility to oxidation of all other amino acids appears to be generally unproblematic, and sometimes advantageous.

1. Introduction

Chemical modifications of biological macromolecules are important challenges to cellular and organismal fitness that are involved in numerous diseases and in the biological aging process [1,2,3]. In general, biological macromolecular persistence may be achieved by different strategies, among them reduced numbers or levels of offenders to structural integrity (such as oxidants), intrinsically enhanced stability (due to altered usage of structural building blocks), or improved mechanisms of replacement and repair (such as reducing enzymes, chaperones, etc.). Regarding the structural stability of proteins to chemical oxidation [4], two opposing effects need to be separated: amino acid oxidation with overall adverse consequences, and amino acid oxidation with overall beneficial consequences. The former may arise when a major or minor route of oxidation leads to irreparable damage of the protein or its surroundings, such as during protein cysteine thiyl radical formation [5,6], which appears to account for the avoidance of mitochondrial membrane cysteines in long-lived animal species [7,8,9]. Cysteine thiyl radicals can relay damage into the interior of proteins [10] and act as chain-transfer catalysts in vivo [5]; they are incompletely scavenged by aqueous antioxidants like ascorbate and glutathione [5,6,11], and no efficient system seems to exist for their scavenging in hydrophobic environments [6,12]. Thus, cysteine one-electron oxidation or a related type of oxidation may account at least in part for the very low degree of conservation of solitary surface cysteine residues, as opposed to paired cysteine residues, in the general proteome [13,14], since disulfides are much more resistant to one-electron oxidation and less damaging in the one-electron oxidized state than free thiol groups [6,11,15].
In contrast, other forms of amino acid oxidation may have systemically beneficial consequences, especially when the individually effectuated chemical modification is structurally rather harmless and, preferentially, readily repairable. Under these conditions, the oxidative event may have prevented potentially much more dangerous chemical damage elsewhere in the cell. For example, protein methionine oxidation to methionine sulfoxide can be repaired enzymatically [16,17] and induces only moderate structural deterioration [16,18,19], but effectively scavenges reactive species that might have triggered a more deleterious reaction cascade elsewhere. For example, respiratory chain complexes have been shown to profoundly accumulate this labile amino acid in parallel with metabolic rate and reactive oxygen species (ROS) production, clearly indicative of an adaptive antioxidant mechanism [18,20,21,22].
In the present study, we have attempted to clarify in a systematic fashion whether the two well-established examples above, mitochondrial cysteine decumulation and mitochondrial methionine accumulation, indeed constituted the major proteomic amino acid adaptations towards oxidative stress in animals, or whether they just represented showcase examples of a much wider and diverse body of potentially adaptive, quantitative changes. Hence, a testing battery of five evolutionary paradigms of persistently increased oxidative burden in animals was devised, and the corresponding proteomes were comparatively evaluated for differences in amino acid usages. The five paradigms of comparison were designated “Peroxisome” (I), “Mitochondrion” (II), “Respiratory Chain” (III), “Longevity” (IV), and “Aerobicity” (V), and they relate to oxidative stress as follows: (i) Peroxisomes have long been recognized for their high levels of hydrogen peroxide production, which locally and sometimes globally involves substantial oxidative stress despite high local concentrations of the antioxidant enzyme, catalase [23,24]. (ii) Mitochondria usually dominate cellular oxidant production due to their potentially high superoxide release from flavoenzymes and ubiquinone oxidoreductases [25,26,27]. (iii) The prime source of superoxide within mitochondria are the respiratory chain complexes of the inner mitochondrial membrane [25,26], which are by themselves in need of special antioxidant protection [18,28]. (iv) Longevity requires particular antioxidant long-term strategies of prevention and repair [1,29,30], several of which are realized through structural amino acid usage changes [7,9,12,20,31,32,33]. (v) As opposed to anaerobic fermentation, aerobic respiration involves cellular oxygen toxicity, which is mediated by the production of reactive oxygen species (ROS) [34,35]. Through their regulatory shut-down of respiration, at least during adulthood, facultatively anaerobic parasitic helminths largely avoid oxidative stress despite the continued histological presence of mitochondria [7,36,37]. Finally, to put the obtained animal data into a wider perspective, a collection of aerobic versus anaerobic prokaryotes with well-defined oxygen requirements was sampled and analyzed accordingly.

2. Materials and Methods

2.1. Inclusion Criteria for the Analyzed Species and Proteins

A basis set for 20 animal species with well-assigned proteomes was defined from an earlier working set of 28 species [7] that had been assembled based on the following criteria: (i) an annotated nuclear-encoded proteome of >105 amino acids in the NCBI RefSeq database (ncbi.nlm.nih.gov/refseq), (ii) availability of reliable lifespan data from an authoritative resource book [38] or from the peer-reviewed zoological literature, and (iii) consideration of only the most comprehensively annotated species per genus to avoid the overrepresentation of single genera. For the current study, the additional criterion of >102 annotated mitochondrially imported protein sequences in NCBI RefSeq was adopted. Partial sequences and isoforms were generally excluded during proteome assembly; however, no other restrictions were applied. Peroxisomal proteins for each species were obtained from the curated Uniprot database (uniprot.org). The resulting species list and the evaluated numbers of proteins are provided in Table 1.
An extended animal species set of 218 species (chordates and arthropods) with sequenced and annotated mitochondrial genomes was compiled as described [7]. Lifespan data of these species were predominantly taken from the same resource as before [38], with the exception of the amphibians, for which another reference work was consulted amino acid usage alterations in five p [39]. Less than 10% of the lifespans were retrieved from varying sources (scholarly articles and textbooks) as detailed in [7].
A list of 86 helminths (predominantly nematodes, platyhelminths and annelids) with sequenced and annotated mitochondrial genomes was assembled as described [8]. These invertebrate species were classified according to their free-living versus parasitic lifestyle as adults; permanent adult parasitism was adopted as the criterion for anaerobic fermentative energy generation, whereas animals with a free-living lifestyle, species feeding on plants, and ectoparasitic species reversibly feeding on alternating animal hosts were classified as aerobic [8,36,37].
Prokaryotic proteomes of the groups archaea, Gram-positive bacteria and Gram-negative bacteria were assembled according to criteria already adopted before [7], especially the characterization of the prokaryotes’ growth conditions (as aerobic or anaerobic) in the authoritative database of the “Deutsche Sammlung von Mikroorganismen und Zellkulturen” (dsmz.de/collection/catalogue/microorganisms/catalogue). Strictly aerobic, microaerophilic, and facultatively anaerobic prokaryotes were collectively classified as ‘aerobic’ because of their fundamental aerotolerance. Only one species or strain per genus was included for analysis to avoid consideration bias. For the current study, sequences were obtained from the NCBI RefSeq database. As regards the archaea, the original basis set was expanded by all further species listed in a reference resource [40], resulting in 16 anaerobes and 11 aerobes. Regarding the Gram-positive bacteria, the original basis set was expanded by additional species listed in a clinical microbiology textbook [41], yielding 18 anaerobes and 16 aerobes. Regarding the Gram-negative bacteria, the original basis set was carried over, which rendered 15 anaerobes and 33 aerobes after consideration of 2 recent species reclassifications in the NCBI RefSeq database.

2.2. Data Analysis

Amino acid frequency data were either analyzed by comparison of group averages, yielding effect size ratios (with ratio = 1 indicating equity and, thus, unaltered usage of an amino acid under two conditions), or by linear correlation (for various analyses of the effects of longevity). Linear correlation analysis was applied to amino acid usage ratios as determined above, calculating log-linear correlation coefficients of all 20 individual amino acid usage ratios from paradigms I, II, and III with log lifespan. Moreover, linear correlation analysis was also applied to direct amino acid frequencies, calculating log-linear correlation coefficients of all 20 individual amino acids in the respiratory chain complexes of 218 animal species with log lifespan. For the unbiased comparative ranking of lifespan associations, comparisons by group were performed in addition to these correlational analyses. To this end, the group of 20 reference animals was arbitrarily divided in two groups of 10 longer-lived (15–100 years) and 10 shorter-lived (0.05–6 years) animals at the median longevity. Amino acid percentages were then directly compared without log-transformation, as in all other paradigms.
All database sequences were numerically evaluated using customized Perl scripts (http://www.perl.org) generated with the BioPerl package (http://www.bioperl.org) for computational biology [42]. Multilinear modeling for the prediction of lifespan (as variable) by the frequencies of all 20 amino acids was performed with IBM SPSS software (version 23), adopting the concomitant inclusion of all 20 amino acids as potential predictors. The resulting fractional amino acid contributions to the final model (q values in %) are provided as an outcome.

2.3. Transmembrane Domain Content Correction

Normalization for identical transmembrane domain (TMD) content was applied to the comparison of peroxisomal proteins versus all proteins (paradigm I), mitochondrial proteins versus all proteins (paradigm II), and mitochondrially encoded respiratory chain proteins versus all mitochondrial proteins (paradigm III). In these analyses, different proteomes with varying TMD contents (ranging from 2% to 53% in Homo sapiens, for example) were compared numerically, potentially introducing bias, especially in regard to the hydrophobic amino acids (which are generally accumulated in TMDs) and the highly polar and charged amino acids (which are generally decumulated in TMDs). To avoid any such bias, the pooled TMDs and the pooled non-TMDs from each proteome were analyzed for their amino acid usage separately and subsequently compared to the pooled TMDs and the pooled non-TMDs from the according reference proteome in the denominator. TMDs were determined with a hidden Markov model as implemented in the software application TMHMM 2.0 [43], which has been validated for the analysis of complete proteomes [44,45]. The output of TMHMM 2.0 was further processed using Perl scripts [7].

2.4. Body Mass and Phylogenetic Corrections

Body mass corrections were done by standard partial linear correlation. Body mass data for vertebrates were predominantly obtained from the AnAge database (genomics.senescence.info/species) [46]. Invertebrate body masses were assembled individually from zoological articles and monographs, or estimated from allometric body length-body mass relationships, as detailed in [7]. Phylogenetic interdependence corrections were realized through the determination of phylogenetically independent contrasts, followed by log-linear correlation of the obtained independent contrasts with log lifespan. Details on the construction of the phylogenetic trees can be found elsewhere [7].

2.5. Statistical Analysis

Quantitative comparisons were evaluated by non-parametric ANOVAs (Kruskal–Wallis tests) because the equal variance criterion for parametric ANOVAs was not met in many of the performed comparisons. ANOVAs were generally conducted with Systat SigmaPlot (version 11). For the linear correlations, Pearson correlation coefficients and p values (from two-sided tests) were determined with SPSS. Regarding the presentation of the data, the results of all 20 proteinogenic amino acids were listed, of which the 4 strongest effects, i.e., the 4 largest deviations from unity irrespective of the direction (up or down), were visually highlighted with color. The according p values of these four amino acids were also highlighted, but only if they were lower than p = 0.001. This low threshold was arbitrarily chosen to account for the multiple testing of 20 amino acids in all comparisons. Through this procedure, both effect size and effect significance were intended to be considered, but with a hierarchy to emphasize mean effects over distributional effects in the presentation of the results.

3. Results and Discussion

3.1. Amino Acid Usage Differences between Cellular Compartments

Within the eukaryotic cell, different compartments with substantially increased production of reactive oxygen species (ROS) can be distinguished. Specifically, peroxisomes are characterized by their notorious generation of hydrogen peroxide [24], whereas mitochondria constitute the predominant source of superoxide in the cell, which is mostly converted to hydrogen peroxide subsequently [26]. Within the mitochondria, the inner membrane-bound respiratory chain complexes are the prevailing local source of superoxide [27,47], providing a rationale for the separate consideration of the mitochondrial proteome as a whole, and of mitochondrially encoded respiratory chain complex subunits. Other sites of functional oxidant production do exist in eukaryotic cells, such as the endoplasmic reticulum with its protein folding-associated flavoprotein Ero1 [48], or the different sites of expression of the various NAPDH oxidases [49], but those were not investigated in the current study.
The aforesaid proteomes were collected for a sample of 20 reference animal species with high-quality proteomic coverage, which in addition were intended to represent a wide range of lifespans and several different phylogenetic groups, even if a predominance of mammals was unavoidable. The final species selection and the sequence numbers associated with each species and compartment are provided in Table 1.
Amino acid frequencies in all protein sets of Table 1 were determined, related to each other within each animal species, and then averaged over all species to calculate generalized usage ratios for each amino acid. The following usage ratios were inquired: peroxisomal proteins versus all proteins (paradigm I), mitochondrial proteins versus all proteins (paradigm II), and mitochondrially encoded respiratory chain proteins versus all mitochondrial proteins (paradigm III). These amino acid usage ratios are listed in Table 2, together with a statistical assessment by non-parametric ANOVA on ranks. In all tables in this work, the four most pronounced deviations from unity (i.e., ratio = 1.0 or correlation coefficient r = 0.0) are color-highlighted, together with the associated p value in case of significance at the p = 0.001 level. These thresholds were arbitrarily chosen, based on the number of comparisons (n = 20 amino acids) and the observed distribution of the data, which was non-normal in a substantial fraction of cases. For each species and paradigm, both protein sets to be compared were adjusted to the same transmembrane domain (TMD) content because the analyzed proteins sets were in some cases characterized by substantially different shares of TMDs. For example, all mitochondrial proteins had about 5% TMDs, whereas mitochondrially encoded proteins had about 55% TMDs, with an accordingly higher content of hydrophobic residues such as leucine and valine. In order to avoid any bias resulting from this disparity, especially any erroneous “accumulation” of hydrophobic amino acids in proteomes that were concomitantly oxidant-exposed and membrane-rich, all proteomes were individually normalized to identical TMD contents.
As evident from Table 2, the amino acid cysteine was among the four most altered amino acids in all three paradigms, and its modulation was always negative and always reached statistical significance. Concordant effects in at least two paradigms were an increase in methionine in paradigms II and III, and a loss of proline and serine in paradigms I and II. Notable effects at the level of a single paradigm were an increase of isoleucine in paradigm I, which was the strongest signal in peroxisomes, and an accumulation of asparagine as well as a loss of arginine in paradigm III (respiratory chain complexes), the latter of which occurred despite correction for equal TMD content.
Since paradigm III investigates proteomes derived from mitochondrial translation systems occasionally using a different genetic code, the associated decoding changes are also provided in Table 2. Qualitatively, methionine accumulation was related to a genetic code change that has been shown to be adaptive and causative to the accumulation [18]. In addition, arginine loss and tryptophan gain were also associated with the deviant genetic decoding in animal mitochondria, which use two codons (UGA and UGG) throughout to encode tryptophan instead of one codon (the canonic UGG codon). Moreover, 13 out of the 20 selected animal species do not encode arginine by AGA and AGG, reducing the number of arginine codons from 6 canonic codons to 4 codons. Both changes were apparently not buffered by compensatory changes in codon usage, as they both contributed to significantly altered amino acid frequencies at the protein level (Table 2). After a long debate on the origin of the modern genetic code variants [50,51], there is compelling evidence now to indicate that deviant genetic codes like those in animal mitochondria primarily evolve to match the amino acid requirements of the encoded proteins [18,51,52] rather than having their foundation at the DNA or RNA level. Interestingly, the loss of the AUA codon for Ile (in 19 out of 20 species) was evidently compensated, and did not result in an overall loss of isoleucine at the protein level. Without compensation, one would have expected such an isoleucine loss if methionine accumulation were the selective driver of the AUA Ile-to-Met transition. No attempts were made in this work to formally adjust for coding differences across species because the final outcome at the protein level was regarded as the relevant biochemical parameter to consider when studying effects of redox reactivity at the protein level. The individual mechanisms by which amino acid changes are installed over evolutionary time in vivo are certainly interesting and multi-dimensional (e.g., individual mutations, collective GC-content shifts, genetic code changes), but without immediate relevance to the interpretation of the final number of reactive chemical groups in proteins, which was the primary aim of the current comparative investigation.
A graphical representation of the quantitative data of Table 2 is provided in Figure 1. Therein, fractional increases or decreases are plotted for each amino acid as a bar stack, in which each bar (as red or green line) represents 1 (out of 20) animal species. The 20 animal species were sorted by decreasing lifespan from left to right analogously, as listed in Table 1. Several aspects in this comparison across paradigms and amino acids are noteworthy. First, amino acid usage changes in mitochondrially encoded proteins (versus all mitochondrial proteins, paradigm III) were generally more pronounced, often reaching a factorial accumulation/depletion by a factor of 2 (or ½). In the other paradigms, the approximate equivalence was a factorial accumulation/depletion by a factor of 1.2 (or 1/1.2). Still, interspecies variation (i.e., the variability within each bar stack) was largely unrelated to the degree of accumulation/depletion, but it appeared to be somewhat higher in peroxisomes, which is likely attributable to the smaller number of proteins in several representatives of this groups (Table 1). The most interesting outcome of the vertical comparison of the three plotted paradigms is yet their similarity, particularly as regards the prototypically redox-active amino acids cysteine, methionine, tyrosine and tryptophan: cysteine was strongly decreased (red color in Figure 1), whereas methionine, tyrosine and tryptophan were increased throughout (green color in Figure 1). As further discussed below, this behavior matches well with the inferred functional consequences of amino acid (one-electron) oxidation in (membrane) proteins as derived from chemical reaction rate constants [12].

3.2. Amino Acid Usage Alterations in Relation to Longevity

In the majority of cases shown in Figure 1, amino acid changes were similar across species and, thus, lifespans; in other words, most bar stacks had a more or less horizontal upper or lower borderline. However, the amino acid cysteine in paradigm III displayed a borderline that pronouncedly inclined to the left, denoting a more severe usage bias in long-lived animal species. This signal has been studied before [7]; however, a formal comparison with all other amino acids and across paradigms has been lacking. Hence, as a fourth paradigm, the correlation of each of the amino acid accumulation/depletion signals with lifespan was investigated.
Table 3 presents the results obtained by a linear correlation analysis of double log-transformed data (log amino acid usage ratio versus log lifespan in years). In the three paradigms I–III, several amino acids reached correlation coefficients of 0.5 or higher in absolute values. However, only one signal in paradigms I and II reached statistical significance, namely asparagine, which was increased in peroxisomal proteins of long-lived species, as directly visible from Figure 1. In all other cases, apparently inclined bar stacks turned out to be insignificantly associated with longevity due to inter-species variation or small effect sizes. In mitochondrially encoded proteins, four amino acid alterations were significantly associated with longevity, namely cysteine loss, which was the strongest correlation, followed by threonine, glutamine and proline (Table 3), all of which were increased with longevity. Of note, selection for threonine has been described before as a longevity-related adaptation in long-lived mammals and explained with structural, non-redox properties of this amino acid in transmembrane domains [32,45]. The same studies have also noted the negative cysteine signal.
Parametric longevity correlations offer the potentially insightful opportunity to compare molecular correlations of interest to other inherent correlatives of longevity such as body mass, which is the undisputed master predictor of lifespan in animals [53], or phylogenetic co-segregation. Both correlatives can be corrected for, either to calculate the degree by which a molecular correlation is better than the master correlation with body mass, or to calculate the degree by which a molecular trait is better correlated with longevity than expected from species phylogeny. These corrections are usually realized by partial correlation, or by phylogenetically independent contrast analysis [54] followed by correlation [7,20,45]. Still, it is important to note that even if a correlation were not significant any more after any of these corrections, it may still be causally and relevantly involved in the limitation of longevity, e.g., in the standard case of molecular co-segregation [55].
To adopt both of these corrections, body masses were collected, and a phylogenetic tree was assembled for the animal species under investigation, as described in [7]. Both corrections were applied to the paradigm with the strongest signals at baseline, namely longevity within the respiratory chain complex signals (denoted paradigm IVc in Table 3). The resulting corrected correlation coefficients and the ensuing p values are provided in Table 4. As a further analysis presented in Table 4, a multilinear predictive modelling of lifespan using all 20 amino acids was performed, followed by inspection of the relative contribution of each amino acid to the generated model (as q value in %). Finally, to better consider effect sizes (which are sometimes obfuscated in correlation analyses), simple binomial ratios were calculated after division of the 20 animal species into 2 groups of 10 species each, encompassing short-lived (0.05–6 years) and long-lived (15–100 years) animals (Table 1).
Body mass correction led to substantially reduced correlation coefficients throughout, turning all single amino acid associations statistically insignificant (Table 4). Phylogenetic correction had a comparable effect, with the exception of cysteine, the negative association of which with longevity was preserved. Notably, cysteine’s correlation coefficient (r = −0.791) still came close to the correlation coefficient of body mass with longevity (r = −0.886), similar to before in the uncorrected data. This is an intriguing result because a single amino acid in a rather small proteome would not be expected to predict a complex trait like lifespan almost equally well as the complex reference trait body mass, which is inevitably linked to lifespan on multiple levels. After all, the building of a large body requires time due to the inherent temporal limits of cell division and bioenergetic uptake. Cysteine also emerged as prime lifespan-associated amino acid from multilinear modelling, contributing 43% to the final correlation, followed by asparagine (12%). In the effect size calculations, cysteine again ranked as the most differentially used amino acid in long-lived versus short-lived species, followed by threonine (increased), glutamine (increased), and phenylalanine (decreased). Of the latter amino acids, glutamine and threonine had also been marked as significantly altered in the same direction in the correlational analysis (Table 4).
The above observations were put on probation by examining a second larger dataset of 218 animal species with known body mass and phylogenetic positioning. Due to the lack of fully sequenced proteome data for the majority of these species, only the smaller mitochondrially encoded proteomes were assessed (i.e., the numerators of the fractions represented by paradigm IVc). The results in Table 5 demonstrate a large overlap with the findings obtained before in the internally controlled 20-species reference set. Specifically, the negative correlation of cysteine with lifespan was reproduced, as were the positive correlations of proline and threonine. In addition, a negative phenylalanine association emerged that had just failed to reach significance in the 20-species cohort. Uniquely, the cysteine-lifespan correlation remained significant after corrections for body mass and lifespan as reported in [7], and hence in all three correlational analyses (Table 5). Proline was the second-best predictor of lifespan in this dataset, as it reached significance in two out of three correlations and ranked second in the multilinear model. Cytosolic proline in yeast [56], as well as several other soluble amino acids in animals, have been shown to be related to longevity or cellular senescence before [57,58]. In plants, cytosolic proline may act as an indirect quencher of different oxidants [59]. Still, it is unclear in how far such data for soluble molecules may be extrapolated to the much less dynamic amino acid contents of proteins. In fact, as exemplified for cysteine, soluble chemical groups may exhibit exactly converse associations with longevity than the same chemical groups in macromolecules because the sacrifice of the former may help protect the latter [6,12].

3.3. Amino Acid Usage Adaptations Linked to an Aerobic Lifestyle

Parasitic helminths have long been known to avoid oxygen-dependent ATP generation and largely depend on anaerobic fermentation after infection of their host, even in cases in which oxygen partial pressures would be more than sufficient for aerobic respiration [36,37]. The most plausible origin of this initially unexpected behavior is the potential toxicity arising from reactive oxygen species (ROS) produced during respiration, to which many parasitic species are poorly shielded when compared to their vertebrate hosts [60,61]. Notably, infections with facultatively anaerobic eukaryotes are in many instances effectively treated with prooxidative drugs such as artemisinin [62], metronidazole [63] or chloroquine [64].
To examine potential differences in the proteomic composition of aerobic and facultatively anaerobic animal species, a collection of mitochondrially encoded proteomes was assembled, encompassing 55 facultatively anaerobic (parasitic) and 31 aerobic (free-living) species predominantly of the nematode and platyhelminth phyla [8]. Table 6 shows the outcome of this comparison. Global cysteine avoidance in aerobes was the most distinct difference in terms of effect size, whereas valine avoidance and glutamine preference were the most significant signals in statistical terms. Beyond these three amino acids, several other amino acids were characterized by usage differences in aerobic versus anaerobic species of more than 50% in the analyzed highly homologous proteomes (i.e., mitochondrial respiratory chain complexes). This unexpected finding prompted us to expand the comparative investigation of the proteomic footprints of aerobicity to larger proteomes. With only four fully sequenced proteomes of aerobic helminths as per the end of 2023, and two of them being of the same genus, Caenorhabditis, total proteomic comparison of helminths was regarded unreliable. Hence, three groups of prokaryotes with well-defined oxygen preferences were collected to augment the data generated in helminths. More specifically, groups of approximately 30–50 archaea (paradigm VI), Gram-positive bacteria (paradigm VII) and Gram-negative bacteria (paradigm VIII) were classified according to their oxygen preferences, followed by analysis of their total proteomic amino acid usage patterns. The score provided in Table 7 indicates that effect sizes and significance levels were smaller than in the helminths, even if certain signals were in fact recapitulated. Cysteine avoidance in Gram-positive bacteria pursuing aerobic respiration was the most pronounced and significant change seen across all groups, followed by glutamine accumulation in Gram-negative bacteria. Aerobic bacteria also elicited increased tryptophan usage. As an entirely new signal absent in animals, lysine avoidance was noted in all groups of aerobic prokaryotes. Still, neither tryptophan nor lysine reached statistical significance in any individual group.
Despite the high diversity of prokaryotes in terms of their metabolic activities, various observations support their usefulness as evolutionary indicators of oxidative stress. For instance, strains of Escherichia coli have been shown to respond to oxidative stress through adaptations in their structural proteomes [65]. Significant differences have also been noted for the proteomic amino acid usage of 295 aerobic versus anaerobic prokaryotes after correction for phylogenetic interdependence, namely increases in tryptophan and threonine and a decrease in cysteine [66]. Comparisons of Escherichia coli with (highly radioresistant) Deinococcus radiodurans have suggested that lower cysteine and methionine contents were decisive for the resistance of the latter species, and that lysine increased the accessibility of neighboring redox-reactive residues to oxidants, leading to its elimination. Tryptophan and tyrosine were likewise significantly enriched in Deinococcus radiodurans [67]. Overall, these patterns clearly resemble the trends listed in Table 7.

3.4. Relation of the Observed Amino Acid Usage Differences to Amino Acid Redox Reactivity

To enable a comparative survey of the detected amino acid usage alterations in the individual paradigms, the measured effect sizes of all changes were ranked from 1–20, with 1 being the most heavily altered amino acid in each comparison and 20 being the least altered amino acid. The resulting score for the animal paradigms I–V is provided in Table 8, with increases in green and decreases in red; the score for the prokaryotes is provided in Table 9. Amino acids that were regarded as potentially altered in two or more paradigms (and thus highlighted in Table 1, Table 2, Table 3, Table 4, Table 5 and Table 6) are further discussed in the following.
Cysteine avoidance in mitochondrial proteins has been described in different datasets before [7,8], but its extent, generalizability and relative relevance compared to other amino acid changes under oxidative stress on evolutionary timescales have not been systematically evaluated. From the present meta-examination, it becomes clear that the general susceptibility of protein cysteine to oxidation is in fact the primary subject of proteomic adaptations to persistently increased levels of reactive oxygen species (ROS). Table 10 provides an overview of the effectuated quantitative amino acid frequency changes for cysteine and the prototypically redox-active amino acids methionine, tryptophan and tyrosine. In three unrelated paradigms (III, V and VII), cysteine was decreased by approximately 50% on a proteome-wide scale, representing the largest effect size of all amino acids.
Independently of specific redox assumptions, a very low degree of conservation of isolated, solitary surface cysteines has been evidenced before, which was contrasted by maximum conservation of potentially disulfide-forming cysteine pairs [13,14]. Even if no separate analysis of solitary versus paired cysteines was attempted in the current investigation, it is highly plausible to suppose that the main component of the redox-related cysteine bias described herein is attributable to solitary cysteines. If this were correct, the global cysteine losses of Table 10 would imply even higher fractional losses when regarding only solitary cysteines. The large frequency changes obtained in paradigms III and IV provide support for this interpretation because both paradigms relate to respiratory chain complexes that do not possess any disulfide bridges and harbor only one conserved cysteine pair that ligates copper in respiratory chain complex IV [7].
Since redox-related cysteine avoidance, as described herein, likely affects solitary surface cysteines [13] and generally unpaired transmembrane cysteines [8] in a selective fashion, it has been proposed that the one-electron oxidation of cysteine thiol groups to thiyl radicals may mediate these features of proteomic adaptation [6]. Detailed analyses of the reaction rate constants of thiyl radicals in the context of lipid peroxidation have supported this view exemplarily [12]. In the current survey, cysteine was the most or second-most heavily affected amino acid in 7 out of 8 paradigms, with Gram-negative bacteria (paradigm VIII) being the only exception (Table 7 and Table 9). As Gram-negative bacteria usually possess the capacity to synthesize ample amounts of dedicated antioxidant thiols like glutathione [68], which contrasts archaea and Gram-negative bacteria, our data strongly suggest that the accumulation of soluble thiols in Gram-negative bacteria enabled the continued high-frequency usage of cysteine as protein structural component despite rising oxygen concentrations or intensified oxygen utilization [69]. This interpretation is chemically plausible, as soluble thiols show very high reactivity towards thiyl radicals to generate mostly harmless, mixed disulfide radicals [6,11].
Adaptive increases in proteomic methionine in oxidant-exposed proteomes have also been observed before, with the conclusion that such increases in methionine essentially mark particularly vulnerable, biologically precious proteins rather than sites with the highest oxidant production [18,21]. This idea fits to the observed antioxidant chemistry and reactivity of methionine, which reacts preferentially with non-radical species by two-electron reduction rather than blunting radical species through one-electron reduction [12]. No association of methionine with longevity was observed in any of our datasets (Table 3, Table 4 and Table 5). This stands in partial contrast to earlier studies focusing on mammals, in which, after phylogenetic correction, an otherwise invisible signal had emerged [19,20]. Potentially, methionine-longevity adaptations occur exclusively in higher vertebrate clades. Notably, negative changes in the soluble level of methionine and various metabolites of the methionine degradation pathway have been found to be prominent biochemical predictors of longevity and stress resistance [33,58,70]. Even if soluble amino acid concentrations should not have any immediate influence on their usage in genetically encoded proteins, a coevolution of both properties is clearly conceivable [33,70].
Serine was significantly reduced in two paradigms and lowered by trend in several other comparisons (Table 2, 5 and 6). Even if serine is not regarded as a typical redox-active amino acid, it has been shown to be one of the fastest amino acids to react with neighboring thiyl radicals, carrying their detrimental reactivity into the interior of the protein [6,10]. Potentially, this secondary effect underlies the avoidance of serine in specific situations. In interesting contrast, the chemically related amino acid threonine was unaltered throughout, except for paradigms investigating respiratory chain complexes and their adaptations towards increased longevity (Table 4 and Table 5). Strikingly, the phenomenon of an adaptive respiratory chain threonine increase in long-lived animals has been characterized in detail before, proposing a structural rather than reactivity origin of this alteration in terms of a selective hydrophobicity modulation [31,45]. The singularity and unusual strength of the current threonine signal underline this view.
Significant avoidance of proline was observed in two paradigms (I and II), but proline was otherwise unchanged or increased in all other comparisons, the latter especially in relation to longevity (Table 4 and Table 5). The interpretation of these findings is unclear for now. Proline shows a relatively high degree of hydroperoxide formation after attack by hydroxyl radicals [4]. However, comparable quality targets of hydroperoxide modification were either unchanged (i.e., Leu, which was the most inert amino acid across all paradigms in this study), or modulated in opposite directions (i.e., Val, showing substantial peaks up (paradigm I and II) and down (paradigm III and V)), or selectively increased (i.e., Ile (paradigm I)). Hence, the changes in proline and the branched-chain hydrophobic amino acids may rather be structural-individual rather than stemming from general chemical reactivity patterns, and may perhaps also be attributable to selective hydrophobicity modulation [31,45]. Of note, together with methionine, valine was among the most meaningful amino acids that, in respiratory chain sites of 324 vertebrates, individually predicted longevity [71]. In any case, the high volatility of valine compared to the structurally related leucine (Table 8) is certainly interesting and requires further study. The same applies to the so-far unexplained negative correlation of the soluble amino acid levels of proline with lifespan [58].
Finally, asparagine and glutamine were significantly increased or showed indicative trends in a number of evaluations, including peroxisome-longevity (paradigm IVa), respiratory chain-longevity (paradigm IVc), and aerobicity (paradigm V) (Table 3, Table 4 and Table 6). By most measures, asparagine and glutamine are hardly oxidation-prone amino acids [4,22,72], such that their accumulation in peroxisomal proteins of longer-lived animals is very likely unrelated to their redox properties. Somewhat paradoxically, though, asparagine is one of the most easily isomerized and hydrolyzed amino acids, resulting in the formation of isoaspartate and D-aspartate in proteins [73]. Both isoaspartate and D-aspartate are well-established markers of protein aging and reflect organismal aging in particularly long-lived proteins such as eye lens crystallins [74]. Hence, one would have rather expected the avoidance of asparagine and perhaps other hydrolysis-prone amino acids in proteins of long-lived species. As a bona fide interpretation of this discrepancy, spontaneous isomerization of asparagine and glutamine may have an only softly adverse effect on organismal longevity that is perhaps readily overcompensated by certain structural advantages of the amino acids themselves. One of these advantages could be their chaotropic solubility-enhancing activity on protein surfaces [22], preventing oxidized or otherwise denatured proteins from aggregation, a common cause of age-related disease.

4. Conclusions

In a multi-panel comparative analysis of global amino acid usage in animal proteomes exposed to different levels and types of intrinsic oxidative stress, we have identified cysteine as the most universally and pronouncedly affected amino acid. In contrast to the observed avoidance of cysteine, the other canonically redox-active amino acids methionine, tyrosine and tryptophan were widely increased, but they only reached significance in a limited number of comparisons, notably methionine in mitochondrial proteomes. As regards all other amino acids, relevant changes in more than one comparison were restricted to increases in glutamine and decreases in serine and proline. For these cases, the literature data and specificity considerations mostly suggest non-redox or only indirectly redox-related causes for their altered occurrence. Cysteine is the only amino acid the presence of which in proteins appears to be evolutionarily suppressed under conditions of oxidative stress.

Author Contributions

Conceptualization, M.S. and B.M.; methodology and tools, M.S.; formal analysis, M.S. and B.M.; data curation, M.S. and B.M.; writing—draft, M.S. and B.M.; writing—editing, B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article. Customized Perl scripts are available from the authors upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Beckman, K.B.; Ames, B.N. The free radical theory of aging matures. Physiol. Rev. 1998, 78, 547–581. [Google Scholar] [CrossRef] [PubMed]
  2. Pamplona, R.; Barja, G. Highly resistant macromolecular components and low rate of generation of endogenous damage: Two key traits of longevity. Ageing Res. Rev. 2007, 6, 189–210. [Google Scholar] [CrossRef] [PubMed]
  3. Golubev, A.; Hanson, A.D.; Gladyshev, V.N. A tale of two concepts: Harmonizing the free radical and antagonistic pleiotropy theories of aging. Antioxid. Redox Signal. 2018, 29, 1003–1017. [Google Scholar] [CrossRef] [PubMed]
  4. Davies, M.J. Protein oxidation and peroxidation. Biochem. J. 2016, 473, 805–825. [Google Scholar] [CrossRef] [PubMed]
  5. Kunath, S.; Schindeldecker, M.; De Giacomo, A.; Meyer, T.; Sohre, S.; Hajieva, P.; von Schacky, C.; Urban, J.; Moosmann, B. Prooxidative chain transfer activity by thiol groups in biological systems. Redox Biol. 2020, 36, 101628. [Google Scholar] [CrossRef] [PubMed]
  6. Moosmann, B.; Hajieva, P. Probing the role of cysteine thiyl radicals in biology: Eminently dangerous, difficult to scavenge. Antioxidants 2022, 11, 885. [Google Scholar] [CrossRef] [PubMed]
  7. Moosmann, B.; Behl, C. Mitochondrially encoded cysteine predicts animal lifespan. Aging Cell 2008, 7, 32–46. [Google Scholar] [CrossRef] [PubMed]
  8. Schindeldecker, M.; Stark, M.; Behl, C.; Moosmann, B. Differential cysteine depletion in respiratory chain complexes enables the distinction of longevity from aerobicity. Mech. Ageing Dev. 2011, 132, 171–179. [Google Scholar] [CrossRef]
  9. Kunath, S.; Moosmann, B. What is the rate-limiting step towards aging? Chemical reaction kinetics might reconcile contradictory observations in experimental aging research. Geroscience 2020, 42, 857–866. [Google Scholar] [CrossRef]
  10. Nauser, T.; Pelling, J.; Schöneich, C. Thiyl radical reaction with amino acid side chains: Rate constants for hydrogen transfer and relevance for posttranslational protein modification. Chem. Res. Toxicol. 2004, 17, 1323–1328. [Google Scholar] [CrossRef]
  11. Nauser, T.; Koppenol, W.H.; Schöneich, C. Protein thiyl radical reactions and product formation: A kinetic simulation. Free Radic. Biol. Med. 2015, 80, 158–163. [Google Scholar] [CrossRef]
  12. Hajieva, P.; Abrosimov, R.; Kunath, S.; Moosmann, B. Antioxidant and prooxidant modulation of lipid peroxidation by integral membrane proteins. Free Radic. Res. 2023, 57, 105–114. [Google Scholar] [CrossRef]
  13. Marino, S.M.; Gladyshev, V.N. Cysteine function governs its conservation and degeneration and restricts its utilization on protein surfaces. J. Mol. Biol. 2010, 404, 902–916. [Google Scholar] [CrossRef] [PubMed]
  14. Marino, S.M.; Gladyshev, V.N. Analysis and functional prediction of reactive cysteine residues. J. Biol. Chem. 2012, 287, 4419–4425. [Google Scholar] [CrossRef] [PubMed]
  15. Winterbourn, C.C. Superoxide as an intracellular radical sink. Free Radic. Biol. Med. 1993, 14, 85–90. [Google Scholar] [CrossRef]
  16. Levine, R.L.; Mosoni, L.; Berlett, B.S.; Stadtman, E.R. Methionine residues as endogenous antioxidants in proteins. Proc. Natl. Acad. Sci. USA 1996, 93, 15036–15040. [Google Scholar] [CrossRef]
  17. Weissbach, H.; Etienne, F.; Hoshi, T.; Heinemann, S.H.; Lowther, W.T.; Matthews, B.; St John, G.; Nathan, C.; Brot, N. Peptide methionine sulfoxide reductase: Structure, mechanism of action, and biological function. Arch. Biochem. Biophys. 2002, 397, 172–178. [Google Scholar] [CrossRef] [PubMed]
  18. Bender, A.; Hajieva, P.; Moosmann, B. Adaptive antioxidant methionine accumulation in respiratory chain complexes explains the use of a deviant genetic code in mitochondria. Proc. Natl. Acad. Sci. USA 2008, 105, 16496–16501. [Google Scholar] [CrossRef]
  19. Aledo, J.C. Methionine in proteins: The Cinderella of the proteinogenic amino acids. Protein Sci. 2019, 28, 1785–1796. [Google Scholar] [CrossRef]
  20. Aledo, J.C.; Li, Y.; de Magalhães, J.P.; Ruíz-Camacho, M.; Pérez-Claros, J.A. Mitochondrially encoded methionine is inversely related to longevity in mammals. Aging Cell 2011, 10, 198–207. [Google Scholar] [CrossRef]
  21. Schindeldecker, M.; Moosmann, B. Protein-borne methionine residues as structural antioxidants in mitochondria. Amino Acids 2015, 47, 1421–1432. [Google Scholar] [CrossRef]
  22. Moosmann, B. Redox biochemistry of the genetic code. Trends Biochem. Sci. 2021, 46, 83–86. [Google Scholar] [CrossRef] [PubMed]
  23. Goel, S.K.; Lalwani, N.D.; Reddy, J.K. Peroxisome proliferation and lipid peroxidation in rat liver. Cancer Res. 1986, 46, 1324–1330. [Google Scholar] [PubMed]
  24. Fransen, M.; Nordgren, M.; Wang, B.; Apanasets, O. Role of peroxisomes in ROS/RNS-metabolism: Implications for human disease. Biochim. Biophys. Acta 2012, 1822, 1363–1373. [Google Scholar] [CrossRef] [PubMed]
  25. Tahara, E.B.; Navarete, F.D.; Kowaltowski, A.J. Tissue-, substrate-, and site-specific characteristics of mitochondrial reactive oxygen species generation. Free Radic. Biol. Med. 2009, 46, 1283–1297. [Google Scholar] [CrossRef]
  26. Brand, M.D. Mitochondrial generation of superoxide and hydrogen peroxide as the source of mitochondrial redox signaling. Free Radic. Biol. Med. 2016, 100, 14–31. [Google Scholar] [CrossRef]
  27. Goncalves, R.L.S.; Watson, M.A.; Wong, H.S.; Orr, A.L.; Brand, M.D. The use of site-specific suppressors to measure the relative contributions of different mitochondrial sites to skeletal muscle superoxide and hydrogen peroxide production. Redox Biol. 2020, 28, 101341. [Google Scholar] [CrossRef]
  28. Suthammarak, W.; Somerlot, B.H.; Opheim, E.; Sedensky, M.; Morgan, P.G. Novel interactions between mitochondrial superoxide dismutases and the electron transport chain. Aging Cell 2013, 12, 1132–1140. [Google Scholar] [CrossRef]
  29. Pamplona, R.; Costantini, D. Molecular and structural antioxidant defenses against oxidative stress in animals. Am. J. Physiol. Regul. Integr. Comp. Physiol. 2011, 301, R843–R863. [Google Scholar] [CrossRef]
  30. Gomez, J.; Mota-Martorell, N.; Jove, M.; Pamplona, R.; Barja, G. Mitochondrial ROS production, oxidative stress and aging within and between species: Evidences and recent advances on this aging effector. Exp. Gerontol. 2023, 174, 112134. [Google Scholar] [CrossRef]
  31. Kitazoe, Y.; Kishino, H.; Hasegawa, M.; Matsui, A.; Lane, N.; Tanaka, M. Stability of mitochondrial membrane proteins in terrestrial vertebrates predicts aerobic capacity and longevity. Genome Biol. Evol. 2011, 3, 1233–1244. [Google Scholar] [CrossRef]
  32. Kitazoe, Y.; Hasegawa, M.; Tanaka, M.; Futami, M.; Futami, J. Mitochondrial determinants of mammalian longevity. Open Biol. 2017, 7, 170083. [Google Scholar] [CrossRef] [PubMed]
  33. Jove, M.; Mota-Martorell, N.; Fernàndez-Bernal, A.; Portero-Otin, M.; Barja, G.; Pamplona, R. Phenotypic molecular features of long-lived animal species. Free Radic. Biol. Med. 2023, 208, 728–747. [Google Scholar] [CrossRef] [PubMed]
  34. Imlay, J.A. Pathways of oxidative damage. Annu. Rev. Microbiol. 2003, 57, 395–418. [Google Scholar] [CrossRef] [PubMed]
  35. Yang, L.; Mih, N.; Anand, A.; Park, J.H.; Tan, J.; Yurkovich, J.T.; Monk, J.M.; Lloyd, C.J.; Sandberg, T.E.; Seo, S.W.; et al. Cellular responses to reactive oxygen species are predicted from molecular mechanisms. Proc. Natl. Acad. Sci. USA 2019, 116, 14368–14373. [Google Scholar] [CrossRef] [PubMed]
  36. Saz, H.J. Energy metabolisms of parasitic helminths: Adaptations to parasitism. Annu. Rev. Physiol. 1981, 43, 323–341. [Google Scholar] [CrossRef] [PubMed]
  37. Barrett, J. The anaerobic end-products of helminths. Parasitology 1984, 88, 179–198. [Google Scholar] [CrossRef] [PubMed]
  38. Carey, J.R.; Judge, D.S. Life Spans of Mammals, Birds, Amphibians, Reptiles, and Fish; Monographs on Population Aging No. 8; Odense University Press: Odense, Denmark, 2000. [Google Scholar]
  39. Bowler, J.K. Longevity of Reptiles and Amphibians in North American Collections as of 1 November 1975; Herpetological Circular No. 6; Society for the Study of Amphibians and Reptiles: Milwaukee, WI, USA, 1977. [Google Scholar]
  40. Oger, P.M.; Cario, A. Adaptation of the membrane in Archaea. Biophys. Chem. 2013, 183, 42–56. [Google Scholar] [CrossRef] [PubMed]
  41. Garcia, L. Anaerobic Gram-Positive Bacilli. In Clinical Microbiology Procedures Handbook, 3rd ed.; ASM Press: Washington, DC, USA, 2010; pp. 770–781. [Google Scholar]
  42. Stajich, J.E.; Block, D.; Boulez, K.; Brenner, S.E.; Chervitz, S.A.; Dagdigian, C.; Fuellen, G.; Gilbert, J.G.; Korf, I.; Lapp, H.; et al. The Bioperl toolkit: Perl modules for the life sciences. Genome Res. 2002, 12, 1611–1618. [Google Scholar] [CrossRef]
  43. Sonnhammer, E.L.; von Heijne, G.; Krogh, A. A hidden Markov model for predicting transmembrane helices in protein sequences. Proc. Int. Conf. Intell. Syst. Mol. Biol. 1998, 6, 175–182. [Google Scholar]
  44. Krogh, A.; Larsson, B.; von Heijne, G.; Sonnhammer, E.L. Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J. Mol. Biol. 2001, 305, 567–580. [Google Scholar] [CrossRef] [PubMed]
  45. Kitazoe, Y.; Kishino, H.; Hasegawa, M.; Nakajima, N.; Thorne, J.L.; Tanaka, M. Adaptive threonine increase in transmembrane regions of mitochondrial proteins in higher primates. PLoS ONE 2008, 3, e3343. [Google Scholar] [CrossRef] [PubMed]
  46. Tacutu, R.; Thornton, D.; Johnson, E.; Budovsky, A.; Barardo, D.; Craig, T.; Diana, E.; Lehmann, G.; Toren, D.; Wang, J.; et al. Human Ageing Genomic Resources: New and updated databases. Nucleic Acids Res. 2018, 46, D1083–D1090. [Google Scholar] [CrossRef]
  47. Fang, J.; Wong, H.S.; Brand, M.D. Production of superoxide and hydrogen peroxide in the mitochondrial matrix is dominated by site IQ of complex I in diverse cell lines. Redox Biol. 2020, 37, 101722. [Google Scholar] [CrossRef]
  48. Ong, G.; Logue, S.E. Unfolding the interactions between endoplasmic reticulum stress and oxidative stress. Antioxidants 2023, 12, 981. [Google Scholar] [CrossRef]
  49. Laddha, A.P.; Kulkarni, Y.A. NADPH oxidase: A membrane-bound enzyme and its inhibitors in diabetic complications. Eur. J. Pharmacol. 2020, 881, 173206. [Google Scholar] [CrossRef] [PubMed]
  50. Osawa, S.; Jukes, T.H.; Watanabe, K.; Muto, A. Recent evidence for evolution of the genetic code. Microbiol. Rev. 1992, 56, 229–264. [Google Scholar] [CrossRef]
  51. Sengupta, S.; Yang, X.; Higgs, P.G. The mechanisms of codon reassignments in mitochondrial genetic codes. J. Mol. Evol. 2007, 64, 662–688. [Google Scholar] [CrossRef]
  52. Swire, J.; Judson, O.P.; Burt, A. Mitochondrial genetic codes evolve to match amino acid requirements of proteins. J. Mol. Evol. 2005, 60, 128–139. [Google Scholar] [CrossRef]
  53. de Magalhães, J.P.; Costa, J.; Church, G.M. An analysis of the relationship between metabolism, developmental schedules, and longevity using phylogenetic independent contrasts. J. Gerontol. Ser. A Biomed. Sci. Med. Sci. 2007, 62, 149–160. [Google Scholar] [CrossRef]
  54. Felsenstein, J. Phylogenies and the comparative method. Am. Nat. 1985, 125, 1–15. [Google Scholar] [CrossRef]
  55. Barja, G. Correlations with longevity and body size: To correct or not correct? J. Gerontol. Ser. A Biomed. Sci. Med. Sci. 2014, 69, 1096–1098. [Google Scholar] [CrossRef]
  56. Nishimura, A.; Yoshikawa, Y.; Ichikawa, K.; Takemoto, T.; Tanahashi, R.; Takagi, H. Longevity regulation by proline oxidation in yeast. Microorganisms 2021, 9, 1650. [Google Scholar] [CrossRef]
  57. Canfield, C.A.; Bradshaw, P. Amino acids in the regulation of aging and aging-related diseases. Transl. Med. Aging 2019, 3, 70–89. [Google Scholar] [CrossRef]
  58. Mota-Martorell, N.; Jové, M.; Berdún, R.; Òbis, È.; Barja, G.; Pamplona, R. Methionine metabolism is down-regulated in heart of long-lived mammals. Biology 2022, 11, 1821. [Google Scholar] [CrossRef] [PubMed]
  59. Kavi Kishor, P.B.; Suravajhala, P.; Rathnagiri, P.; Sreenivasulu, N. Intriguing role of proline in redox potential conferring high temperature stress tolerance. Front. Plant Sci. 2022, 13, 867531. [Google Scholar] [CrossRef] [PubMed]
  60. Mehlotra, R.K. Antioxidant defense mechanisms in parasitic protozoa. Crit. Rev. Microbiol. 1996, 22, 295–314. [Google Scholar] [CrossRef] [PubMed]
  61. Guevara-Flores, A.; Martínez-González, J.J.; Rendón, J.L.; Del Arenal, I.P. The architecture of thiol antioxidant systems among invertebrate parasites. Molecules 2017, 22, 259. [Google Scholar] [CrossRef] [PubMed]
  62. Li, W.; Mo, W.; Shen, D.; Sun, L.; Wang, J.; Lu, S.; Gitschier, J.M.; Zhou, B. Yeast model uncovers dual roles of mitochondria in the action of artemisinin. PLoS Genet. 2005, 1, e36. [Google Scholar] [CrossRef] [PubMed]
  63. Löfmark, S.; Edlund, C.; Nord, C.E. Metronidazole is still the drug of choice for treatment of anaerobic infections. Clin. Infect. Dis. 2010, 50 (Suppl. S1), S16–S23. [Google Scholar] [CrossRef] [PubMed]
  64. Kavishe, R.A.; Koenderink, J.B.; Alifrangis, M. Oxidative stress in malaria and artemisinin combination therapy: Pros and Cons. FEBS J. 2017, 284, 2579–2591. [Google Scholar] [CrossRef] [PubMed]
  65. Mih, N.; Monk, J.M.; Fang, X.; Catoiu, E.; Heckmann, D.; Yang, L.; Palsson, B.O. Adaptations of Escherichia coli strains to oxidative stress are reflected in properties of their structural proteomes. BMC Bioinform. 2020, 21, 162. [Google Scholar] [CrossRef]
  66. Vieira-Silva, S.; Rocha, E.P. An assessment of the impacts of molecular oxygen on the evolution of proteomes. Mol. Biol. Evol. 2008, 25, 1931–1942. [Google Scholar] [CrossRef] [PubMed]
  67. Chang, R.L.; Stanley, J.A.; Robinson, M.C.; Sher, J.W.; Li, Z.; Chan, Y.A.; Omdahl, A.R.; Wattiez, R.; Godzik, A.; Matallana-Surget, S. Protein structure, amino acid composition and sequence determine proteome vulnerability to oxidation-induced damage. EMBO J. 2020, 39, e104523. [Google Scholar] [CrossRef]
  68. Copley, S.D.; Dhillon, J.K. Lateral gene transfer and parallel evolution in the history of glutathione biosynthesis genes. Genome Biol. 2002, 3, 0025. [Google Scholar] [CrossRef]
  69. Moosmann, B.; Schindeldecker, M.; Hajieva, P. Cysteine, glutathione and a new genetic code: Biochemical adaptations of the primordial cells that spread into open water and survived biospheric oxygenation. Biol. Chem. 2020, 401, 213–231. [Google Scholar] [CrossRef] [PubMed]
  70. Lewis, K.N.; Rubinstein, N.D.; Buffenstein, R. A window into extreme longevity; the circulating metabolomic signature of the naked mole-rat, a mammal that shows negligible senescence. Geroscience 2018, 40, 105–121. [Google Scholar] [CrossRef]
  71. Mariadassou, M.; Pellay, F.X. Identification of amino acids in mitochondrially encoded proteins that correlate with lifespan. Exp. Gerontol. 2014, 56, 53–58. [Google Scholar] [CrossRef]
  72. Stadtman, E.R. Protein oxidation and aging. Free Radic. Res. 2006, 40, 1250–1258. [Google Scholar] [CrossRef]
  73. Geiger, T.; Clarke, S. Deamidation, isomerization, and racemization at asparaginyl and aspartyl residues in peptides. Succinimide-linked reactions that contribute to protein degradation. J. Biol. Chem. 1987, 262, 785–794. [Google Scholar] [CrossRef]
  74. Fujii, N.; Takata, T.; Fujii, N.; Aki, K.; Sakaue, H. D-Amino acids in protein: The mirror of life as a molecular index of aging. Biochim. Biophys. Acta Proteins Proteom. 2018, 1866, 840–847. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Graphical representation of the factorial accumulation (positive values, green color) or depletion (negative values, red color) of each amino acid in paradigms I, II and III. Each of the 20 investigated animal species in these paradigms is represented by a single line within every 20-line stack per amino acid. Animals were sorted by decreasing longevity from left to right, as per Table 1 (i.e., from Homo sapiens, left, to Caenorhabditis elegans, right).
Figure 1. Graphical representation of the factorial accumulation (positive values, green color) or depletion (negative values, red color) of each amino acid in paradigms I, II and III. Each of the 20 investigated animal species in these paradigms is represented by a single line within every 20-line stack per amino acid. Animals were sorted by decreasing longevity from left to right, as per Table 1 (i.e., from Homo sapiens, left, to Caenorhabditis elegans, right).
Antioxidants 13 00267 g001
Table 1. The investigated animals and their lifespans, body masses, and sampled sequence numbers.
Table 1. The investigated animals and their lifespans, body masses, and sampled sequence numbers.
SpeciesPhylogenetic GroupLifespan (Years)Body Mass (kg)Analyzed Sequences
Full proteomeMitochondrionPeroxisome
Homo sapiensMammalia1007015,072978238
Equus caballusMammalia62100016,148283100
Pan troglodytesMammalia604514,349376168
Macaca mulattaMammalia35815,939314145
Bos taurusMammalia3075028,590749110
Canis familiarisMammalia304014,903429227
Gallus gallusAves302.615,525301101
Strongylocentrotus purpuratusEchinoidea200.125,69630929
Oryctolagus cuniculusMammalia181.814,80826990
Cavia porcellusMammalia150.7316,41920459
Monodelphis domesticaMammalia60.1116,59327346
Mus musculusMammalia60.02121,5511186188
Apis melliferaInsecta59 × 10−5852920419
Danio rerioPisces55 × 10−425,04958862
Rattus norvegicusMammalia40.326,060799158
Ciona intestinalisAscidiae21 × 10−411,65320325
Tribolium castaneumInsecta22.5 × 10−6934019355
Drosophila melanogasterInsecta0.151 × 10−6974816245
Anopheles gambiaeInsecta0.072.5 × 10−711,20012636
Caenorhabditis elegansNematoda0.052 × 10−923,86521939
Table 2. Numerical evaluation of the amino acid usage ratios in paradigms I, II, and III.
Table 2. Numerical evaluation of the amino acid usage ratios in paradigms I, II, and III.
Amino AcidPeroxisome (Paradigm I)Mitochondrion (Paradigm II)Respiratory Chain (Paradigm III)
RatiopRatiopRatiopDecoding *
A1.093 × 10−41.153 × 10−50.571 × 10−7
C0.863 × 10−50.741 × 10−70.411 × 10−7
D0.918 × 10−50.934 × 10−40.742 × 10−7
E0.913 × 10−70.922 × 10−60.636 × 10−8
F1.122 × 10−61.085 × 10−51.281 × 10−2
G1.099 × 10−51.124 × 10−40.731 × 10−6+AGA, +AGG (1/20)
H0.971 × 10−20.903 × 10−61.363 × 10−5
I1.191 × 10−41.122 × 10−31.124 × 10−3−AUA (19/20)
K0.977 × 10−21.042 × 10−10.794 × 10−5−AAA (1/20)
L1.051 × 10−21.066 × 10−31.112 × 10−7
M1.112 × 10−41.161 × 10−62.142 × 10−7+AUA (19/20)
N1.032 × 10−10.908 × 10−31.939 × 10−8+AAA (1/20)
P0.863 × 10−40.863 × 10−41.283 × 10−3
Q0.889 × 10−60.934 × 10−40.812 × 10−5
R0.922 × 10−31.032 × 10−20.496 × 10−8−AGA, −AGG (13/20)
S0.872 × 10−70.826 × 10−81.293 × 10−6+AGA, +AGG (6/20)
T0.996 × 10−10.971 × 10−11.487 × 10−3
V1.141 × 10−61.101 × 10−60.536 × 10−8
W1.102 × 10−31.053 × 10−21.816 × 10−8+UGA (20/20)
Y1.104 × 10−31.061 × 10−21.336 × 10−6
p Values derive from non-parametric ANOVA on ranks analyses. Red color highlights the 4 most pronounced effect sizes among the 20 amino acids and their significance. The according p values were also highlighted if a predefined threshold of p < 0.001 was reached. Amino acids are given in single-letter code. * When amino acids were encoded in mitochondria by more or fewer codons than in the standard genetic code, those codons are given in this column. The numbers in brackets indicate how many of the examined 20 animal species use the non-standard decoding.
Table 3. Quantitative correlation with lifespan of the signals in paradigms I, II and III (paradigm IV).
Table 3. Quantitative correlation with lifespan of the signals in paradigms I, II and III (paradigm IV).
Amino AcidLongevity within Paradigm I
(Paradigm IVa)
Longevity within Paradigm II
(Paradigm IVb)
Longevity within Paradigm III (Paradigm IVc)
rprprp
A−0.5985 × 10−30.2204 × 10−10.3701 × 10−1
C−0.1665 × 10−10.5232 × 10−2−0.8882 × 10−7
D0.1595 × 10−10.2204 × 10−10.1904 × 10−1
E0.0648 × 10−1−0.5052 × 10−20.2154 × 10−1
F0.3651 × 10−1−0.1755 × 10−1−0.6512 × 10−3
G−0.6552 × 10−3−0.0827 × 10−1−0.4365 × 10−2
H−0.0967 × 10−10.3501 × 10−10.5352 × 10−2
I0.5936 × 10−30.6164 × 10−30.5022 × 10−2
K0.3491 × 10−1−0.3751 × 10−1−0.1316 × 10−1
L−0.4933 × 10−2−0.6363 × 10−3−0.5541 × 10−2
M0.1565 × 10−10.4853 × 10−2−0.1635 × 10−1
N0.7124 × 10−40.3531 × 10−10.0728 × 10−1
P−0.0189 × 10−10.0209 × 10−10.6691 × 10−3
Q0.6712 × 10−3−0.1545 × 10−10.6868 × 10−4
R0.0997 × 10−10.0538 × 10−10.3212 × 10−1
S0.5102 × 10−2−0.3691 × 10−1−0.4843 × 10−2
T−0.0121 × 10−00.4903 × 10−20.6976 × 10−4
V−0.3521 × 10−10.2194 × 10−1−0.5621 × 10−2
W0.2363 × 10−10.2413 × 10−1−0.4743 × 10−2
Y−0.2573 × 10−10.2233 × 10−1−0.1615 × 10−1
p Values derive from linear correlations of double log-transformed data. Red color highlights the 4 strongest correlations among the 20 amino acids and, if applicable, their significance.
Table 4. Physiological analyses on paradigm IVc: correlation of longevity with amino acid usage in respiratory chain proteins versus all mitochondrial proteins.
Table 4. Physiological analyses on paradigm IVc: correlation of longevity with amino acid usage in respiratory chain proteins versus all mitochondrial proteins.
Amino AcidRawBody Mass-CorrectedPhylogeny-CorrectedEffect Size
rpq [%]rprpRatiop
A0.3701 × 10−11−0.0997 × 10−10.0638 × 10−11.222 × 10−2
C−0.8882 × 10−743−0.1625 × 10−1−0.7916 × 10−50.533 × 10−4
D0.1904 × 10−11−0.1505 × 10−1−0.1615 × 10−11.005 × 10−1
E0.2154 × 10−19−0.1316 × 10−10.2812 × 10−11.042 × 10−1
F−0.6512 × 10−330.2343 × 10−1−0.1705 × 10−10.742 × 10−2
G−0.4365 × 10−250.0847 × 10−1−0.4346 × 10−20.843 × 10−1
H0.5352 × 10−24−0.2862 × 10−1−0.0111 × 10−01.173 × 10−2
I0.5022 × 10−24−0.0209 × 10−10.6195 × 10−31.063 × 10−1
K−0.1316 × 10−14−0.2214 × 10−1−0.3751 × 10−10.925 × 10−1
L−0.5541 × 10−200.1047 × 10−1−0.5192 × 10−20.975 × 10−1
M−0.1635 × 10−13−0.2852 × 10−1−0.3791 × 10−10.882 × 10−1
N0.0728 × 10−112−0.5202 × 10−2−0.1416 × 10−11.046 × 10−1
P0.6691 × 10−310.0598 × 10−10.4128 × 10−21.223 × 10−2
Q0.6868 × 10−42−0.0817 × 10−10.2922 × 10−11.271 × 10−2
R0.3212 × 10−110.3252 × 10−10.4734 × 10−21.026 × 10−1
S−0.4843 × 10−21−0.0071 × 10−00.1775 × 10−10.934 × 10−1
T0.6976 × 10−40−0.2473 × 10−10.1456 × 10−11.527 × 10−3
V−0.5621 × 10−240.0199 × 10−1−0.5771 × 10−20.881 × 10−1
W−0.4743 × 10−220.1336 × 10−1−0.0289 × 10−10.917 × 10−2
Y−0.1615 × 10−10−0.3691 × 10−1−0.5641 × 10−20.988 × 10−1
Body mass0.9332 × 10−90.8865 × 10−7
p Values derive from linear correlations of double log-transformed data. For ratios, p values derive from ANOVA on ranks analyses. Red color highlights the 4 strongest correlations or the 4 largest effects among the 20 amino acids and, if applicable, their significance.
Table 5. Mitochondrially encoded amino acid usage and longevity in 218 species.
Table 5. Mitochondrially encoded amino acid usage and longevity in 218 species.
Amino AcidRawBody Mass-CorrectedPhylogeny-Corrected
rpq [%]rprp
A0.5695 × 10−2000.1306 × 10−2−0.0228 × 10−1
C−0.7372 × 10−3839−0.3222 × 10−6−0.3251 × 10−6
D0.0396 × 10−11−0.0743 × 10−10.0882 × 10−1
E0.5164× 10−1650.3874 × 10−90.1198 × 10−2
F−0.7405 × 10−391−0.2261 × 10−3−0.1335 × 10−2
G−0.0366 × 10−110.0059 × 10−10.0041 × 10−0
H0.7064× 10−3480.2502 × 10−40.2013 × 10−3
I−0.3472 × 10−70−0.1296 × 10−20.0574 × 10−1
K−0.2161 × 10−370.1493 × 10−20.0594 × 10−1
L0.4162 × 10−1000.0396 × 10−10.0336 × 10−1
M−0.2552 × 10−41−0.0982 × 10−1−0.0079 × 10−1
N−0.3512 × 10−77−0.3094 × 10−6−0.1041 × 10−1
P0.7681 × 10−43170.3502 × 10−70.1632 × 10−2
Q0.6102 × 10−2300.0962 × 10−10.0059 × 10−1
R0.4752 × 10−1340.0773 × 10−1−0.0524 × 10−1
S−0.6202 × 10−240−0.1051 × 10−1−0.0366 × 10−1
T0.7474 × 10−4010.2365 × 10−40.1671 × 10−2
V−0.1837 × 10−31−0.0227 × 10−10.0247 × 10−1
W0.3297 × 10−700.1111 × 10−10.1721 × 10−2
Y−0.5183 × 10−169−0.3933 × 10−9−0.0643 × 10−1
Body mass0.8682 × 10−67---0.5556 × 10−19
p Values derive from linear correlations of double log-transformed data. Red color highlights the 4 strongest correlations or the 4 largest effects among the 20 amino acids and, if applicable, their significance.
Table 6. Respiratory chain amino acid usage in aerobic versus facultatively anaerobic helminths (paradigm V).
Table 6. Respiratory chain amino acid usage in aerobic versus facultatively anaerobic helminths (paradigm V).
Amino AcidAerobicity (Paradigm V)
Ratiop
A1.776 × 10−5
C0.439 × 10−10
D0.841 × 10−5
E0.996 × 10−1
F0.828 × 10−5
G0.838 × 10−6
H1.354 × 10−11
I1.161 × 10−3
K1.395 × 10−5
L1.016 × 10−1
M1.512 × 10−6
N1.161 × 10−2
P1.635 × 10−9
Q1.779 × 10−12
R1.042 × 10−1
S0.935 × 10−5
T1.681 × 10−9
V0.583 × 10−12
W1.036 × 10−1
Y0.682 × 10−11
p Values derived from ANOVA on ranks analyses. Red color highlights the 4 most pronounced effect sizes among the 20 amino acids and their significance.
Table 7. Amino acid usage in aerobic versus obligately anaerobic bacteria and archaea.
Table 7. Amino acid usage in aerobic versus obligately anaerobic bacteria and archaea.
Amino AcidArchaea
(Paradigm VI)
Gram-Positive Bacteria
(Paradigm VII)
Gram-Negative Bacteria
(Paradigm VIII)
RatiopRatiopRatiop
A1.315 × 10−21.066 × 10−11.112 × 10−1
C0.732 × 10−20.531 × 10−60.947 × 10−1
D1.234 × 10−10.962 × 10−11.019 × 10−1
E0.942 × 10−10.991 × 10+00.875 × 10−3
F0.912 × 10−10.999 × 10−10.951 × 10−1
G1.102 × 10−21.046 × 10−10.997 × 10−1
H1.073 × 10−11.061 × 10−11.118 × 10−3
I0.772 × 10−20.926 × 10−10.901 × 10−1
K0.603 × 10−30.895 × 10−10.853 × 10−1
L0.973 × 10−11.064 × 10−31.031 × 10−1
M0.907 × 10−20.872 × 10−20.964 × 10−1
N0.773 × 10−20.946 × 10−11.028 × 10−1
P1.098 × 10−21.106 × 10−11.046 × 10−1
Q1.072 × 10−11.055 × 10−11.338 × 10−5
R1.243 × 10−21.067 × 10−10.998 × 10−1
S1.038 × 10−10.941 × 10−11.001 × 10+0
T1.165 × 10−21.071 × 10−21.024 × 10−1
V1.122 × 10−21.043 × 10−11.018 × 10−1
W1.066 × 10−11.113 × 10−11.182 × 10−3
Y0.925 × 10−10.914 × 10−10.914 × 10−1
p Values derived from ANOVA on ranks analyses. Red color highlights the 4 largest effect sizes among the 20 amino acids and, if applicable, their significance.
Table 8. Ranking of amino acid usage alterations in five paradigms of oxidative stress in animals.
Table 8. Ranking of amino acid usage alterations in five paradigms of oxidative stress in animals.
Amino AcidPeroxisome (I)Mitochondrion (II)Respiratory Chain (III)Longevity (IV)Aerobicity (V)
A145752 *
C #2 *1 *1 *1 *1 *
D1013122013
E91181619
F712153 *11
G13710712
H17911810
I1 *6191415
K191817129
L1616201720
M #83 *2 *107
N1884 *1514
P #3 *4 *1666
Q #614184 *3 *
R15203 *1817
S #4 *2 *141316
T201992 *5
V510594 *
W121761118
Y111513198
Ranks were determined by sorting the averaged amino acid usage ratios of all 20 animal species (86 species in paradigm V) by decreasing distance to unity (factor 1, indicating unchanged amino acid usage), as individually determined in Table 2, Table 4 and Table 6. Distances higher than 1 or lower than 1 were treated equally, but the former were labeled in green (denoting increased amino acid usage), whereas the latter were labeled in red (denoting decreased amino acid usage). Aspects of statistical significance were not considered in this merely relative ranking. Asterisks indicate that a specific amino acid was among the four most altered amino acids within a paradigm. Hash signs in the left column highlight amino acids that belonged to the latter asterisk-labeled group in at least two paradigms, namely the amino acids C, M, P, Q, and S.
Table 9. Ranking of amino acid usage alterations in bacteria and archaea.
Table 9. Ranking of amino acid usage alterations in bacteria and archaea.
Amino AcidArchaea (VI)Gram-Positive Bacteria (VII)Gram-Negative Bacteria (VIII)
A3 *117
C #2 *1 *9
D71716
E17204 *
F121910
G111818
H1595
I4 *76
K #1 *3 *3 *
L191413
M102 *12
N51314
P13511
Q16151 *
R61017
S201220
T8815
V91619
W #184 *2 *
Y1468
Ranks were determined by sorting the averaged usage ratios of all sampled prokaryotes by decreasing distance to unity (factor 1, indicating unchanged amino acid usage), as individually determined in Table 7. Asterisks indicate that a specific amino acid was among the four most altered amino acids within a paradigm. Hash signs in the left column highlight amino acids that belonged to the latter asterisk-labeled group in at least two paradigms, namely the amino acids C, K, and W. Further details are provided in the legend of Table 8.
Table 10. Frequency changes of the four primary redox-active amino acids in all investigated paradigms (in %).
Table 10. Frequency changes of the four primary redox-active amino acids in all investigated paradigms (in %).
Amino AcidPerox (I)Mito (II)RC (III)Long (IV) Aero (V)Arch (VI)GPB (VII)GNB (VIII)
C−14−26−59−47−57−27−47−6
M+11+16+114−12+51−10−13−4
W+10+5+81−9+3+6+11+18
Y+10+6+33−2−32−8−9−9
Abbreviations denote: Perox, peroxisome; Mito, mitochondrion; RC, respiratory chain; Long, longevity; Aero, aerobicity; Arch, archaea; GPB, Gram-positive bacteria; GNB, Gram-negative bacteria.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Schindeldecker, M.; Moosmann, B. Cysteine Is the Only Universally Affected and Disfavored Proteomic Amino Acid under Oxidative Conditions in Animals. Antioxidants 2024, 13, 267. https://doi.org/10.3390/antiox13030267

AMA Style

Schindeldecker M, Moosmann B. Cysteine Is the Only Universally Affected and Disfavored Proteomic Amino Acid under Oxidative Conditions in Animals. Antioxidants. 2024; 13(3):267. https://doi.org/10.3390/antiox13030267

Chicago/Turabian Style

Schindeldecker, Mario, and Bernd Moosmann. 2024. "Cysteine Is the Only Universally Affected and Disfavored Proteomic Amino Acid under Oxidative Conditions in Animals" Antioxidants 13, no. 3: 267. https://doi.org/10.3390/antiox13030267

APA Style

Schindeldecker, M., & Moosmann, B. (2024). Cysteine Is the Only Universally Affected and Disfavored Proteomic Amino Acid under Oxidative Conditions in Animals. Antioxidants, 13(3), 267. https://doi.org/10.3390/antiox13030267

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop