MINI REVIEW article

Front. Cell Dev. Biol., 04 June 2020

Sec. Cell Death and Survival

Volume 8 - 2020 | https://doi.org/10.3389/fcell.2020.00397

Cysteine Oxidations in Mitochondrial Membrane Proteins: The Case of VDAC Isoforms in Mammals

  • 1. Section of Molecular Biology, Department of Biological, Geological and Environmental Sciences, University of Catania, Catania, Italy

  • 2. Section of Biology and Genetics, Department of Biomedical and Biotechnological Sciences, University of Catania, Catania, Italy

  • 3. Organic Mass Spectrometry Laboratory, Department of Chemical Sciences, University of Catania, Catania, Italy

Abstract

Cysteine residues are reactive amino acids that can undergo several modifications driven by redox reagents. Mitochondria are the source of an abundant production of radical species, and it is surprising that such a large availability of highly reactive chemicals is compatible with viable and active organelles, needed for the cell functions. In this work, we review the results highlighting the modifications of cysteines in the most abundant proteins of the outer mitochondrial membrane (OMM), that is, the voltage-dependent anion selective channel (VDAC) isoforms. This interesting protein family carries several cysteines exposed to the oxidative intermembrane space (IMS). Through mass spectrometry (MS) analysis, cysteine posttranslational modifications (PTMs) were precisely determined, and it was discovered that such cysteines can be subject to several oxidization degrees, ranging from the disulfide bridge to the most oxidized, the sulfonic acid, one. The large spectra of VDAC cysteine oxidations, which is unique for OMM proteins, indicate that they have both a regulative function and a buffering capacity able to counteract excess of mitochondrial reactive oxygen species (ROS) load. The consequence of these peculiar cysteine PTMs is discussed.

Genome–Proteome Discordance: an Overview on Protein Posttranslational Modifications

The mammalian proteome is vastly more complex than the related genome: an imbalance between one million different proteins against approximately 25,000 genes is estimated in humans (International Human Genome Sequencing Consortium, 2004; Clamp et al., 2007). The reasons for this inconsistency reside both in the molecular mechanisms that allow a single gene to encode for numerous proteins (i.e., alternative splicing, genomic recombination, alternative promoters, and termination sites) and in the posttranslational modifications (PTMs) which represent a wide range of chemical changes that proteins may undergo after synthesis. They include the specific cleavage of protein precursors, the covalent addition or removal of low-molecular groups, and the formation of disulfide bonds or other redox modifications (Wang et al., 2014; Duan and Walther, 2015). PTMs are crucial for several cellular processes such as protein turnover and signaling and are commonly mediated by enzymatic activity. Many of them have been extensively characterized, and phosphorylation is perhaps the most known example. Conversely, only in the last two decades have reactive oxygen species (ROS)/reactive nitrogen species (RNS) been identified as physiological regulators of intracellular signaling pathways through the covalent modification of specific cysteine residues within redox-sensitive proteins (Chung et al., 2013). Although the “oxidative posttranslational modifications” (Ox-PTMs) still represent a little explored field, the increasing number of tools aimed at identifying and quantifying them [e.g., high-throughput mass spectrometry (MS) analysis] continues to broaden the knowledge of redox regulation (Spickett and Pitt, 2012; Shakir et al., 2017; Anjo et al., 2019).

Peculiarity of Cysteine Ox-PTMs and Their Correlation With Pathological Conditions

Because of their redox-reactive thiol (–SH) side chain, cysteine residues are likely subjected to various Ox-PTMs including S-nitrosylation (or S-nitrosation, SNO), sulfhydration (SSH), S-acylation, S-glutathionylation, disulfide bonds (RS-SR), sulfenylation (SOH), sulfinic acid (SO2H), and sulfonic acid (SO3H) (Forrester and Stamler, 2007; Murray and Van Eyk, 2012; Alcock et al., 2018). Except for sulfonic acid, all the reported Ox-PTMs are readily reversible and ruled by specific enzymatic activities. Sulfiredoxin (Srx1), for example, acts on oxidative states up to sulfinic acid by reducing them back to thiol in an ATP-dependent manner (Biteau et al., 2003; Chang et al., 2004). The biological significance of most of these reversible modifications has been amply investigated under physiological and non-physiological conditions (Bechtel and Weerapana, 2017). On the contrary, knowledge about the impact of irreversible Ox-PTMs on cell physiology is quite restricted.

Reversible Thiol Modification: Addition of Groups to Cysteine

S-nitrosylation, the covalent attachment of a nitrous (NO) moiety to –SH functional groups, protects proteins against further cellular oxidative and nitrosative stress. Specific enzymes, called nitrosylases, are responsible for NO group transfer in either metal-to-cysteine or cysteine-to-cysteine mechanisms (Anand and Stamler, 2012). Treatments of mouse hearts with NO donors to increase S-nitrosylation, especially in mitochondria, prior to an ischemic insult, reduce indeed infarct size by avoiding critical cysteine irreversible oxidation (Duranski et al., 2005). Under physiological conditions, it represents an important modulator of signal transduction pathways: mitochondrial SNO proteins inhibit respiratory complex I (CI) to modulate mitochondrial ROS production, promote the selective import of mitochondrial protein, enhance mitochondrial fission, and affect the mitochondrial permeability transition pore (MPTP) opening (Piantadosi, 2011; Fernando et al., 2019). However, aging or environmental toxins that increase NO production lead to aberrant S-nitrosylation reactions that contribute to neurodegenerative diseases, including Alzheimer’s disease (AD) and Parkinson’s disease (PD) (Nakamura et al., 2013).

Sulfhydration (SSH) has lately been acknowledged as a PTM analogous to nitrosylation that consists in the conversion of a –SH group to a –SSH or a persulfide group. Hydrogen sulfide (H2S) represents a ubiquitous gaseous signaling molecule with important physiological vasorelaxant properties (Iciek et al., 2015; Zhang et al., 2017) that, in mammals, is enzymatically generated by three enzymes: cystathionine β-synthase (CBS), cystathionine γ-lyase (CTH or CSE), and 3-mercaptopyruvate sulfurtransferase (3MST) (Rose et al., 2017). The reason that proteins undergo this type of modification is not known, although they have been identified by LC–MS/MS. Recently, Tonks et al. suggested that H2S, produced by CSE consequently to endoplasmic reticulum (ER) stress, sulfhydrates protein tyrosine phosphatase 1B (PTP1B) (Krishnan et al., 2011). This event, in turn, causes the ER kinase PERK activation during the response to ER stress. Interestingly, anomalous sulfhydration has been linked to several pathological conditions ranging from heart diseases to neurodegeneration (i.e., PD) (Paul and Snyder, 2018).

S-acylation is a highly conserved PTM that takes place in all eukaryotic organisms and is regulated by the same enzyme families from yeast to humans. It consists in the covalent attachment of an acyl chain to a cysteine residue and is the only fully reversible posttranslational lipid modification of proteins. Because of the weak nature of the thioester bonds within the intracellular environment, most S-acylated proteins rapidly undergo S-acylation and deacylation cycles (Chamberlain and Shipston, 2015; Feldman et al., 2015; Pedro et al., 2017). S-acylation increases protein hydrophobicity that, in turn, can affect several properties ranging from structure to assembly, maturation, and function. Although the major lipid incorporated into endogenous proteins is palmitate (C 16:0), fatty acids such as oleate (C 18:1) and stearate (C 18:0) can also be added. In particular, S-acylation is detected in transmembrane proteins where it is required for stable membrane binding. In this regard, the blockage in the S-acylation of specific residues increases protein ubiquitination and degradation (Valdez Taubas and Pelham, 2005). The recent development of S-acylation proteomic profiling together with the identification of enzymes that regulate protein S-acylation and deacylation has brought new interest into the physiological function of this peculiar Ox-PTM. To date, 47 acylated proteins have been detected in the yeast Saccharomyces cerevisiae (Roth et al., 2006) whereas mammalian cells contain hundreds of these modified proteins (Wan et al., 2007; Yang et al., 2010). Interestingly, there is no evidence for S-acylation occurrence in prokaryotes. On the contrary, S-acylation of many viral proteins is catalyzed by the host cell machinery (Kordyukova et al., 2010).

S-glutathionylation (SSG) consists in the addition of the tripeptide glutathione (GSH), the main low-molecular-weight antioxidant of both prokaryotes and eukaryotes, to protein cysteine residues through the establishment of a covalent linkage (Dalle-Donne et al., 2009). This reversible thiol modification is promoted by oxidative and/or nitrosative stress and acts as a repository for reduced glutathione, since the oxidized form (GSSG) is either rapidly excreted from the cells or reduced back to GSH via NADPH-dependent glutathione reductase (Lushchak, 2012). Anyhow, S-glutathionylation also occurs under unstressed conditions (Dalle-Donne et al., 2009). Not all cysteines can undergo S-glutathionylation: it has been proposed that low pKa values and, possibly, the three-dimensional proximity to His, Lys, and Arg residues are key factors that make specific Cys appropriate targets for such PTMs (Grek et al., 2013). Numerous molecular mechanisms have been suggested to explain protein S-glutathionylation, ranging from the direct thiol–disulfide exchange reaction of a cysteinyl residue with GSSG to the activation of a thiol group to reactive sulfur intermediates (sulfenyl amide, sulfenic acid, or S-nitrosyl) which then react with GSH (Allen and Mieyal, 2012). The enzymes responsible for GSH conjugation are called glutathione S-transferases, while the reverse reaction is catalyzed by glutaredoxins (Grxs). Certainly, the main role of S-glutathionylation is to prevent thiol irreversible oxidation, though it can affect protein function by altering charge and molecular mass, by inducing conformational changes, or by sterically blocking catalytic sites (Grek et al., 2013). The synergy between S-glutathionylation and S-deglutathionylation thus contributes to the regulation of redox homeostasis, cytoskeletal assembly, protein folding, and stability. Moreover, several evidences have reported their involvement in the control of cell signaling pathways associated with viral infections and apoptosis induced by tumor necrosis factor (Dalle-Donne et al., 2009). Mitochondria contain several S-glutathionylation targets thanks to the thiol transferase glutaredoxin 2 (Grx2), which is considered an important redox sensor within the organelle (Xiong et al., 2011).

Reversible Thiol Modification: Redox Modification of Cysteine Sulfur

Disulfide bonds (RS-SR) are essential PTMs involved in protein folding and in the stabilization of their tertiary and quaternary structures. Disulfide formation depends on the spatial proximity to another cysteine and can also occur through a reaction with sulfenic acid in the presence of high concentrations of ROS. They indeed convert SOH groups into thiol radicals (RS⋅) which, in turn, react with other thiolates to form a disulfide bond (reviewed in Summa et al., 2007; Depuydt et al., 2011; Herrmann and Riemer, 2012). In eukaryotic cells, specific enzymes catalyze disulfide exchange within the ER and the mitochondrial intermembrane space (IMS). In yeast, for example, the ER contains the sulfhydryl oxidase endoplasmic oxidoreductin 1 (Ero1) that transfers oxidizing equivalents first to the disulfide isomerase (PDI) and then to secretory proteins. Disulfide formation in the IMS of mitochondria is entrusted, instead, to the MIA system, which will be described in depth in the following sections (Stojanovski et al., 2008). Noteworthy is that some RS-SR are dynamic and modulate protein signaling: that is, in mammalian cells, the stress sensor NPGPx transmits oxidative stress signals through its Cys57–Cys86 disulfide bond. Once oxidized, NPGPx binds to GRP78 (glucose-regulated protein), giving rise to covalent intermediates between NPGPx-Cys86 and GRP78-Cys41/Cys420 that culminate in the enhancement of the GRP78 chaperone activity. The knockout (KO) of NPGPx gene increases the intracellular ROS content and impairs GRP78 chaperone activity, leading to the accumulation of misfolded proteins (Wei et al., 2012). An involvement of S–S bridges in neurodegenerative disorders has been recently proposed. Some evidence suggests that intermolecular disulfide bonds in mutant superoxide dismutase 1 (SOD1) may have a role in its aggregation, a hallmark of some familial amyotrophic lateral sclerosis (ALS) variants (Furukawa et al., 2006). Moreover, intramolecular and intermolecular disulfide bonds seem to be involved in the pathological Tau isoforms typical of AD and other neurodegenerative disorders (Walker et al., 2012; Kim et al., 2015).

Sulfenylation (RSOH) is the first step of cysteinyl thiol oxidation. Because of their high reactivity and instability, sulfenic acids (–SOH) have been reputed for a long time as intermediates to additional cysteine modifications (Lo Conte and Carroll, 2013). They can be generated following the condensation of a thiol with a hydroxyl radical, the hydrolysis of S-nitrosothiols, and the reaction of a thiol or a thiolate anion with a high concentration of hydrogen peroxide and with thiosulfinates (Hall et al., 2010). Moreover, sulfenic acids can result from the reaction of a thiol or thiolate anion with hypochlorous acid, which represents a powerful bactericidal compound generated within neutrophils during inflammation. The reactivity of the thiol side chain is modulated by the surrounding protein environment that drives the fate of sulfenic acid to an irreversible chemical modification or to a “protective” disulfide bond (Ferrer-Sueta et al., 2011; Roos and Messens, 2011). Accumulating evidence proposes cysteine S-sulfenylation as a redox-based signal transduction mechanism in living cells, which is also involved in transcription regulation (Poole and Nelson, 2008; Roos and Messens, 2011; Heppner et al., 2018). The functional domains of several protein phosphatases, acetyltransferases, kinases, deubiquitinases, and deacetylases contain S-sulfenylated sites, suggesting the existence of a regulatory cross-talk between S-sulfenylation and other major PTM events. Tandem MS revealed –SOH modification also in the catalytic cycle of many enzymes as well as peroxiredoxin (Prx) and NADH peroxidase (Poole et al., 2004; Portillo et al., 2017). Dysregulated protein sulfenylation has been associated with several human pathologies including hypertension (Camargo et al., 2018) and aggressive cancer phenotypes (Truong et al., 2016).

Sulfinic acid (RSO2H) is the further oxidized step of –SOH achieved in the presence of excess oxidants (Griesser et al., 2018; Urmey and Zondlo, 2020). The biological function of sulfinic acid has been largely debated, due to the general belief that this Ox-PTM could be an artifact triggered by protein isolation procedures. However, this idea was totally discredited after the discovery that the catalytic cysteine of human Prx1 is selectively oxidized to Cys-SO2H but not to Cys-SO3H, as sulfinic acid is re-reduced to Cys-SH by sulfiredoxin (SRX). In mammals, intracellular sulfinic acid levels depend upon cysteine dioxygenase (CDO), a metabolic enzyme that regulates cysteine homeostasis according to the dietary intakes of sulfur amino acids (Stipanuk et al., 2009). In the last decade, a functional role for –SO2H modification has been reported for a growing list of proteins (Lo Conte and Carroll, 2013), including the copper–zinc SOD1, the PD protein DJ-1, and the D-amino acid oxidase (DAO) (Blackinton et al., 2009). Sulfinic acid modulates indeed protein metal binding features: in matrix metalloproteases, the oxidation of a cysteine residue to –SO2H activates protease function and is essential for catalytic activity. Despite its association with diseases linked to oxidative stress, including cancer and neurodegenerative disorders, the rationale of sulfinic acid is essentially unknown because of the difficulties of its detection (Neumann et al., 2003; Canet-Avilés et al., 2004). Antibodies against sulfinic acid have limited affinity and specificity. Moreover, MS analysis struggles with the existence of other modifications with the same nominal mass shift, like S-sulfhydration (–SSH) (Akter et al., 2018).

Irreversible Thiol Modification: A Way With No Return

Sulfonic acid (RSO3H) results from the oxidation of –SO2H by strong oxidizers such as ONOOH, H2O2, or HOX, and it is considered a dead-end product since no biological pathway has been discovered for –SO3H reduction. Due to their relative stability compared to the reversible thiol modifications, they can be detected directly by MS, although with some difficulties. Sulfinic and sulfonic acid fragmentations are indeed nearly indistinguishable from the fragmentation patterns of other PTMs and necessitate high-resolution MS instruments to differentiate them (Gaillot et al., 2020). To date, the role of cysteine sulfonic acid in protein function/conformation is still controversial since it is not completely clarified if Cys-SO3H can be counted as an oxidative damage or a signaling modification. Literature contains just a few examples in this regard. Lim et al. (2008) proposed that the sulfonic acid formation of yeast Prx1 cysteine in an active site enhances the protein chaperone activity, representing a marker of the cumulative oxidative stress in cells. Lately, Wu et al. (2017) refined this concept, concluding that the two indispensable cysteines of human Prx1, cysteine in the active site (CP), responsible for peroxide reduction, and resolving cysteine (CR), which deals with CP regeneration, exhibit different propensities to sulfonation. Accordingly, the most frequent overoxidation of human Prx1 hits CR: sulfonation of Prx1 fosters its activity as a chaperone, preventing proteins from unfolding and irreversible aggregation (Wu et al., 2017). The existence of cysteine residues susceptible to sulfonation appears as an exclusive evolutionary feature of eukaryotic Prxs involved in cell survival during oxidative stress and argues in favor of Cys-SO3H as a signaling PTM rather than a simple oxidative damage (Delaunay et al., 2002; Jang et al., 2004; Jarvis et al., 2012; Sobotta et al., 2015). Interestingly, overoxidation of Prxs cysteines follows a precise circadian rhythm in human blood cells that has been identified as a transcription-independent marker (O’Neill and Reddy, 2011; Edgar et al., 2012). The Cys149 residue of the active site of GAPDH is also rapidly overoxidized to sulfonic acid by H2O2, and this modification is associated with glycolysis inhibition and with numerous side activities including the induction of apoptosis (Colussi et al., 2000). Cys149 irreversible oxidation is regulated by S-glutathionylation and disulfide bond formation with the residue Cys154 (Muronetz et al., 2020). In human mutant SOD1, the peroxidation of Cys111 to sulfonic acid has been implied in the pathology of the degeneration/death of familial ALS motor neurons. Immunohistochemical analysis proved indeed the accumulation of mutant SOD1 with overoxidized Cys111 residue in the Lewy body-like hyaline inclusions and vacuole rims of the spinal cord of familial ALS-transgenic mice (Fujiwara et al., 2007). A pathophysiological relevance of sulfonation has been observed also in PD: an increase in the level of sulfonated parkin correlates with its insolubility in human PD brains (Chakraborty et al., 2017). In endocrinology, sulfonation represents a particularly explored biological phenomenon (Wallace et al., 2010). Cys-SO3H modification of neuroendocrine peptides strongly affects their receptor binding (Wallace et al., 2010). The oxidation of an N-terminal cysteine residue to sulfonic acid, such as the Cys2 residue of GTPase-activating proteins, has been described as a marker for ubiquitin-dependent protein degradation (Arg/N–degron pathway) (Barros and McStay, 2020). The acetylation of the N-terminal amino group also represents a protective mechanism against irreversible oxidation and seems to act as a degradation signal in the Ac/N–degron pathway. Sulfonic acid is present in mammalian cells as taurine, a low-molecular-weight antioxidant, also an osmoregulator and participant in calcium signaling pathways. Taurine is particularly abundant in the heart, skeletal muscle tissue, retinas, and central nervous system.

Mitochondrial Proteins PTMs Involving Cysteine Oxidation

Cysteine residues constitute only 2% of the total amino acid content, but the mitochondrial proteome is rich in protein thiols. Mitochondria are dynamic organelles responsible for maintaining redox cellular homeostasis thanks to the continuous generation and wasting of ROS/RNS. The production of ATP and the formation of the proton gradient are, however, imperfect, as the premature leak of electrons from complexes I and III leads to the formation of O2. MnSOD rapidly dismutates O2 to hydrogen peroxide more stably than other ROSs. The electron transport chain, in addition to the nitrite reductase activity, has been proposed as responsible for the NO⋅ generation within mitochondria: the existence of a mitochondrial nitric oxide synthase (mtNOS) is, indeed, still debated. H2O2 is the most relevant radical species in redox signaling that reversibly oxidizes protein cysteine thiols (–SH) to sulfenic acid (–SOH) through a mechanism called “redox switch” (Spadaro et al., 2010). However, O2 also exerts a signaling role by disassembling Fe–S clusters in several enzymes of the TCA cycle and in respiratory complexes (Mailloux et al., 2014). At high concentrations, O2 and H2O2 can lead to oxidative stress due to the production of the hydroxyl radical (⋅OH), the ROS that is most reactive. In addition, O2 can combine with nitric oxide (NO) to form the powerful oxidant peroxynitrite (ONOO) (Fridovich, 1995; Radi, 2013). Under physiological conditions, ROS/RNS generated in mitochondria mainly regulate mitochondrial protein interactions, localization, and activity (Ray et al., 2012; Mailloux et al., 2014) through the oxidative modification of redox-sensitive cysteine residues. It has been reported that the sulfenylation of the Cys253 residue of uncoupling protein 1 (UCP1) within the brown adipose mitochondria of mammalian cells participates in the thermogenic regulation of energy expenditure (Chouchani et al., 2016). In vivo in mice, the selective S-nitrosation of the Cys39 residue on the ND3 subunit of CI was instead proposed to be cardioprotective (Chouchani et al., 2013). If unregulated, ROS can oxidatively damage biomolecules (proteins, DNA, and lipids), ultimately leading to neurological disorders, metabolic diseases, aging (Bulteau et al., 2006), and other pathologies. Hence, the control of the cellular ROS balance mainly relies on two redox regulatory systems: the reduced/oxidized thioredoxin and the glutathione (GSH)/glutathione disulfide (GSSG) (Holmgren et al., 2005; Jones, 2006).

Mitochondrial Compartments Possess Different Redox Potentials

The outer mitochondrial membrane (OMM), non-specifically permeable to low-molecular-weight solutes, and the inner mitochondrial membrane (IMM), selectively permeable, define two distinct aqueous compartments within the organelle: the matrix and the IMS. The mitochondrial matrix, that is, the innermost compartment, exhibits a pH of about 8.0 (Llopis et al., 1998) and is the site of numerous enzymatic reactions including organellar DNA replication, transcription, and protein biosynthesis. Components of the mitochondrial IMS are involved, instead, in several processes ranging from the transport of proteins, metal ions, or electrons to the inner membrane proteins assembling up to cellular respiration. The IMS also sequesters several apoptotic components released into the cytosol to trigger programmed cell death. The IMS pH value is about 0.8 units more acidic than the matrix: this asymmetry creates an electrochemical gradient across the inner membrane, exploited in ATP synthesis (Holoubek et al., 2007). Interestingly, the pH of the matrix and IMS differently affects their redox potentials. ROS production is strongly influenced by pH (Shu et al., 1997): in hepatocytes, the extracellular pH controls hydroxyl radical (⋅OH) and NO production (Shu et al., 1997). Likewise, the more acidic pH of the mitochondrial IMS makes it a strong oxidizing environment (Hu et al., 2008; Kojer et al., 2012, 2015). Besides, the IMS of eukaryotic cells contains a dedicated machinery, called the MIA/CHCHD4 complex or disulfide relay-dependent import machinery, responsible for the import and oxidative folding of substrates with conserved cysteine motifs (CX3C and CX9C). The oxidoreductase Mia40 and the FAD-dependent sulfhydryl oxidase essential for respiration and viability 1 (Erv1) are the crucial components of the MIA complex in yeast. Initially, Mia40 forms a transient intermolecular disulfide bond with imported precursor proteins. Afterward, the isomerization of this disulfide bond causes the introduction of an intramolecular S–S bridge within the imported precursor protein. Following oxidation of the substrate, Mia40 is reoxidized by Erv1 (Lee et al., 2000; Dabir et al., 2007; Banci et al., 2009; Tienson et al., 2009; Bien et al., 2010). A stepwise evolution of the MIA pathway has been proposed in Allen et al. (2008). Accordingly, the ancestral IMS import pathway only included an Erv1 protein, as nowadays in Leishmania and Trypanosoma. At this stage, Erv1 can still work alone, suggesting that Mia40 behaves as a kind of import chaperone. In fungi, Erv1 necessarily required Mia40 to oxidize substrate proteins because of its inability to do it by itself. The metazoan homolog of yeast Mia40 is the CHCHD4 protein, with the same role as the core component of the disulfide relay-dependent import machinery. They share the same ability to oxidize pairs of reduced cysteines to disulfide bridges, provided that they are present in the motif CX3C or CX9C. While the yeast Mia40 is a transmembrane protein, metazoan CHCHD4 is smaller, soluble, and moving inside the IMS. Since Mia/CHCHD4 are involved in the proper folding of Cys-containing proteins, they are also considered key components of the quality control system that impacts the IMS targeted proteins (Chatzi and Tokatlidis, 2013). To date, several IMS proteins have been detected with structural disulfide bridges, and the list is constantly increasing (Cavallaro, 2010; Gross et al., 2011; Vögtle et al., 2012). A few examples are the intermitochondrial space Atp23 protein, which is a strongly oxidized protease associated with the IMM (Weckbecker et al., 2012); the protein anamorsin, implicated in Fe/S cluster assembly (Banci et al., 2011); and the Ca2+ uniporters MICU1 and MICU2 (Petrungaro et al., 2015). Disulfide bonds within IMS proteins have also lately been considered as part of a quality control system: Longen et al. (2014) proposed that the presence of these Ox-PTMs in the eukaryotic mitochondrial ribosomal Mrp10 protein is sufficient to avoid its proteolytic degradation in the IMS before moving to the mitochondrial matrix. Nevertheless, the presence of special reductases (i.e., thioredoxin and glutaredoxin; Holmgren, 2000) within the IMS causes the reduction of cysteine residues in many proteins and protein domains contained therein. Most of these cysteines are coordination sites for metal ions (Sco1, Sco2, Cox2, Cox11, and Cox17), heme cofactors (cytochrome c and cytochrome c1), and iron–sulfur clusters (Rieske Fe/S protein Rip1). The mitochondrial matrix possesses several redox repair systems: reduced/oxidized thioredoxin and the GSH/GSSG are both involved in disulfide bonds or sulfenic acid reduction (Santo-Domingo et al., 2015).

Redox Modification of OMM Proteins: The VDAC Case

In spite of plenty of reports about the redox modifications of the mitochondrial matrix, IMS, and inner membrane proteins, very little is known about the outer membrane proteins. This is possibly due to underrepresentation of membrane proteins in proteome profiles (protein heterogeneity, hydrophobicity, limited solubility, restricted enzyme accessibility, and low abundance are responsible for it) that downsizes their detection and identification. In order to overcome these restrictions, numerous strategies have been developed to increase the enrichment, solubilization, separation, and cleavage of membrane proteins (Schmitt et al., 2006; Niemann et al., 2013; Bak et al., 2017). Voltage-dependent anion selective channels (VDACs) are surely the most represented proteins of the OMM involved in cellular processes that range from metabolism regulation to cell death control (Benz, 1994; Colombini et al., 1996). Alterations of VDAC expression and activity correlate with several pathologies such as cancer and neurodegenerative disorders (Shoshan-Barmatz et al., 2010; Magrì et al., 2016, 2018; Reina and De Pinto, 2017). In higher eukaryotes, VDAC exists in three different isoforms called VDAC1, VDAC2, and VDAC3 according to the chronological order of their uncovering. Evolutionary analyses suggest that VDAC3 is the oldest protein and VDAC1 is the most recent mitochondrial porin (Young et al., 2007; Messina et al., 2012). Structurally, VDAC forms an aqueous channel arranged in a transmembrane β-barrel linked to an α-helix moiety at the N-terminus (Bayrhuber et al., 2008; Hiller et al., 2008; Ujwal et al., 2008; Schredelseker et al., 2014). In each β-strand, hydrophilic and hydrophobic residues alternate regularly: the former point to the lumen of the pore while the latter interact with the membrane environment. Interestingly, all β-strands follow a steady array of antiparallel units with the exception of β-strands 1 and 19 that, running parallel, make both the N-terminus and the C-terminus point toward the IMS (Tomasello et al., 2013). The amphipathic α-helix tail seems to be placed inside the pore, but, to date, the exact position and local structure of this segment are still debated (De Pinto et al., 2007). However, although the general organization of the three isoforms is similar, different molecular weights and amino acid compositions may be associated to distinct roles in the cell. In particular, VDAC isoforms differ in the number and distribution of cysteines so that they have been the subject of specific investigations to unravel a peculiar role for these residues in protein activities (De Pinto et al., 2016; Reina et al., 2016). In humans, VDAC1, VDAC2, and VDAC3 have two, nine, and six cysteines, respectively. Considering VDAC3 is the oldest isoform, it is tempting to speculate that evolution has decreased cysteine number in VDAC1, according to an ubiquitous expression and proapoptotic function, while increasing cysteine content in VDAC2, possibly to favor its involvement as an antiapoptotic protein (De Pinto et al., 2016). VDACs are exposed to both cytosol and IMS (Figure 1). Interestingly, many of the cysteine residues are located in protein loops exposed to IMS and its oxidative potential. In addition, VDAC has been proposed as the only escape route from the mitochondria of the unreacted superoxide radical (O2–•), a by-product of complex III of the respiratory chain (Han et al., 2003). At a variance from hydrogen peroxide (H2O2), O2–• is very active but has a short life and cannot freely diffuse trough the membrane (Bleier et al., 2015). Our research group was the first to report a detailed profile of the oxidative state of VDAC cysteines and methionines (Saletti et al., 2017, 2018) by affinity chromatography purification and UHLC/high-resolution n-ESI-MS/MS. What we noticed is that a cysteine redox modification pattern is conserved during evolution. In both human (Pittalà et al., 2020) and rat VDACs (Saletti et al., 2017, 2018), indeed, residues are in the reduced form (i.e., free thiol groups available to disulfide bond formation), while others are constantly and irreversibly overoxidized to sulfonic acid. We thus hypothesized that a specific modification of each cysteine may correlate to a specific structural or functional role. Very recently, a “shot-gun” cysteine-targeted MS analysis that identified ∼1500 reactive cysteine residues on ∼450 mitochondrial proteins in HEK293T cells was able to detect cysteines highly sensitive to S-nitrosoglutathione (GSNO) also in VDAC proteins (Bak et al., 2017).

FIGURE 1

Cysteine Residues in VDAC Isoforms

VDAC1

VDAC1 is the best-characterized VDAC isoform. It provided the first description of typical voltage dependence and ion selectivity features that are the fingerprint of these pores. VDAC1 has been thoroughly characterized by electrophysiological studies (Blachly-Dyson et al., 1993; Xu et al., 1999) and is considered one of the main actors of bioenergetic metabolism and cell death regulation (Maldonado et al., 2010; Camara et al., 2011; Maldonado and Lemasters, 2012). Its 3D structure was elucidated (Bayrhuber et al., 2008; Hiller et al., 2008; Ujwal et al., 2008). A large number of reports describe VDAC1 as a docking site for proapoptotic (Bax, Bak, and Bim) and antiapoptotic (Bcl2 and Bcl-xL) factors (Shimizu et al., 2001; Vander Heiden et al., 2001; Arbel et al., 2012; Huang et al., 2013; Liu et al., 2015), as well as for cytosolic enzymes [hexokinases (HK) I and II] (Fiek et al., 1982; Abu-Hamad et al., 2008). The interaction with HK makes VDAC1 a key player in cancer and other pathologies, including neurodegenerative diseases (Maldonado and Lemasters, 2012; Leanza et al., 2014; Reina and De Pinto, 2017; Magrì et al., 2018), that share mitochondrial dysfunction and oxidative stress. VDAC1 contains only two cysteines (Cys127 and Cys232). Cys127 residue (β-strand 8) protrudes in phospholipidic hydrophobic milieu, while Cys232 is located in β-strand 16 and faces the water-accessible side of the channel (Bayrhuber et al., 2008; Hiller et al., 2008; Ujwal et al., 2008). Still, there is no evidence that oxidation can rule VDAC1 activity through the Ox-PTMs of its sulfhydryl groups. Nevertheless, already in De Pinto et al. (1991), the authors established that the two conserved cysteine residues of VDAC1 from bovine heart mitochondria are sensitive to oxidation. Both the oxidized and reduced forms of VDAC1, however, showed the same biophysical properties into artificial lipid membranes, disclaiming any redox regulation of channel activity (De Pinto et al., 1991). Later, Aram et al. (2010) confirmed that the pore is not affected by VDAC1 cysteine deletion and observed that ROS-producing agents such as selenite, As2O3, or H2O2 do not change apoptosis in human cells overexpressing rat VDAC1. Moreover, cross-linking experiments let them conclude that VDAC1 cysteines are not involved in intermolecular S–S bridge and thus do not take part in oligomerization (Aram et al., 2010). Overall, reports available so far agree that cysteine residues have no functional role in VDAC1 (Teijido et al., 2012). Recent findings revealed a conserved pattern in oxidative status, conserved among human and rat VDAC1: Cys127 was always detected in the trioxidized form of –SO3H and Cys232 exclusively in the reduced and carboxyamidomethylated form (Pittalà et al., 2020).

VDAC2

Analogous to isoform 1, VDAC2 participates in several cellular processes, from metabolite exchange through the OMM (Maurya and Mahalakshmi, 2016; Naghdi and Hajnóczky, 2016) to calcium homeostasis regulation (Subedi et al., 2011; Shimizu et al., 2015). Although ubiquitous, it is predominantly expressed in the heart and testes and, at the same time, poorly represented in the liver (Yamamoto et al., 2006; Calvo et al., 2016). The N-terminal region, which contains the voltage sensor of the channel, carries an extension of 11 additional residues only in VDAC2. The only available crystal structure is zebrafish VDAC2 (zfVDAC2), whose sequence is only about 90% similar to mammalian VDAC2 (Schredelseker et al., 2014), and confirms substantial analogy with the published structures of mouse and human VDAC1. Mouse, human, zebrafish, and bovine isoform 2 easily insert into bilayer membranes as canonical VDAC pores with similar single-channel conductance, voltage gating, and anion selectivity (Blachly-Dyson et al., 1993; Xu et al., 1999; Komarov et al., 2005; Menzel et al., 2009; Maurya and Mahalakshmi, 2015). VDAC2 is considered antiapoptotic, a distinction from VDAC1 and VDAC3: VDAC2 complexes with the proapoptotic factors BAK (Cheng et al., 2003) and BAX (Chin et al., 2018), thereby hampering programmed cell death. Accordingly, VDAC2 KO is lethal in utero (Cheng et al., 2003), and conditional KO animals display severe dysfunction in the apoptotic pathways (Ren et al., 2009). VDAC isoform 2 also possesses the highest cysteine content (nine in humans and 11 in rats, Figure 2). This intriguing hallmark has prompted numerous studies aimed at identifying the role of such residues in protein structure and function. Maurya and Mahalakshmi (2013, 2014, 2015) demonstrated that VDAC2 cysteines strengthen β-barrel association to OMM, while lowering structure stability. Moreover, it has been reported that deletion (Majumder and Fisk, 2013) or modification (Piroli et al., 2016) of human VDAC2 cysteine residues can impair channel function. For instance, VDAC2 Cys47 and Cys76 succination has been linked to reduced ATP synthesis within mitochondria of the Leigh syndrome mouse model (Piroli et al., 2016). In Naghdi et al. (2015), cysteines are considered unessential for VDAC2-mediated Bak recruitment and tBid−induced cytochrome c release. A role for cysteines in ROS regulation has also been speculated (Maurya and Mahalakshmi, 2015), although no experimental evidence has yet proven it. Saletti et al. (2018) and Pittalà et al. (2020) first supplied a detailed profile of VDAC2 Cys Ox-PTMs that reveals a precise and evolutively conserved scheme in thiol oxidation. The N-terminal cysteine residues 8 and 13, together with the unique thiol group located in a loop exposed toward the cytosol (i.e., Cys 227 in Figure 1), were reported as totally reduced. Cysteines 47, 76, 103, and 210, all located in IMS loops, were instead identified as partially oxidated to –SO3H. Surprisingly, VDAC2 Cys138, which faces the lipid environment of the OMM (Figure 1), was found fully trioxidized to sulfonic acid as the homologous residue in VDAC1 (Cys127, Figures 1, 2).

FIGURE 2

VDAC3

VDAC3 is the least known isoform. Its 3D structure has not yet been experimentally resolved, albeit several homology modeling predictions have confirmed a β-barrel organization almost identical to VDAC1 and VDAC2 (Bayrhuber et al., 2008; Hiller et al., 2008; Ujwal et al., 2008; De Pinto et al., 2010; Schredelseker et al., 2014). Unlike VDAC1 and VDAC2, ubiquitously expressed in all eukaryotes, VDAC3 prevails in the cerebral cortex, liver, heart, testes, and spermatozoa (Sampson et al., 2001). Herein, it participates in sperm mobility since VDAC3-null mice show a single microtubule doublet loss within the epididymal axoneme (Sampson et al., 2001) and in disassembly of cilia during the cell cycle by targeting Mps1 protein kinase to centrosomes (Majumder and Fisk, 2013). S. cerevisiae cells devoid of endogenous porin (Δpor1) have been extensively exploited for functional studies on VDAC proteins: as reported, the heterologous expressions of VDAC1 and VDAC2 fully complement the growth defect associated with gene loss, while that of VDAC3 only partially did so (Sampson et al., 1997; Reina et al., 2010). VDAC3 pore-forming activity was more difficult to be detected and did not show the homogeneous activity usually found with VDAC1 and VDAC2. Recombinant VDAC3 resulted in poor insertion in bilayers and hectic pore formation (Xu et al., 1999; Checchetto et al., 2014). In Checchetto et al. (2014), we demonstrated that, under non-reducing conditions, the majority of inserted VDAC3 forms very small pores with a lower conductance compared to the other isoforms. The modified activity pattern of VDAC3 pushed to look for differences between the VDAC3 structure and, in particular, VDAC1. The most striking was the cysteine content (De Pinto et al., 2016). Human VDAC3 contains indeed six cysteines (residues 2, 8, 36, 65, 122, and 229; Cys2 is the actual N-terminal residue, since Met1 is removed after biosynthesis) (Saletti et al., 2017, 2018; Pittalà et al., 2020). Except for Cys8, located within the pore, all cysteine residues protrude toward the IMS. The first clue about the relevance of VDAC3 cysteines goes back to the year 2010, when “swapping” experiments of the VDAC3 N-terminus (containing two cysteines) with the corresponding region of VDAC1 (lacking cysteines) totally changed isoform 3 pore activity (Reina et al., 2010). Next, electrophysiological data reported VDAC3 as able to form typical voltage-dependent pores (although with a remarkable lower insertion rate compared to VDAC1 and VDAC2), when refolded under highly reducing conditions (Queralt-Martín et al., 2020). Deletion of cysteines in engineered human VDAC3 molecules strongly affected both electrophysiological parameters and its ability to revert the mutant phenotype of Δpor1 yeast (Okazaki et al., 2015; Reina et al., 2016). The most relevant residue(s) for protein function appears to be 2, 8, and 122 in human VDAC3, since their mutation always restored formation of large pores and Δpor1 growth (Reina et al., 2016). Furthermore, the simultaneous mutation of any VDAC3 cysteines to alanine had a similar effect (Queralt-Martín et al., 2020). From MS analysis (Saletti et al., 2018; Pittalà et al., 2020) the preferred redox state of cysteines is conserved between rat and human VDAC3 analyzed. Accordingly, Cys2, Cys8, Cys122, and Cys229 were entirely identified in the reduced form, with the N-terminal Cys2 found acetylated. Cys36 and Cys65 were detected in a reduced form and trioxidized to sulfonic acid.

Conclusion

Hypotheses About the Role of Cys Ox-PTMs in VDAC Isoforms

Based on MS data reported for VDAC1, VDAC2, and VDAC3 in humans and rats, we found that the oxidation state of each specific residue is the same (Saletti et al., 2018; Pittalà et al., 2020). Thus, we speculated about the relationships between location of residue and oxidation state. In Figure 1, the position of the cysteines with respect to the water/phospholipidic phases is shown. Figure 2, instead, exhibits, on an analytical scheme, the oxidative modifications borne by the residues with respect to their location in structural elements of the pore. We thus noticed that cysteines whose lateral residue is embedded in the hydrophobic moiety of the membrane are trioxidized to –SO3H. This finding is very surprising, since the charged modified sulfur residues protrude into the most hydrophobic part of the membrane bilayer, an apparent oddity, if it has no reason to be, unless it has a function as a spot for a structure intended to dock the OMM. These residues are present only in VDAC1 and VDAC2.

Cysteines in the N-terminal domain sequence steadily emerge as totally reduced. To be detected in MS as carboxymethylated sulfur, they indeed have to be protected in vivo by irreversible oxidation events. VDAC2 and VDAC3 alone have cysteines in the N-terminus. This group of totally reduced cysteines, exposed to IMS or residing inside the channel, could form disulfide bonds in VDACs as proposed in Okazaki et al. (2015) and Reina et al. (2016) and with the aim of changing pore useful diameter for modulating water-soluble molecule flow. It is not known how the formation of these hypothetical disulfide bridges may happen. It is tempting to speculate that the MIA/CHCHD4 complex could be involved. Although VDACs cannot be the subject or component of this specific quality control system because of a different import pathway, it is suggested to imagine that in metazoans CHCHD4 might interact with the VDAC2 and VDAC3 cysteines protruding toward IMS in reduced form (Pittalà et al., 2020). VDAC2 shows a motif CX4C and VDAC3 a motif CX5C in the N-terminal end. The possibility that CHCHD4 can oxidize cysteine pairs in a non-conventional motif has been highlighted for mitochondrial p53, which contains CX5C and CX1C motifs (Zhuang et al., 2013; Park et al., 2016), and it is true also for Erv1, containing a CX16C motif (Terziyska et al., 2007). CHCHD4, indeed, could represent a direct link of VDAC with apoptosis-inducing factor (AIF) (reviewed in Reinhardt et al., 2020): CHCHD4 is the first AIF mitochondrial interactor that was discovered. A CHCHD4 cysteine pair is reoxidized by ALR and docked to AIF-NADH. Another intriguing connection between VDACs and AIF could be found in the essential role of AIF in import and folding of CI subunits (Vahsen et al., 2004): it is known that CI is the most important ROS producer within the IMS, and these ROS are candidate to oxidize VDAC cysteines.

A ring of residues exposed or predicted to be exposed to the IMS is the third group of cysteines with distinct modifications as found by MS spectra: it turns out that they are in either oxidized or carboxyamidomethylated form, thus inclined to variable and possible reversible oxidations. This girdle of cysteine residues located in the turns connecting the β-strands and exposed to IMS can thus work as an oxidation buffer located on the inner surface of the OMM with the function of counteracting excess of ROS produced by the OXPHOS (Reina et al., 2016).

Our interpretation of such differences, experimentally determined by MS, is that types of cysteine residues may have specific roles in the operation of VDAC isoforms.

Thus, the pattern of cysteine oxidation could also be the key to understanding the evolution of the three VDAC isoforms: VDAC1 has only two, one of them in sulfonic state and exposed to the hydrophobic layer of the OMM; VDAC2 has 9–11 cysteines, in the three categories we described above; thus, it has also the overoxidized one; VDAC3 has again six to seven cysteines but lacks the category of overoxidized one.

Okazaki et al. (2015) proposed that the transient formation of a S–S bridge between VDAC3 residues placed in the protein N-terminus with those at the bottom of the pore could change its permeability, albeit the nature of such conductance hindrance is not clear yet. Molecular dynamic simulations have shown that disulfide formation between the nearest cysteines available did not affect the channel diameter but rather changed the electrical charge disposition on the protein surface (Amodeo et al., 2014; Guardiani et al., 2016). There are only indications of experimentally determined disulfide bridges in VDAC (Reina et al., 2016). Interestingly, the insertion of negative charges by –SO3H formation can modify the protein conformation by electric repulsions inside the chain or toward phospholipids. Some authors suggested that these conformational changes can initiate protein incorporation into mitochondria-derived vesicles (MDVs), later targeted to lysosomes. MDVs, whose production is induced by mitochondrial stress (Sugiura et al., 2014), contain numerous oxidized proteins derived mainly from the OMM. VDAC1 appears among these proteins (Soubannier et al., 2012), whereas no evidence about the VDAC3 cysteine status or its presence within MDVs has been reported.

Posttranslational removal of N-terminal methionine leaves VDAC3 Cys2 as the first amino acid. All proteomic analysis revealed this residue as completely acetylated, a PTM that, beside significantly influencing protein stability, activity, folding, and localization, triggers specific protein degradation by the Ac/N–degron pathway. This peculiar branch of the N-rule pathway controls several aspects of protein quality control (Eldeeb and Fahlman, 2016). Any dysregulation of N-terminal acetylation, which protects cysteine against irreversible oxidation, leads to serious pathological conditions including neurodegenerative diseases, cancers, hypertension, and X-linked genetic disorders (Rope et al., 2011; Aksnes et al., 2015, 2016).

Possible Pathological Implications of Redox Cysteine Modifications in VDACs

Despite VDAC being involved in several mitochondria-associated pathologies such as hypertension (Alvira et al., 2012; Zhu et al., 2019), diabetes (Sasaki et al., 2012; Zhang et al., 2019), cancer (Maldonado and Lemasters, 2012; Reina and De Pinto, 2017; Magrì et al., 2018), cardiovascular diseases (Karachitos et al., 2017), and neurodegenerative disorders, such as ALS (Magrì et al., 2016), PD (Rostovtseva et al., 2015), and AD (Manczak and Reddy, 2012), very little is known about the correlation of its Cys-Ox-PTMs with the aforementioned diseases. This is undoubtedly attributable to the few data available on these VDAC modifications, whose identification and analysis results are challenging. To quote some examples, Zahid et al. reported alterations in the S-nitrosylation pattern of VDAC2 in AD (Zhao et al., 2015); in Lam et al. (2010), the authors listed VDAC1, VDAC2, and VDAC3 among the S-nitrosylated targets in the prostate epithelial cell line.

The role of VDAC cysteine disulfide bridges or of its oxidation to sulfonic acid has yet to be discovered. It is suggested that the degree of VDAC oxidation on the mitochondrial surface could function as a signal of ROS load within the mitochondrial network, even if much more needs to be done to truly understand their biological meaning.

Statements

Author contributions

SR conceived the review organization and wrote most of it. MP selected the references and contributed to the oxidized cysteine discovery. FG, AM, and VD wrote the discussion section. SF and RS reviewed and revised the text and the Mass Spectroscopy analysis.

Funding

This research was partly supported by a grant from “Piano della Ricerca di Ateneo 2016–2018” of the University of Catania, Italy, to all the participants and by a MIUR PNR “Proof of Concept 2018” grant (codex: PEPSLA) to AM.

Acknowledgments

The authors gratefully acknowledge the Bio-Nanotech Research and Innovation Tower (BRIT; PON project financed by the Italian Ministry for Education, University and Research MIUR) and Dr. Salvatore A.M. Cubisino (University of Catania) for modeling human VDAC structures.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    Abu-HamadS.ZaidH.IsraelsonA.NahonE.Shoshan-BarmatzV. (2008). Hexokinase-I protection against apoptotic cell death is mediated via interaction with the voltage-dependent anion channel-1: mapping the site of binding.J. Biol. Chem.2831348213490. 10.1074/jbc.M708216200

  • 2

    AksnesH.DrazicA.ArnesenT. (2015). (Hyper)tension release by N-terminal acetylation.Trends Biochem. Sci.40422424. 10.1016/j.tibs.2015.05.003

  • 3

    AksnesH.DrazicA.MarieM.ArnesenT. (2016). First things first: vital protein marks by N-terminal acetyltransferases.Trends Biochem. Sci.41746760. 10.1016/j.tibs.2016.07.005

  • 4

    AkterS.FuL.JungY.Lo ConteM.LawsonJ. R.LowtherW. T.et al (2018). Chemical proteomics reveals new targets of cysteine sulfinic acid reductase.Nat. Chem. Biol.149951004. 10.1038/s41589-018-0116-2

  • 5

    AlcockL. J.PerkinsM. V.ChalkerJ. M. (2018). Chemical methods for mapping cysteine oxidation.Chem. Soc. Rev.47231268. 10.1039/c7cs00607a

  • 6

    AllenE. M.MieyalJ. J. (2012). Protein-thiol oxidation and cell death: regulatory role of glutaredoxins.Antioxid. Redox Signal.1717481763. 10.1089/ars.2012.4644

  • 7

    AllenJ. W.FergusonS. J.GingerM. L. (2008). Distinctive biochemistry in the trypanosome mitochondrial intermembrane space suggests a model for stepwise evolution of the MIA pathway for import of cysteine-rich proteins.FEBS Lett.58228172825. 10.1016/j.febslet.2008.07.015

  • 8

    AlviraC. M.UmeshA.HustedC.YingL.HouY.LyuS. C.et al (2012). Voltage-dependent anion channel-2 interaction with nitric oxide synthase enhances pulmonary artery endothelial cell nitric oxide production.Am. J. Respir. Cell Mol.47669678. 10.1165/rcmb.2011-0436OC

  • 9

    AmodeoG. F.ScorciapinoM. A.MessinaA.De PintoV.CeccarelliM. (2014). Charged residues distribution modulates selectivity of the open state of human isoforms of the voltage dependent anion-selective channel.PLoS One9:e103879. 10.1371/journal.pone.0103879

  • 10

    AnandP.StamlerJ. S. (2012). Enzymatic mechanisms regulating protein S-nitrosylation: implications in health and disease.J. Mol. Med.90233244. 10.1007/s00109-012-0878-z

  • 11

    AnjoS. I.MeloM. N.LoureiroL. R.SabalaL.CastanheiraP.GrãosM.et al (2019). oxSWATH: An integrative method for a comprehensive redox-centered analysis combined with a generic differential proteomics screening.Redox Biol.22:101130. 10.1016/j.redox.2019.101130

  • 12

    AramL.GeulaS.ArbelN.Shoshan-BarmatzV. (2010). VDAC1 cysteine residues: topology and function in channel activity and apoptosis.Biochem. J.427445454. 10.1042/BJ20091690

  • 13

    ArbelN.Ben-HailD.Shoshan-BarmatzV. (2012). Mediation of the antiapoptotic activity of Bcl-xL protein upon interaction with VDAC1 protein.J. Biol. Chem.2872315223161. 10.1074/jbc.M112.345918

  • 14

    BakD. W.PizzagalliM. D.WeerapanaE. (2017). Identifying functional cysteine residues in the mitochondria.ACS Chem. Biol.12947957. 10.1021/acschembio.6b01074

  • 15

    BanciL.BertiniI.CefaroC.Ciofi-BaffoniS.GalloA.MartinelliM.et al (2009). MIA40 is an oxidoreductase that catalyzes oxidative protein folding in mitochondria.Nat. Struct. Mol. Biol.16198206. 10.1038/nsmb.1553

  • 16

    BanciL.BertiniI.Ciofi-BaffoniS.BoscaroF.ChatziA.MikolajczykM.et al (2011). Anamorsin is a [2Fe-2S] cluster-containing substrate of the Mia40-dependent mitochondrial protein trapping machinery.Chem. Biol.18794804. 10.1016/j.chembiol.2011.03.015

  • 17

    BarrosM. H.McStayG. P. (2020). Modular biogenesis of mitochondrial respiratory complexes.Mitochondrion5094114. 10.1016/j.mito.2019.10.008

  • 18

    BayrhuberM.MeinsT.HabeckM.BeckerS.GillerK.VillingerS.et al (2008). Structure of the human voltage-dependent anion channel.Proc. Natl. Acad. Sci. U.S.A.1051537015375. 10.1073/pnas.0808115105

  • 19

    BechtelT. J.WeerapanaE. (2017). From structure to redox: the diverse functional roles of disulfides and implications in disease.Proteomics17:10.1002/pmic.201600391. 10.1002/pmic.201600391

  • 20

    BenzR. (1994). Permeation of hydrophilic solutes through mitochondrial outer membranes: review on mitochondrial porins.Biochim. Biophys. Acta1197167196. 10.1016/0304-4157(94)90004-3

  • 21

    BienM.LongenS.WagenerN.ChwallaI.HerrmannJ. M.RiemerJ. (2010). Mitochondrial disulfide bond formation is driven by intersubunit electron transfer in Erv1 and proofread by glutathione.Mol. Cell37516528. 10.1016/j.molcel.2010.01.017

  • 22

    BiteauB.LabarreJ.ToledanoM. B. (2003). ATP-dependent reduction of cysteine-sulphinic acid by S. cerevisiae sulphiredoxin.Nature425980984. 10.1038/nature02075

  • 23

    Blachly-DysonE.ZambroniczE. B.YuW. H.AdamsV.McCabeE. R.AdelmanJ.et al (1993). Cloning and functional expression in yeast of two human isoforms of the outer mitochondrial membrane channel, the voltage-dependent anion channel.J. Biol. Chem.26818351841.

  • 24

    BlackintonJ.LakshminarasimhanM.ThomasK. J.AhmadR.GreggioE.RazaA. S.et al (2009). Formation of a stabilized cysteine sulfinic acid is critical for the mitochondrial function of the parkinsonism protein DJ-1.J. Biol. Chem.28464766485. 10.1074/jbc.M806599200

  • 25

    BleierL.WittigI.HeideH.StegerM.BrandtU.DröseS. (2015). Generator-specific targets of mitochondrial reactive oxygen species.Free Radic. Biol. Med.78110. 10.1016/j.freeradbiomed.2014.10.511

  • 26

    BulteauA. L.SzwedaL. I.FriguetB. (2006). Mitochondrial protein oxidation and degradation in response to oxidative stress and aging.Exp. Gerontol.41653657. 10.1016/j.exger.2006.03.013

  • 27

    CalvoS. E.ClauserK. R.MoothaV. K. (2016). MitoCarta2.0: an updated inventory of mammalian mitochondrial proteins.Nucleic Acids Res.4412511257. 10.1093/nar/gkv1003

  • 28

    CamaraA. K.BienengraeberM.StoweD. F. (2011). Mitochondrial approaches to protect against cardiac ischemia and reperfusion injury.Front. Physiol.2:13. 10.3389/fphys.2011.00013

  • 29

    CamargoL. L.HarveyA. P.RiosF. J.TsiropoulouS.Da SilvaR. N. O.CaoZ.et al (2018). Vascular Nox (NADPH Oxidase) compartmentalization, protein hyperoxidation, and endoplasmic reticulum stress response in hypertension.Hypertension72235246. 10.1161/HYPERTENSIONAHA.118.10824

  • 30

    Canet-AvilésR. M.WilsonM. A.MillerD. W.AhmadR.McLendonC.BandyopadhyayS.et al (2004). The Parkinson’s disease protein DJ-1 is neuroprotective due to cysteine-sulfinic acid-driven mitochondrial localization.Proc. Natl. Acad. Sci. U.S.A.10191039108. 10.1073/pnas.0402959101

  • 31

    CavallaroG. (2010). Genome-wide analysis of eukaryotic twin CX9C proteins.Mol. Biosyst.624592470. 10.1039/c0mb00058b

  • 32

    ChakrabortyJ.BassoV.ZivianiE. (2017). Post translational modification of Parkin.Biol. Direct12:6. 10.1186/s13062-017-0176-3

  • 33

    ChamberlainL. H.ShipstonM. J. (2015). The physiology of protein S-acylation.Physiol. Rev.95341376. 10.1152/physrev.00032.2014

  • 34

    ChangT. S.JeongW.WooH. A.LeeS. M.ParkS.RheeS. G. (2004). Characterization of mammalian sulfiredoxin and its reactivation of hyperoxidized peroxiredoxin through reduction of cysteine sulfinic acid in the active site to cysteine.J. Biol. Chem.2795099451001. 10.1074/jbc.m409482200

  • 35

    ChatziA.TokatlidisK. (2013). The mitochondrial intermembrane space: a hub for oxidative folding linked to protein biogenesis.Antioxid. Redox Signal.195462. 10.1089/ars.2012.4855

  • 36

    ChecchettoV.ReinaS.MagrìA.SzaboI.De PintoV. (2014). Recombinant human voltage dependent anion selective channel isoform 3 (hVDAC3) forms pores with a very small conductance.Cell. Physiol. Biochem.34842853. 10.1159/000363047

  • 37

    ChengE. H.SheikoT. V.FisherJ. K.CraigenW. J.KorsmeyerS. J. (2003). VDAC2 inhibits BAK activation and mitochondrial apoptosis.Science301513517. 10.1126/science.1083995

  • 38

    ChinH. S.LiM. X.TanI. K. L.NinniS. R. L.DewsonG. (2018). VDAC2 enables BAX to mediate apoptosis and limit tumor development.Nat. Commun.9:4976. 10.1038/s41467-018-07309-4

  • 39

    ChouchaniE. T.KazakL.JedrychowskiM. P.LuG. Z.EricksonB. K.SzpytJ.et al (2016). Mitochondrial ROS regulate thermogenic energy expenditure and sulfenylation of UCP1.Nature532112116. 10.1038/nature17399

  • 40

    ChouchaniE. T.MethnerC.NadtochiyS. M.LoganA.PellV. R.DingS.et al (2013). Cardioprotection by S-nitrosation of a cysteine switch on mitochondrial complex I.Nat. Med.19753759. 10.1038/nm.3212

  • 41

    ChungH. S.WangS. B.VenkatramanV.MurrayC. I.Van EykJ. E. (2013). Cysteine oxidative post translational modifications: emerging regulation in the cardiovascular system.Circ. Res.112382392. 10.1161/CIRCRESAHA.112.268680

  • 42

    ClampM.FryB.KamalM.XieX.CuffJ.LinM. F.et al (2007). Distinguishing protein-coding and noncoding genes in the human genome.Proc. Natl. Acad. Sci. U.S.A.1041942819433. 10.1073/pnas.0709013104

  • 43

    ColombiniM.Blachly-DysonE.ForteM. (1996). VDAC, a channel in the outer mitochondrial membrane.Ion Channels4169202. 10.1007/978-1-4899-1775-1_5

  • 44

    ColussiC.AlbertiniM. C.CoppolaS.RovidatiS.GalliF.GhibelliL. (2000). H2O2-induced block of glycolysis as an active ADP-ribosylation reaction protecting cells from apoptosis.FASEB J.1422662276. 10.1096/fj.00-0074com

  • 45

    DabirD. V.LeverichE. P.KimS. K.TsaiF. D.HirasawaM.KnaffD. B.et al (2007). A role for cytochrome c and cytochrome c peroxidase in electron shuttling from Erv1.EMBO J.2648014811. 10.1038/sj.emboj.7601909

  • 46

    Dalle-DonneI.RossiR.ColomboG.GiustariniD.MilzaniA. (2009). Protein S-Glutathionylation: a regulatory device from bacteria to humans.Trends Biochem. Sci.348596. 10.1016/j.tibs.2008.11.002

  • 47

    De PintoV.Al JamalJ. A.BenzR.GenchiG.PalmieriF. (1991). Characterization of SH groups in porin of bovine heart mitochondria. Porin cysteines are localized in the channel walls.Eur. J. Biochem.202903911. 10.1111/j.1432-1033.1991.tb16450.x

  • 48

    De PintoV.GuarinoF.GuarneraA.MessinaA.ReinaS.TomaselloF. M.et al (2010). Characterization of human VDAC isoforms: a peculiar function for VDAC3?Biochim. Biophys. Acta179712681275. 10.1016/j.bbabio.2010.01.031

  • 49

    De PintoV.ReinaS.GuptaA.MessinaA.MahalakshmiR. (2016). Role of cysteines in mammalian VDAC isoforms’ function.Biochim. Biophys. Acta185712191227. 10.1016/j.bbabio.2016.02.020

  • 50

    De PintoV.TomaselloF.MessinaA.GuarinoF.BenzR.La MendolaD.et al (2007). Determination of the conformation of the human VDAC1 N-terminal peptide, a protein moiety essential for the functional properties of the pore.Chem. Biochem.8744756. 10.1002/cbic.200700009

  • 51

    DelaunayA.PfliegerD.BarraultM. B.VinhJ.ToledanoM. B. (2002). A thiol peroxidase is an H2O2 receptor and redoxtransducer in gene activation.Cell111471481. 10.1016/s0092-8674(02)01048-6

  • 52

    DepuydtM.MessensJ.ColletJ. F. (2011). How proteins form disulfide bonds.Antioxid. Redox Signal.154966. 10.1089/ars.2010.3575

  • 53

    DuanG.WaltherD. (2015). The roles of post-translational modifications in the context of protein interaction networks.PLoS Comput. Biol.11:e1004049. 10.1371/journal.pcbi.1004049

  • 54

    DuranskiM. R.GreerJ. J.DejamA.JaganmohanS.HoggN.LangstonW.et al (2005). Cytoprotective effects of nitrite during in vivo ischemia-reperfusion of the heart and liver.J. Clin. Invest.11512321240. 10.1172/jci22493

  • 55

    EdgarR. S.GreenE. W.ZhaoY.van OoijenG.OlmedoM.QinX.et al (2012). Peroxiredoxins are conserved markers of circadian rhythms.Nature485459464. 10.1038/nature11088

  • 56

    EldeebM.FahlmanR. (2016). The-N-end rule: the beginning determines the end.Protein Pept. Lett.23343348. 10.2174/0929866523666160108115809

  • 57

    FeldmanJ. L.Dittenhafer-ReedK. E.KudoN.ThelenJ. N.ItoA.YoshidaM.et al (2015). Kinetic and structural basis for Acyl-group selectivity and NAD (+) dependence in sirtuin-catalyzed deacylation.Biochemistry5430373050. 10.1021/acs.biochem.5b00150

  • 58

    FernandoV.ZhengX.WaliaY.SharmaV.LetsonJ.FurutaS. (2019). S-Nitrosylation: an emerging paradigm of redox signaling.Antioxidants8:E404. 10.3390/antiox8090404

  • 59

    Ferrer-SuetaG.MantaB.BottiH.RadiR.TrujilloM.DenicolaA. (2011). Factors affecting protein thiol reactivity and specificity in peroxide reduction.Chem. Res. Toxicol.24434450. 10.1021/tx100413v

  • 60

    FiekC.BenzR.RoosN.BrdiczkaD. (1982). Evidence for identity between the hexokinase-binding protein and the mitochondrial porin in the outer membrane of rat liver mitochondria.Biochim. Biophys. Acta688429440. 10.1016/0005-2736(82)90354-6

  • 61

    ForresterM. T.StamlerJ. S. (2007). A classification scheme for redox-based modifications of proteins.Am. J. Respir. Cell Mol. Biol.36135137. 10.1165/rcmb.2006-001ed

  • 62

    FridovichI. (1995). Superoxide radical and superoxide dismutases.Annu. Rev. Biochem.6497112. 10.1146/annurev.bi.64.070195.000525

  • 63

    FujiwaraN.NakanoM.KatoS.YoshiharaD.OokawaraT.EguchiH.et al (2007). Oxidative modification to cysteine sulfonic acid of Cys111 in human copper-zinc superoxide dismutase.J. Biol. Chem.2823593335944. 10.1074/jbc.m702941200

  • 64

    FurukawaY.FuR.DengH. X.SiddiqueT.O’HalloranT. V. (2006). Disulfide cross-linked protein represents a significant fraction of ALS-associated Cu, Zn-superoxide dismutase aggregates in spinal cords of model mice.Proc. Natl. Acad. Sci. U.S.A.10371487153. 10.1073/pnas.0602048103

  • 65

    GaillotC.DelolmeF.FabreL.CharreyreM. T.LadavièreC.FavierA. (2020). Taking advantage of oxidation to characterize thiol-containing polymer chains by MALDI-TOF mass spectrometry.Anal. Chem.9238043809. 10.1021/acs.analchem.9b05207

  • 66

    GrekC. L.ZhangJ.ManevichY.TownsendD. M.TewK. D. (2013). Causes and consequences of cysteine S-Glutathionylation.J. Biol. Chem.2882649726504. 10.1074/jbc.R113.461368

  • 67

    GriesserM.ChauvinJ. R.PrattD. A. (2018). The hydrogen atom transfer reactivity of sulfinic acids.Chem. Sci.972187229. 10.1039/c8sc02400f

  • 68

    GrossD. P.BurgardC. A.ReddehaseS.LeitchJ. M.CulottaV. C.HellK. (2011). Mitochondrial Ccs1 contains a structural disulfide bond crucial for the import of this unconventional substrate by the disulfide relay system.Mol. Biol. Cell2237583767. 10.1091/mbc.E11-04-0296

  • 69

    GuardianiC.LeggioL.ScorciapinoM. A.De PintoV.CeccarelliM. (2016). A computational study of ion current modulation in hVDAC3 induced by disulfide bonds.Biochim. Biophys. Acta1858:813823. 10.1016/j.bbamem.2016.01.013

  • 70

    HallA.ParsonageD.PooleL. B.KarplusP. A. (2010). Structural evidence that peroxiredoxin catalytic power is based on transition-state stabilization.J. Mol. Biol.402194209. 10.1016/j.jmb.2010.07.022

  • 71

    HanD.AntunesF.CanaliR.RettoriD.CadenasE. (2003). Voltage-dependent anion channels control the release of the superoxide anion from mitochondria to cytosol.J. Biol. Chem.27855575563. 10.1074/jbc.m210269200

  • 72

    HeppnerD. E.DustinC. M.LiaoC.HristovaM.VeithC.LittleA. C.et al (2018). Direct cysteine sulfenylation drives activation of the Src kinase.Nat. Commun.9:4522. 10.1038/s41467-018-06790-1

  • 73

    HerrmannJ. M.RiemerJ. (2012). Mitochondrial disulfide relay: redox-regulated protein import into the intermembrane space.J. Biol. Chem.28744264433. 10.1074/jbc.R111.270678

  • 74

    HillerS.GarcesR. G.MaliaT. J.OrekhovV. Y.ColombiniM.WagnerG. (2008). Solution structure of the integral human membrane protein VDAC-1 in detergent micelles.Science32112061210. 10.1126/science.1161302

  • 75

    HolmgrenA. (2000). Antioxidant function of thioredoxin and glutaredoxin systems.Antioxid. Redox Signal.2811820. 10.1089/ars.2000.2.4-811

  • 76

    HolmgrenA.JohanssonC.BerndtC.LönnM. E.HudemannC.LilligC. H. (2005). Thiol redox control via thioredoxin and glutaredoxin systems.Biochem. Soc. Trans.3313751377. 10.1042/bst0331375

  • 77

    HoloubekA.VečeřJ.SiglerK. (2007). Monitoring of the proton electrochemical gradient in reconstituted vesicles: quantitative measurements of both transmembrane potential and intravesicular pH by ratiometric fluorescent probes.J. Fluoresc.17201213. 10.1007/s10895-007-0159-3

  • 78

    HuJ.DongL.OuttenC. E. (2008). The redox environment in the mitochondrial intermembrane space is maintained separately from the cytosol and matrix.J Biol. Chem.2832912629134. 10.1074/jbc.M803028200

  • 79

    HuangH.HuX.EnoC. O.ZhaoG.LiC.WhiteC. (2013). An interaction between Bcl-xL and the voltage-dependent anion channel (VDAC)promotes mitochondrial Ca2+ uptake.J. Biol. Chem.2881987019881. 10.1074/jbc.M112.448290

  • 80

    IciekM.Kowalczyk-PachelD.Bilska-WilkoszA.KwiecieńI.GórnyM.WłodekL. (2015). S-sulfhydration as a cellular redox regulation.Biosci. Rep.36:e00304. 10.1042/BSR20150147

  • 81

    International Human Genome Sequencing Consortium. (2004). Finishing the euchromatic sequence of the human genome.Nature431931945. 10.1038/nature03001

  • 82

    JangH. H.LeeK. O.ChiY. H.JungB. G.ParkS. K.ParkJ. H.et al (2004). Two enzymes in one; two yeast peroxiredoxins display oxidative stress-dependent switching from a peroxidase to a molecular chaperone function.Cell117625635.

  • 83

    JarvisR. M.HughesS. M.LedgerwoodE. C. (2012). Peroxiredoxin 1 functions as a signal peroxidase to receive, transduce, and transmit peroxide signals in mammalian cells.Free Radic. Biol. Med.5315221530. 10.1016/j.freeradbiomed.2012.08.001

  • 84

    JonesD. P. (2006). Redefining oxidative stress.Antioxid. Redox Signal.818651879. 10.1089/ars.2006.8.1865

  • 85

    KarachitosA.JordanJ.KmitaH. (2017). VDAC-targeted drugs affecting cytoprotection and mitochondrial physiology in cerebrovascular and cardiovascular diseases.Curr. Med. Chem.2444194434. 10.2174/0929867324666170530073238

  • 86

    KimD.LimS.HaqueM.RyooN.HongH. S.RhimH.et al (2015). Identification of disulfide cross-linked tau dimer responsible for tau propagation.Sci. Rep.5:15231. 10.1038/srep15231

  • 87

    KojerK.BienM.GangelH.MorganB.DickT. P.RiemerJ. (2012). Glutathione redox potential in the mitochondrial intermembrane space is linked to the cytosol and impacts the Mia40 redox state.EMBO J.3131693182. 10.1038/emboj.2012.165

  • 88

    KojerK.PelehV.CalabreseG.HerrmannJ. M.RiemerJ. (2015). Kinetic control by limiting glutaredoxin amounts enables thiol oxidation in the reducing mitochondrial intermembrane space.Mol. Biol. Cell26195204. 10.1091/mbc.E14-10-1422

  • 89

    KomarovA. G.DengD.CraigenW. J.ColombiniM. (2005). New insights into the mechanism of permeation through large channels.Biophys. J.8939503959. 10.1529/biophysj.105.070037

  • 90

    KordyukovaL. V.SerebryakovaM. V.BaratovaL. A.VeitM. (2010). Site-specific attachment of palmitate or stearate to cytoplasmic versus transmembrane cysteines is a common feature of viral spike proteins.Virology3984956. 10.1016/j.virol.2009.11.039

  • 91

    KrishnanN.FuC.PappinD. J.TonksN. K. (2011). H2S-Induced sulfhydration of the phosphatase PTP1B and its role in the endoplasmic reticulum stress response.Sci. Signal.4:ra86. 10.1126/scisignal.2002329

  • 92

    LamY. W.YuanY.IsaacJ.BabuC. V.MellerJ.HoS. M. (2010). Comprehensive identification and modified-site mapping of S-nitrosylated targets in prostate epithelial cells.PLoS One5:e9075. 10.1371/journal.pone.0009075

  • 93

    LeanzaL.ZorattiM.GulbinsE.SzabòI. (2014). Mitochondrial ion channels as oncological targets.Oncogene3355695581. 10.1038/onc.2013.578

  • 94

    LeeJ.HofhausG.LisowskyT. (2000). Erv1p from Saccharomyces cerevisiae is a FAD-linked sulfhydryl oxidase.FEBS Lett.4776266. 10.1016/s0014-5793(00)01767-1

  • 95

    LimJ. C.ChoiH. I.ParkY. S.NamH. W.WooH. A.KwonK. S.et al (2008). Irreversible oxidation of the active-site cysteine of peroxiredoxin to cysteine sulfonic acid for enhanced molecular chaperone activity.J. Biol. Chem.2832887328880. 10.1074/jbc.M804087200

  • 96

    LiuZ.LuoQ.GuoC. (2015). Bim and VDAC1 are hierarchically essential for mitochondrial ATF2 mediated cell death.Cancer Cell Int.15:34. 10.1186/s12935-015-0188-y

  • 97

    LlopisJ.McCafferyJ. M.MiyawakiA.FarquharM. G.TsienR. Y. (1998). Measurement of cytosolic, mitochondrial, and Golgi pH in single living cells with green fluorescent proteins.Proc. Natl. Acad. Sci. U.S.A.9568036808. 10.1073/pnas.95.12.6803

  • 98

    Lo ConteM.CarrollK. S. (2013). The redox biochemistry of protein sulfenylation and sulfinylation.J. Biol. Chem.2882648026488. 10.1074/jbc.R113.467738

  • 99

    LongenS.WoellhafM. W.PetrungaroC.RiemerJ.HerrmannJ. M. (2014). The disulfide relay of the intermembrane space oxidizes the ribosomal subunit mrp10 on its transit into the mitochondrial matrix.Dev. Cell283042. 10.1016/j.devcel.2013.11.007

  • 100

    LushchakV. I. (2012). Glutathione homeostasis and functions: potential targets for medical interventions.J. Amino Acids2012:736837. 10.1155/2012/736837

  • 101

    MagrìA.BelfioreR.ReinaS.TomaselloM. F.Di RosaM. C.GuarinoF.et al (2016). Hexokinase I N-terminal based peptide prevents the VDAC1-SOD1 G93A interaction and re-establishes ALS cell viability.Sci Rep.6:34802. 10.1038/srep34802

  • 102

    MagrìA.ReinaS.De PintoV. (2018). VDAC1 as pharmacological target in cancer and neurodegeneration: focus on its role in apoptosis.Front. Chem.6:108. 10.3389/fchem.2018.00108

  • 103

    MaillouxR. J.JinX.WillmoreW. G. (2014). Redox regulation of mitochondrial function with emphasis on cysteine oxidation reactions.Redox Biol.2123139. 10.1016/j.redox.2013.12.011

  • 104

    MajumderS.FiskH. A. (2013). VDAC3 and Mps1 negatively regulate ciliogenesis.Cell Cycle12849858. 10.4161/cc.23824

  • 105

    MaldonadoE. N.LemastersJ. J. (2012). Warburg revisited: regulation of mitochondrial metabolism by voltage-dependent anion channels in cancer cells.J. Pharmacol. Exp. Ther.342637641. 10.1124/jpet.112.192153

  • 106

    MaldonadoE. N.PatnaikJ.MullinsM. R.LemastersJ. J. (2010). Free tubulin modulates mitochondrial membrane potential in cancer cells.Cancer Res.701019210201. 10.1158/0008-5472.CAN-10-2429

  • 107

    ManczakM.ReddyP. H. (2012). Abnormal interaction of VDAC1 with amyloid beta and phosphorylated tau causes mitochondrial dysfunction in Alzheimer’s disease.Hum. Mol. Genet.2151315146. 10.1093/hmg/dds360

  • 108

    MauryaS. R.MahalakshmiR. (2013). Modulation of human mitochondrial voltage-dependent anion channel 2 (hVDAC-2) structural stability by cysteine-assisted barrel-lipid interactions.J. Biol. Chem.2882558425592. 10.1074/jbc.M113.493692

  • 109

    MauryaS. R.MahalakshmiR. (2014). Influence of protein – micelle ratios and cysteine residues on the kinetic stability and unfolding rates of human mitochondrial VDAC-2.PLoS One9:e87701. 10.1371/journal.pone.0087701

  • 110

    MauryaS. R.MahalakshmiR. (2015). N-helix and cysteines inter-regulate human mitochondrial VDAC-2 function and biochemistry.J. Biol. Chem.2903024030252. 10.1074/jbc.M115.693978

  • 111

    MauryaS. R.MahalakshmiR. (2016). VDAC-2: Mitochondrial outer membrane regulator masquerading as a channel?FEBS J.28318311836. 10.1111/febs.13637

  • 112

    MenzelV. A.CassaráM. C.BenzR.De PintoV.MessinaA.CunsoloV.et al (2009). Molecular and functional characterization of VDAC2 purified from mammal spermatozoa.Biosci. Rep.29351362. 10.1042/BSR20080123

  • 113

    MessinaA.ReinaS.GuarinoF.De PintoV. (2012). VDAC isoforms in mammals.Biochim. Biophys. Acta181814661476. 10.1016/j.bbamem.2011.10.005

  • 114

    MuronetzV. I.MelnikovaA. K.SasoL.SchmalhausenE. V. (2020). Influence of oxidative stress on catalytic and non-glycolytic functions of glyceraldehyde-3-phosphate dehydrogenase.Curr. Med. Chem.2720402058. 10.2174/0929867325666180530101057

  • 115

    MurrayC. I.Van EykJ. E. (2012). Chasing cysteine oxidative modifications: proteomic tools for characterizing cysteine redox status.Circ. Cardiovasc. Genet.5:591. 10.1161/circgenetics.111.961425

  • 116

    NaghdiS.HajnóczkyG. (2016). VDAC2-specific cellular functions and the underlying structure.Biochim. Biophys. Acta186325032514. 10.1016/j.bbamcr.2016.04.020

  • 117

    NaghdiS.VárnaiP.HajnóczkyG. (2015). Motifs of VDAC2 required for mitochondrial Bak import and tBid-induced apoptosis.Proc. Natl. Acad. Sci. U.S.A.112E5590E5599. 10.1073/pnas.1510574112

  • 118

    NakamuraT.TuS.AkhtarM. W.SunicoC. R.OkamotoS.LiptonS. A. (2013). Aberrant protein s-nitrosylation in neurodegenerative diseases.Neuron78596614. 10.1016/j.neuron.2013.05.005

  • 119

    NeumannC. A.KrauseD. S.CarmanC. V.DasS.DubeyD. P.AbrahamJ. L.et al (2003). Essential role for the peroxiredoxin Prdx1 in erythrocyte antioxidant defence and tumour suppression.Nature424561565. 10.1038/nature01819

  • 120

    NiemannM.WieseS.ManiJ.ChanfonA.JacksonC.MeisingerC.et al (2013). Mitochondrial outer membrane proteome of Trypanosoma brucei reveals novel factors required to maintain mitochondrial morphology.Mol. Cell. Proteomics12515528. 10.1074/mcp.M112.023093

  • 121

    OkazakiM.KurabayashiK.AsanumaM.SaitoY.DodoK.SodeokaM. (2015). VDAC3 gating is activated by suppression of disulfide-bond formation between the N-terminal region and the bottom of the pore.Biochim. Biophys. Acta184831883196. 10.1016/j.bbamem.2015.09.017

  • 122

    O’NeillJ. S.ReddyA. B. (2011). Circadian clocks in human red blood cells.Nature469498503. 10.1038/nature09702

  • 123

    ParkJ. H.ZhuangJ.LiJ.HwangP. M. (2016). p53 as guardian of the mitochondrial genome.FEBS Lett.590924934. 10.1002/1873-3468.12061

  • 124

    PaulB. D.SnyderS. H. (2018). Gasotransmitter hydrogen sulfide signaling in neuronal health and disease.Biochem. Pharmacol.149101109. 10.1016/j.bcp.2017.11.019

  • 125

    PedroM. P.VilcaesA. A.GomezG. A.DaniottiJ. L. (2017). Individual S-acylated cysteines differentially contribute to H-Ras endomembrane trafficking and acylation/deacylation cycles.Mol. Biol. Cell28962974. 10.1091/mbc.E16-08-0603

  • 126

    PetrungaroC.ZimmermannK. M.KüttnerV.FischerM.DengjelJ.BogeskiI.et al (2015). The Ca2+-dependent release of the Mia40-Induced MICU1-MICU2 dimer from MCU regulates mitochondrial Ca2+ uptake.Cell Metab.22721733. 10.1016/j.cmet.2015.08.019

  • 127

    PiantadosiC. (2011). Regulation of mitochondrial processes by protein S-nitrosylation.Biochim. Biophys. Acta1820712721. 10.1016/j.bbagen.2011.03.008

  • 128

    PiroliG. G.ManuelA. M.ClapperA. C.WallaM. D.BaatzJ. E.PalmiterR. D.et al (2016). Succination is Increased on Select Proteins in the Brainstem of the NADH dehydrogenase (ubiquinone) Fe-S protein 4 (Ndufs4) Knockout Mouse, a Model of Leigh Syndrome.Mol. Cell. Proteomics15445461. 10.1074/mcp.M115.051516

  • 129

    PittalàM. G. G.SalettiR.ReinaS.CunsoloV.De PintoV.FotiS. (2020). A high resolution mass spectrometry study reveals the potential of disulfide formation in human mitochondrial voltage-dependent anion selective channel isoforms (hVDACs).Int. J. Mol. Sci.21:E1468. 10.3390/ijms21041468

  • 130

    PooleL. B.KarplusP. A.ClaiborneA. (2004). Protein sulfenic acids in redox signaling.Annu. Rev. Pharmacol. Toxicol.44325347.

  • 131

    PooleL. B.NelsonK. J. (2008). Discovering mechanisms of signaling-mediated cysteine oxidation.Curr. Opin. Chem. Biol.121824. 10.1016/j.cbpa.2008.01.021

  • 132

    PortilloS. D.Ferrer-SuetaG.CoitiñoL. (2017). Insights into the mechanism of peroxiredoxin 6 sulfenic acid reduction by ascorbate.Free Radic. Biol. Med.112:31. 10.1016/j.freeradbiomed.2017.10.036

  • 133

    Queralt-MartínM.BergdollL.TeijidoO.MunshiN.JacobsD.KuszakA. J.et al (2020). A lower affinity to cytosolic proteins reveals VDAC3 isoform-specific role in mitochondrial biology.J. Gen. Physiol.152:e201912501. 10.1085/jgp.201912501

  • 134

    RadiR. (2013). Peroxynitrite, a stealthy biological oxidant.J. Biol. Chem.2882646426472. 10.1074/jbc.R113.472936

  • 135

    RayP. D.HuangB. W.TsujiY. (2012). Reactive oxygen species (ROS) homeostasis and redox regulation in cellular signaling.Cell Signal.24981990. 10.1016/j.cellsig.2012.01.008

  • 136

    ReinaS.ChecchettoV.SalettiR.GuptaA.ChaturvediD.GuardianiC.et al (2016). VDAC3 as a sensor of oxidative state of the intermembrane space of mitochondria: the putative role of cysteine residue modifications.Oncotarget722492268. 10.18632/oncotarget.6850

  • 137

    ReinaS.De PintoV. (2017). Anti-cancer compounds targeted to VDAC: potential and perspectives.Curr. Med. Chem.2444474469. 10.2174/0929867324666170530074039

  • 138

    ReinaS.PalermoV.GuarneraA.GuarinoF.MessinaA.MazzoniC.et al (2010). Swapping of the N-terminus of VDAC1 with VDAC3 restores full activity of the channel and confers anti-aging features to the cell.FEBS Lett.58428372844. 10.1016/j.febslet.2010.04.066

  • 139

    ReinhardtC.ArenaG.NedaraK.EdwardsR.BrennerC.TokatlidisK.et al (2020). AIF meets the CHCHD4/Mia40-dependent mitochondrial import pathway.Biochim. Biophys. Acta1866:165746. 10.1016/j.bbadis.2020.165746

  • 140

    RenD.KimH.TuH. C.WestergardT. D.FisherJ. K.RubensJ. A.et al (2009). The VDAC2-BAK rheostat controls thymocyte survival.Sci. Signal.2:ra48. 10.1126/scisignal.2000274

  • 141

    RoosG.MessensJ. (2011). Protein sulfenic acid formation: From cellular damage to redox regulation.Free Radic. Biol. Med.51314326. 10.1016/j.freeradbiomed.2011.04.031

  • 142

    RopeA. F.WangK.EvjenthR.XingJ.JohnstonJ. J.SwensenJ. J.et al (2011). Using VAAST to identify an X-linked disorder resulting in lethality in male infants due to N-terminal acetyltransferase deficiency.Am. J. Hum. Genet.892843. 10.1016/j.ajhg.2011.05.017

  • 143

    RoseP.MooreP. K.ZhuY. Z. (2017). H2S biosynthesis and catabolism: new insights from molecular studies.Cell. Mol. Life Sci.7413911412. 10.1007/s00018-016-2406-8

  • 144

    RostovtsevaT. K.GurnevP. A.ProtchenkoO.HoogerheideD. P.YapT. L.PhilpottC. C.et al (2015). α-synuclein shows high affinity interaction with voltage-dependent anion channel, suggesting mechanisms of mitochondrial regulation and toxicity in Parkinson disease.J. Biol. Chem.2901846718477. 10.1074/jbc.M115.641746

  • 145

    RothA. F.WanJ.BaileyA. O.SunB.KucharJ. A.GreenW. N.et al (2006). Global analysis of protein palmitoylation in yeast.Cell12510031013. 10.1016/j.cell.2006.03.042

  • 146

    SalettiR.ReinaS.PittalàM. G. G.BelfioreR.CunsoloV.MessinaA.et al (2017). High resolution mass spectrometry characterization of the oxidation pattern of methionine and cysteine residues in rat liver mitochondria voltage-dependent anion selective channel 3 (VDAC3).Biochim. Biophys. Acta Biomembr.1859301311. 10.1016/j.bbamem.2016.12.003

  • 147

    SalettiR.ReinaS.PittalàM. G. G.MagrìA.CunsoloV.FotiS.et al (2018). Post-translational modifications of VDAC1 and VDAC2 cysteines from rat liver mitochondria.Biochim. Biophys. Acta Bioenerg.1859806816. 10.1016/j.bbabio.2018.06.007

  • 148

    SampsonM. J.DeckerW. K.BeaudetA. L.RuitenbeekW.ArmstrongD.HicksM. J.et al (2001). Immotile sperm and infertility in mice lacking mitochondrial voltage-dependent anion channel type 3.J. Biol. Chem.2763920639212. 10.1074/jbc.m104724200

  • 149

    SampsonM. J.LovellR. S.CraigenW. J. (1997). The murine voltage-dependent anion channel gene family. Conserved structure and function.J. Biol. Chem.2721896618973. 10.1074/jbc.272.30.18966

  • 150

    Santo-DomingoJ.WiederkehrA.De MarchiU. (2015). Modulation of the matrix redox signaling by mitochondrial Ca2+.World J. Biol. Chem.6310323. 10.4331/wjbc.v6.i4.310

  • 151

    SasakiK.DonthamsettyR.HeldakM.ChoY. E.ScottB. T.MakinoA. (2012). VDAC: old protein with new roles in diabetes.Am. J. Physiol. Cell. Physiol.303C1055C1060. 10.1152/ajpcell.00087.2012

  • 152

    SchmittS.ProkischH.SchlunckT.CampD. G.AhtingU.WaizeneggerT.et al (2006). Proteome analysis of mitochondrial outer membrane from Neurospora crassa.Proteomics67280.

  • 153

    SchredelsekerJ.PazA.LópezC. J.AltenbachC.LeungC. S.DrexlerM. K.et al (2014). High resolution structure and double electron-electron resonance of the zebrafish voltage-dependent anion channel 2 reveal an oligomeric population.J. Biol. Chem.2891256612577. 10.1074/jbc.M113.497438

  • 154

    ShakirS.VinhJ.ChiappettaG. (2017). Quantitative analysis of the cysteine redoxome by iodoacetyl tandem mass tags.Anal. Bioanal. Chem.40938213830. 10.1007/s00216-017-0326-6

  • 155

    ShimizuH.SchredelsekerJ.HuangJ.LuK.NaghdiS.LuF.et al (2015). Mitochondrial Ca2+ uptake by the voltage-dependent anion channel 2 regulates cardiac rhythmicity.eLife4:e04801.

  • 156

    ShimizuS.MatsuokaY.ShinoharaY.YonedaY.TsujimotoY. (2001). Essential role of voltage-dependent anion channel in various forms of apoptosis in mammalian cells.J. Cell Biol.152237250. 10.1083/jcb.152.2.237

  • 157

    Shoshan-BarmatzV.De PintoV.ZweckstetterM.RavivZ.KeinanN.ArbelN. (2010). VDAC, a multi-functional mitochondrial protein regulating cell life and death.Mol. Aspects Med.31227285. 10.1016/j.mam.2010.03.002

  • 158

    ShuZ.JungM.BegerH. G.MarzinzigM.HanF.ButzerU.et al (1997). pH-dependent changes of nitric oxide, peroxynitrite, and reactive oxygen species in hepatocellular damage.Am. J. Physiol.27311181126. 10.1152/ajpgi.1997.273.5.G1118

  • 159

    SobottaM. C.LiouW.StockerS.TalwarD.OehlerM.RuppertT.et al (2015). Peroxiredoxin-2 and STAT3 form a redox relay for H2O2 signaling.Nat. Chem. Biol.116470. 10.1038/nchembio.1695

  • 160

    SoubannierV.RippsteinP.KaufmanB. A.ShoubridgeE. A.McBrideH. M. (2012). Reconstitution of mitochondria derived vesicle formation demonstrates selective enrichment of oxidized cargo.PLoS One7:e52830. 10.1371/journal.pone.0052830

  • 161

    SpadaroD.YunB. W.SpoelS. H.ChuC.WangY. Q.LoakeG. J. (2010). The redox switch: dynamic regulation of protein function by cysteine modifications.Physiol. Plant.138360371. 10.1111/j.1399-3054.2009.01307.x

  • 162

    SpickettC. M.PittA. R. (2012). Protein oxidation: role in signalling and detection by mass spectrometry.Amino Acids42521. 10.1007/s00726-010-0585-4

  • 163

    StipanukM. H.UekiI.DominyJ. E.Jr.SimmonsC. R.HirschbergerL. L. (2009). Cysteine dioxygenase: a robust system for regulation of cellular cysteine levels.Amino Acids375563. 10.1007/s00726-008-0202-y

  • 164

    StojanovskiD.MüllerJ. M.MilenkovicD.GuiardB.PfannerN.ChacinskaA. (2008). The MIA system for protein import into the mitochondrial intermembrane space.Biochim. Biophys. Acta1783610617. 10.1016/j.bbamcr.2007.10.004

  • 165

    SubediK. P.KimJ. C.KangM.SonM. J.KimY. S.WooS. H. (2011). Voltage-dependent anion channel 2 modulates resting Ca2+ sparks, but not action potential-induced Ca2+ signaling in cardiac myocytes.Cell Calcium49136143. 10.1016/j.ceca.2010.12.004

  • 166

    SugiuraA.McLellandG. L.FonE. A.McBrideH. M. (2014). A new pathway for mitochondrial quality control: mitochondrial-derived vesicles.EMBO J.3321422156. 10.15252/embj.201488104

  • 167

    SummaD.SpigaO.BerniniA.VendittiV.PrioraR.FrosaliS.et al (2007). Protein-thiol substitution or protein dethiolation by thiol/disulfide exchange reactions: the albumin model.Proteins69369378. 10.1002/prot.21532

  • 168

    TeijidoO.UjwalR.HillerdalC. O.KullmanL.RostovtsevaT. K.AbramsonJ. (2012). Affixing N-terminal α-helix to the wall of the voltage-dependent anion channel does not prevent its voltage gating.J. Biol. Chem.2871143711445. 10.1074/jbc.M111.314229

  • 169

    TerziyskaN.GrumbtB.BienM.NeupertW.HerrmannJ. M.HellK. (2007). The sulfhydryl oxidase Erv1 is a substrate of the Mia40-dependent protein translocation pathway.FEBS Lett.58110981102. 10.1016/j.febslet.2007.02.014

  • 170

    TiensonH. L.DabirD. V.NealS. E.LooR.HassonS. A.BoontheungP.et al (2009). Reconstitution of the mia40-erv1 oxidative folding pathway for the small tim proteins.Mol. Biol. Cell2034813490. 10.1091/mbc.E08-10-1062

  • 171

    TomaselloF. M.GuarinoF. M.ReinaS.MessinaA.De PintoV. (2013). The voltage-dependent anion selective channel 1 (VDAC1) topography in the mitochondrial outer membrane as detected in intact cell.PLoS One8:e81522. 10.1371/journal.pone.0081522

  • 172

    TruongT. H.UngP. M.PaldeP. B.PaulsenC. E.SchlessingerA.CarrollK. S. (2016). Molecular basis for redox activation of epidermal growth factor receptor kinase.Cell Chem. Biol.23837848. 10.1016/j.chembiol.2016.05.017

  • 173

    UjwalR.CascioD.ColletiercJ. P.FahamaS.ZhangaJ.TorodL.et al (2008). The crystal structure of mouse VDAC1 at 2.3 A resolution reveals mechanistic insights into metabolite gating.Proc. Natl. Acad. Sci. U.S.A.1051774217747. 10.1073/pnas.0809634105

  • 174

    UrmeyA. R.ZondloN. J. (2020). Structural preferences of cysteine sulfinic acid: the sulfinate engages in multiple local interactions with the peptide backbone.Free Radic. Biol. Med.14896107. 10.1016/j.freeradbiomed.2019.12.030

  • 175

    VahsenN.CandeC.BriereJ. J.BenitP.JozaN.LarochetteN.et al (2004). AIF deficiency compromises oxidative phosphorylation.EMBO J.2346794689. 10.1038/sj.emboj.7600461

  • 176

    Valdez TaubasJ.PelhamH. (2005). Swf1-dependent palmitoylation of the SNARE Tlg1 prevents its ubiquitination and degradation.EMBO J.2425242532. 10.1038/sj.emboj.7600724

  • 177

    Vander HeidenM. G.LiX. X.GottleibE.HillR. B.ThompsonC. B.ColombiniM. (2001). Bcl-xL promotes the open configuration of the voltage-dependent anion channel and metabolite passage through the outer mitochondrial membrane.J. Biol. Chem.2761941419419. 10.1074/jbc.m101590200

  • 178

    VögtleF. N.BurkhartJ. M.RaoS.GerbethC.HinrichsJ.MartinouJ. C.et al (2012). Intermembrane space proteome of yeast mitochondria.Mol. Cell. Proteomics1118401852. 10.1074/mcp.M112.021105

  • 179

    WalkerS.UllmanO.StultzC. M. (2012). Using intramolecular disulfide bonds in tau protein to deduce structural features of aggregation-resistant conformations.J. Biol. Chem.28795919600. 10.1074/jbc.M111.336107

  • 180

    WallaceD. C.FanW.ProcaccioV. (2010). Mitochondrial energetics and therapeutics.Annu. Rev. Pathol.5297348. 10.1146/annurev.pathol.4.110807.092314

  • 181

    WanJ.RothA. F.BaileyA. O.DavisN. G. (2007). Palmitoylated proteins: purification and identification.Nat. Protoc.215731584. 10.1038/nprot.2007.225

  • 182

    WangY.PetersonS.LoringJ. (2014). Protein post-translational modifications and regulation of pluripotency in human stem cells.Cell Res.24143160. 10.1038/cr.2013.151

  • 183

    WeckbeckerD.LongenS.RiemerJ.HerrmannJ. M. (2012). Atp23 biogenesis reveals a chaperone-like folding activity of Mia40 in the IMS of mitochondria.EMBO J.3143484358. 10.1038/emboj.2012.263

  • 184

    WeiP. C.HsiehY. H.SuM. I.JiangX.HsuP. H.LoW. T.et al (2012). Loss of the oxidative stress sensor NPGPx compromises GRP78 chaperone activity and induces systemic disease.Mol. Cell48747759. 10.1016/j.molcel.2012.10.007

  • 185

    WuC.DaiH.YanL.LiuT.CuiC.ChenT.et al (2017). Sulfonation of the resolving cysteine in human peroxiredoxin 1: a comprehensive analysis by mass spectrometry.Free Radic. Biol. Med.108785792. 10.1016/j.freeradbiomed.2017.04.341

  • 186

    XiongY.UysJ. D.TewK. D.TownsendD. M. (2011). S-Glutathionylation: from molecular mechanisms to health.Antioxid. Redox Signal.15233270. 10.1089/ars.2010.3540

  • 187

    XuX.DeckerW.SampsonM. J.CraigenW. J.ColombiniM. (1999). Mouse VDAC isoforms expressed in yeast: channel properties and their roles in mitochondrial outer membrane permeability.J. Membr. Biol.17089102. 10.1007/s002329900540

  • 188

    YamamotoT.YamadaA.WatanabeM.YoshimuraY.YamazakiN.YoshimuraY.et al (2006). VDAC1, having a shorter N-terminus than VDAC2 but showing the same migration in an SDS-polyacrylamide gel, is the predominant form expressed in mitochondria of various tissues.J. Proteome Res.533363344. 10.1021/pr060291w

  • 189

    YangW.Di VizioD.KirchnerM.SteenH.FreemanM. R. (2010). Proteome scale characterization of human S-acylated proteins in lipid raft-enriched and non-raft membranes.Mol. Cell. Proteomics95470. 10.1074/mcp.M800448-MCP200

  • 190

    YoungM. J.BayD. C.HausnerG.CourtD. A. (2007). The evolutionary history of mitochondrial porins.BMC Evol. Biol.7:31. 10.1186/1471-2148-7-31

  • 191

    ZhangD.DuJ.TangC.HuangY.JinH. (2017). H2S-induced sulfhydration: biological function and detection methodology.Front. Pharmacol.8:608. 10.3389/fphar.2017.00608

  • 192

    ZhangE.Al-AmilyI. M.MohammedS.LuanC.AsplundO.AhmedM.et al (2019). Preserving insulin secretion in diabetes by inhibiting VDAC1 overexpression and surface translocation in β cells.Cell Metab.296477. 10.1016/j.cmet.2018.09.008

  • 193

    ZhaoQ.YuJ.TanL. (2015). S-Nitrosylation in Alzheimer’s disease.Mol. Neurobiol.51268280. 10.1007/s12035-014-8672-2

  • 194

    ZhuC.TianL.YangH.ChenP.LiY.LiuY. (2019). Mitochondrial outer membrane voltage-dependent anion channel is involved in renal dysfunction in a spontaneously hypertensive rat carrying transfer RNA mutations.Eur. J. Pharm.865:172622. 10.1016/j.ejphar.2019.172622

  • 195

    ZhuangJ.WangP. Y.HuangX.ChenX.KangJ. G.HwangP. M. (2013). Mitochondrial disulfide relay mediates translocation of p53 and partitions its subcellular activity.Proc. Natl. Acad. Sci. U.S.A.1101735617361. 10.1073/pnas.1310908110

Summary

Keywords

cysteine overoxidation, outer mitochondrial membrane, voltage-dependent anion selective channel isoforms, Orbitrap Fusion Tribrid, posttranslational modification, ROS

Citation

Reina S, Pittalà MGG, Guarino F, Messina A, De Pinto V, Foti S and Saletti R (2020) Cysteine Oxidations in Mitochondrial Membrane Proteins: The Case of VDAC Isoforms in Mammals. Front. Cell Dev. Biol. 8:397. doi: 10.3389/fcell.2020.00397

Received

11 February 2020

Accepted

29 April 2020

Published

04 June 2020

Volume

8 - 2020

Edited by

Pablo M. Garcia-Roves, University of Barcelona, Spain

Reviewed by

Johannes M. Herrmann, University of Kaiserslautern, Germany; Catherine Brenner, Délégation Paris-Villejuif-01 (CNRS), France; Esther Imperlini, Institute of Research and Medical Care (IRCCS) SDN, Italy

Updates

Copyright

*Correspondence: Simona Reina,

This article was submitted to Cell Death and Survival, a section of the journal Frontiers in Cell and Developmental Biology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics